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Preface 


Preface 


“Essential Electrodynamics” and my previous book “Essential Electromagnetism” (also pub¬ 
lished by Ventus Publishing ApS) are intended to be resources for students taking electromag¬ 
netism courses while pursuing undergraduate studies in physics and engineering. Due to limited 
space available in this series, it is not possible to go into the material in great depth, so I have 
attempted to encapsulate what I consider to be the essentials. This book does not aim to replace 
existing textbooks on these topics of which there are many excellent examples, several of which 
are listed in the bibliography. Nevertheless, if appropriately supplemented, this book and my 
other book “Essential Electromagnetism” could together serve as a textbook for 2nd and 3rd 
year electromagnetism courses at Australian and British universities, or for junior/senior level 
electromagnetism courses at American universities/colleges. 

The book assumes a working knowledge of partial differential equations, vectors and vector cal¬ 
culus as would normally be acquired in mathematics courses taken by physics and engineering 
students. It also assumes knowledge of electromagnetism at the level of “Essential Electro¬ 
magnetism”, which also contains very brief introductions to vectors, vector calculus and index 
notation. Some of the mathematical derivations have been relegated to the appendices, and 
some of those are carried out using index notation, but elsewhere in the book manipulation of 
equations involving vector differential calculus is done using standard vector calculus identities 
given in the appendices. 

“Essential Electrodynamics” starts with the electromotive force and Faraday’s law, the displace¬ 
ment current, Maxwell’s equations and conservation laws. It then discusses the wave equation, 
electromagnetic waves on lossless transmission lines, in empty space, and in linear dielectrics 
(including reflection and transmission at an interface). This is followed by electromagnetic 
waves in dispersive media including dielectrics, conductors and diffuse plasmas, as well as in 
waveguides. The book ends with radiation and scattering, using first an heuristic approach to 
derive Larmor’s formula, and then apply it to simple problems before taking up a more formal 
approach using the retarded potentials in the far zone to discuss antenna radiation. 

Each chapter is followed by several exercise problems, and solutions to these problems are pub¬ 
lished separately by Ventus as “Essential Electrodynamics - Solutions”. I suggest you attempt 
these exercises before looking at the solutions. 

I hope you find this book useful. If you find typos or errors I would appreciate you letting me 
know so that I can fix them in the next edition. Suggestions for improvement are also welcome 
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- please email them to me at protheroe.essentialphysics@gmail.com. 

I am grateful to thank Professors Anita Reimer and Todor Stanev for kindly reading a draft of 
the manuscript. However, all errors are entirely due to me. This book was mainly written in 
the evenings and I would like to thank my family for their support and forbearance. 

This book is dedicated to the memory of my parents, who nurtured my interest in science. 

Raymond John Protheroe, 

School of Chemistry and Physics, The University of Adelaide, Australia 


Adelaide, May 2013 
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Electrodynamics and conservation laws 


1 Electrodynamics and conservation laws 

Learning objectives 

— To learn that a non-conservative electric field is present inside a source of an electro-motive 
force (emf) such as a battery or an electric generator, and that the total electric field in 
general has both conservative (electrostatic) Ees and non-conservative (electro-motive) 
Eem parts, such that the emf is 

£= E E m • dr = (p E E m • dr = (p (E E m + E E s) * dr. 

J inside source J via source J via source 

— To understand that a changing magnetic field produces a non-conservative electric field 
such that the emf around circuit T is minus the rate of change of magnetic flux through 
surface S bounded by the circuit, 

/ E • dr = ~[b-<IS = —/ A -dr. 

Jr dt Js dt J r 

— To understand the concept of inductance and be able to calculate the mutual and self 
inductance in simple problems. 

— To know that the change in energy density stored in a magnetic field is in general 

du — H • <£B, 

and for linear materials the energy density is 

u = B = lnH 2 = ^-B 2 . 

2 2 P 2 (i 

— To know that by introducing the concept of the displacement current, and modifying 
Ampere’s Law, Maxwell was able to unify electricity and magnetism in the form of his 
four equations which give the divergence and curl of the electric and magnetic fields. 

— To know and understand the various terms in Poynting’s theorem which expresses con¬ 
servation of energy in electrodynamics, and that the rate of energy flow is described by 
the Poynting vector S = E x H. 
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To know that the electromagnetic momentum density is g = eo'E x B = /i q£oE x H = 
S/c 2 . 


1.1 Electro-motive force 


The electro-motive force (abbreviation emf, symbol £, unit V) is what drives a current around 
a circuit. Examples of sources of emf include batteries, piezoelectric crystals, solar cells and 
electrical generators (alternators and dynamos). In electrical generators an engine fuelled, for 
example by coal, oil or nuclear fission, or a water turbine in a hydroelectric plant or a wind 
turbine in a wind farm, moves conductors through a magnetic field to produce electricity. 

The purpose of a source of emf is to maintain a potential difference across its terminals. Within 
a source of emf there is a non-conservative “electro-motive” force Fem acting on positive charges 
and pushing them towards the positive terminal “A”, and on negative charges pushing them 
towards the negative terminal “B” (see Fig. 1.1a). Integrating the electro-motive force per unit 
charge Eem = Fem /q from the negative terminal through the source of emf to the positive 
terminal gives the emf, 

£ = f F E m • dr. (1.1) 

J inside source 


(a) 


source of emf 


/ 



(b) 


magnetic field B 


y 


C 

0 


0 0 

D 

0 0 0 


0 0 

B 

0 0 

A ^ 

0 0 

0 j 0 

/ 

Ri 

0 0 

B 


Figure 1.1: (a) Source of emf: inside the source the electro-motive “force” pushes positive 
charge towards the positive terminal; outside the source the electro-motive “force” is zero and 
the electrostatic force pushes positive charge towards the negative terminal, (b) Motional emf: 
the circuit is pulled to the right at constant velocity v through a uniform magnetic field B 
(pointing into the page) which causes current I to flow. 


The electro-motive force exists only inside the source of emf. Within the source it is in the 
opposite direction to the (conservative) electrostatic force #Ees which pushes positive charges 
away from the positive terminal, and the electro-motive and the electrostatic forces cancel each 
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other out. Thus, if we integrate the vector sum of the two forces per unit charge, E tota i = 
Eem + Ees around a closed loop from the negative terminal through the source of emf to the 
positive terminal then outside the source back to the negative terminal we also get the emf 
because Eem = 0 outside the source and (p Ees • dr = 0, and so 


£ = <f E tota i • dr (1.2) 

J via source 

From now I shall represent E to tai simply by E and remember it is in general the vector sum of 
the electrostatic field and the non-conservative electric field associated with a source of emf. 


1.1.1 Motional emf 


Consider a wire moving at velocity v through a magnetic field B as shown in Fig. 1.1(b). A 
charge q in the wire experiences a force qv x B in the direction from D to C. Integrating the 
magnetic force per unit charge v x B around the circuit, we obtain the motional emf 
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£ = 


v x B • dr = vBy , 


(1.3) 


which is entirely due to section D-C as other parts of the wire either have v x B • dr = 0 
(B-D and C-A) or are outside the magnetic field region. The emf drives a current through the 
resistor which dissipates as heat at a rate £ 2 / R equal to the work done by someone pulling the 
wire through the magnetic field (assuming 100% efficiency). This is the basic principle behind 
electrical generators. 


1.1.2 Electromagnetic induction and Faraday's law 


English physicist and chemist Michael Faraday (1791-1867) discovered that the emf produced 
was proportional to the rate of change of magnetic flux through the circuit. For the circuit 
shown in Fig. 1.1(b) the magnetic flux through the circuit is 


— / B • dS — Byx. 

Js 


(1.4) 


As the loop moves, the flux decreases: 


= By^ = -Byv. 
dt y dt y 


(1.5) 


Apart from the sign (because x is decreasing) this is just the emf we found before, 


d&B 

dt 


(1.6) 


Faraday conducted experiments of the type shown in Fig. 1.1(b), as well as experiments keeping 
the circuit fixed and moving the magnet instead of the circuit, and experiments in which the 
circuit was fixed but the strength of the magnetic field was varied. In all cases a current flowed, 
and the emf was described by Faraday’s law (Eq. 1.6) which is sometimes called the universal 
flux law , i.e. a changing magnetic flux induces an emf proportional to its rate of change. The 
minus sign is there because induced currents produce a magnetic field which tends to oppose 
the change in magnetic flux ( Lenz’ law , after Russian physicist Heinrich Friedrich Emil Lenz 
1804-1865) - nature resists change! 
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It remained a puzzle why this same formula (Eq. 1.6) applied to the three quite different types 
of experiment carried out by Faraday, until Einstein developed his special theory of relativity. 


We can write Faraday’s law as 


E • dr = [ B • dS = - [ ^ • dS 

dt Js Js dt 


(1.7) 


where loop V bounds surface S. But, using Stokes’ theorem 


/.< 


E • dr = / (V x E) • dS, 


(1.8) 


we obtain Faraday’s law in differential form 


V x E = - 


<9B 

~dt' 


(1.9) 


We see that a changing magnetic field produces an induced electric field which has V x E as 
above and V • E = 0. This is in contrast to the electrostatic field for which V x E = 0 and 
V • E = p/so — note that in the static case, i.e. dH/dt = 0, we have V x E = 0 as expected. 


1.1.3 Faraday's law in terms of the vector potential 

The magnetic flux through a surface S bounded by closed curve F is 


§b = f B(r) • dS = / V x A(r) • dS = / A(r) • dr, 
j s J s j r 


where we have used Stokes’ theorem to get the final result above. Hence, 

(f E • dr — —d- [ B ■ dS — ^- (£ A • dr = — (f • dr 
Jr dtj s dtj r Jr dt 

and we can write Faraday’s law in terms of the vector potential 


( 1 . 10 ) 


( 1 . 11 ) 
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Hence, the total electric field, including both the electrostatic part and the electric field due to 
a changing magnetic field, is 

E = -W- —. (1.12) 

d t v ; 

1.1.4 Mutual inductance 


Since electric currents produce magnetic fields, time-varying currents produce time-varying 
magnetic fields. Hence, a time-varying current in a circuit will produce a time-varying mag¬ 
netic flux through a nearby circuit loop causing an emf in that loop. The mutual inductance 
was introduced in Chapter 5 of “Essential Electromagnetism” where it was defined as follows. 
Suppose you have two coils of wire at rest and you run current I\ around Coil 1 the magnetic 
flux due to Coil 1 through Coil 2 divided by I\ is the mutual inductance M of the two coils. 
For the case of time-varying currents, from Faraday’s law, the emf in Coil 2 due to a changing 
current in Coil 1, and the emf in Coil 1 due to a changing current in Coil 2 are 


£2 = 


d&B2 

dt 



£1 = 


d&Bi 

dt 



(1.13) 
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As an example of calculating the mutual inductance we consider a solenoid comprising N\ 
coils tightly wound around a cylindrical rod made of a material with magnetic permeability /i, 
and a second solenoid comprising N 2 coils loosely wound around the first solenoid as shown 
in Fig. 1.2(a). It is easier to calculate the flux due to the tightly-wound solenoid through the 
loosely-wound solenoid, than the flux due to the loosely-wound solenoid through the tightly- 
wound solenoid. The magnitude of the magnetic field inside the tightly-wound solenoid is 
approximately 




(1.14) 


and so the magnetic flux threading the N 2 coils of the loosely-wound solenoid is 


$B2(t) = N 2 AfJ, — Ii(t) 


(1.15) 


giving 


M = 


NiN 2 A 


(1.16) 



Figure 1.2: (a) Example of mutual inductance of a loosely-wound solenoid around a tightly- 
wound solenoid, (b) Circuit symbols for transformers. 

A circuit such as that shown in Fig. 1.2(a) could serve as a transformer in AC circuits as the 
output voltage would differ from the input voltage. Generally if the two coils share the same 
magnetic flux per individual coil element or “turn” the ratio of the voltage across the secondary 
coils ( N s turns) to that across the primary coils (N p turns) is 


Vs = Ns 
V p N p 


(1.17) 


In the example above the tightly-wound solenoid is the primary (p) and the loosely wound one is 
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the secondary (s). To ensure the two coils share the same magnetic flux per turn, transformers 
may have a core made of a high permeability material, i.e. one having [i /iq such as soft 
iron. Circuit symbols for air-core and iron-core transformers are shown in Fig. 1.2(b). The core 
would often be in the shape of a torus to minimise magnetic flux leakage. 


1.1.5 Self inductance 


A changing current not only induces an emf in a nearby circuit, it also induces an emf in the 
source circuit. As an example, we consider a variable voltage source V(t) connected to a solenoid 
wound around a cylindrical rod made of a material with permeability /a. as in Fig. 1.3(a). 



Figure 1.3: (a) Solenoid connected to a variable voltage source; arrow shows direction of the 
back-emf induced due to an increasing current, (b)—Circuit symbols for an inductor. 


The magnetic flux through the solenoid is proportional to the current = LI where L is 
the self-inductance (or simply inductance ) of the loop which depends only on its geometry. In 
the case of the solenoid shown, with current /(£), the magnitude of magnetic field inside the 
solenoid is approximately 

B{t) = n^I{t), (1-18) 

and so the magnetic flux threading the N coils of the solenoid is 


$s(t) = NAn—Ut) 


(1.19) 


giving 


N 2 A 

L = “— 


( 1 . 20 ) 
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If the current changes, the emf around the loop is 

a- 

Thus, changing the current through a circuit induces a “back emf” which tends to oppose the 
change in magnetic flux (Lenz’ law) as shown in Fig. 1.3(a). Circuit elements designed to 
have significant self-inductance are called inductors and would usually have a core made of a 
magnetically soft material, e.g. soft iron and, as with transformers, the core would often be in 
the shape of a torus to minimise magnetic flux leakage. Circuit symbols are shown in Fig. 1.3(b). 


B ^dl 

dt dt 


1.1.6 Energy stored in magnetic fields 

A certain amount of energy is needed to start a current in a circuit - work is done against the 
back emf to get the current going and set up the magnetic field. The work done to move unit 
charge against the back emf in one trip around the circuit is £. The amount of charge per unit 
time passing down the wire is /, so the work done per unit time is 
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dW 

dt 


-SI 



The total work done in establishing the magnetic field is 


W = L j Idl = i LI 2 . 


( 1 . 22 ) 


(1.23) 


We can show that the work W done in establishing the magnetic field is stored as magnetic 
field energy U. In establishing the field, the power source does work at rate dW/dt = dU/dt = 
—IS = — f JE, for current density J(r) = V x H(r) and induced electric field E(r) = —dA/dt. 
Then 



(1.24) 

(1.25) 

(1.26) 
(1.27) 


where the product rule identity for V • (a x b) (Eq. E.6) is used to get Eq. 1.25 and Gauss’ 
theorem to get Eq. 1.26. As the volume integrals are over all space, the surface integral is for 
the surface at r = oo, and with the integrand dropping faster than dS ~ r~ 2 we arrive at 
Eq. 1.27. This tells us that the change in the magnetic energy density due to a small change 
dB (r) in the magnetic field is 


d[u(r)] = H(r) • dB(r). 


(1.28) 


For linear materials, 


<*Kr)] = H(r)-4H(r)] = \d [H(r) • /xH(r)] = l -d [H(r) • B(r)]. 


(1.29) 
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Hence, the magnetic energy density in linear materials is 


^( r ) = = 


H(r)B(r) 


l*m r)] S 


2 fi 


(1.30) 


1.2 Maxwell's equations 


Scottish theoretical physicist James Clerk Maxwell (1831-1879) unified electricity and mag¬ 
netism into classical electrodynamics. A logical development of his equations follows. In going 
from electrostatics (V x E = 0) to electrodynamics, E is no longer conservative, and V x E 
is given instead by Faraday’s law. In a similar way, Ampere’s law needs modifying. To show 
this, consider the second-derivative V • (V x a) which must be zero for any vector field a(r). 
Substitution shows that E(r) satisfies this condition in electrodynamics, but that B(r) does not 
unless Ampere’s law is modified. 

The divergence of Ampere’s law 

V • (V x B) = V • (/i 0 J) (1.31) 

is non-zero in electrodynamics because charge conservation requires that 

vj - -% < L32 > 
Eq. 1.32 is the continuity equation. 

We can gain further insight to the problem by considering the charging a capacitor as in Fig. 1.4. 
Ampere’s law in integral form 

j) B -dr = /i() j J-dS (1.33) 

should apply for any surface S bounded by closed loop T, such as surfaces Si, S 2 or S 3 . The 
right hand side is /iqI for surfaces Si and S 3 , but is zero for surface S 2 which goes between the 
capacitor plates. 
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Figure 1.4: The problem with Ampere’s law as illustrated by the charging of a parallel plate 
capacitor. 


1.2.1 Displacement current 


We can use Gauss’ law to replace p = • E in the continuity equation (Eq. 1.32) to get 


V J + 


c>(e 0 V • E) 
dt 


= 0. 


(1.34) 
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V • (J + J d ) — 0 


(1.35) 


where 


Jd = £o 



(1.36) 


is the displacement current density , and has the same dimensions as the current density. Since 
is caused by a changing electric field, it will be non-zero between the capacitor plates in 
Fig. 1.4. 


Using Eq. 1.35, Ampere’s Law can then be modified to make V • (V x B) = 0, i.e. 

dE 

V x B = /j, 0 J + n 0 £ 0 — (1-37) 

at 

which is Ampere’s law as modified by Maxwell. From this we may conclude that currents and 
displacement currents are on an equal footing in electrodynamics, and that a changing electric 
field induces a magnetic field. 


Returning to the problem of the charging capacitor (Fig. 1.4), and assuming the electric field 
is uniform between the plates and zero elsewhere, Gauss’ law gives 


E = 


a 


A 


(1.38) 


where a is the surface charge density, Q is the charge and A the plate area. Thus between the 
plates the displacement current is 


J D 




<9E 

dt 


1_dQ _ I_ 
A dt A ’ 


.'. Id 




= I. 


(1.39) 


Hence, the total displacement current Id between the plates is identical to the current / in the 
wires charging the plates. 
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1.2.2 Maxwell's equations in vacuum and in matter 

Maxwell modified Ampere’s law in 1865 to make it apply to electrodynamics by adding the 
displacement current to the current. Replacing £oE with D and B//io with H we obtain 
versions applicable to matter. The sets of four equations of electrodynamics thus modified are 
known as Maxwell’s equations: 


Vacuum: 


Matter: 



V -E 

P_ 

£o 

V D 

= Pf 

(Gauss’ law) 

V B 

= 0 

V B 

= 0 

(no magnetic charge) 

V x E 

3B 

- ~~dt 

V x E 

3B 

~~dt 

(Faraday’s law) 

V x B 

9E 

= IM)J + //o^o-T^r 

V x H 

Q 1 ^ 

is 1^ 
+ 

>“5 

(Ampere’s law) 


(1.40) 

These four equations together with the Lorentz force law F = q{ E + v x B) encapsulate classical 
electrodynamics. 


1.3 Conservation of energy 


In 1884 English physicist John Henry Poynting (1852-1914) published his theorem, which is an 
expression of the law of conservation of energy in electrodynamics. We shall obtain Poynting’s 
theorem from Maxwell’s equations by taking the dot product of E with Ampere’s Law and the 
dot product of H with Faraday’s Law to obtain 


E-(VxH) = E-(j + ^), 

H'(VXE) = H-(-f). 

Subtracting Eq. 1.42 from Eq. 1.41 gives 


(1.41) 

(1.42) 


E • (V x H) - H • (V x E) 


E J + 



<9D 

~dt 


+ H 



(1.43) 
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Then using the product rule V • (e x h) = h • (V x e) — e • (V x h) we obtain 

E- — + H-— J. (1.44) 

Taking the volume integral of both sides, and then using Gauss’ theorem on the left hand side 
we get Poynting’s Theorem: 


- j> (ExH)-rfS 


X 


E J d 3 r + 



where dE is the surface element, and surface S bounds volume V. 


(1.45) 


Take a look at the two terms on the right hand side of Eq. 1.45: 

/ E-J d 3 r — this is the power dissipated in volume V as kinetic energy (e.g. Joule heating); 

Jv 

E • + H • ) d 3 r — this is the rate of increase of field energy in volume V. 


X 

Hence j) (E x H) • (—dE) must be the rate of flow of energy through area E into volume V 
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(into because of the minus sign). Thus we define the Poynting Vector 

SeExH (1.46) 

which gives the energy flux (W m -2 ) carried by the electromagnetic field, and points in the 
direction of energy flow. 


1.4 Conservation of momentum 


One can derive the equation describing conservation of momentum from the force on a particle 
of charge g, as given by Lorentz force equation F = (gE + qv x B). For the charges inside 
volume V this becomes 

< Pj^= f (pE + J x B)d 3 r (1.47) 

at Jy 

where P pa rt is the total momentum of the particles inside V , and we have used pv = J. Using 
all four of Maxwell’s equations, the integrand can be written 


I = (eqV • E)E + 


-(VxB)-eo^ 
Mo at 


1 


= (e 0 V • E)E + —(V x B) x B - 


Mo 


£o 


x B, 

<9(E x B) 


dt 


„ SB' 

- e 0 E x — 
at 


= (e 0 V • E)E + —(V x B) x B - x B ) _ £qE x (V x E), 


(1.48) 

(1.49) 

(1.50) 


S(g 0 E x B) 
dt 


1 


+ s 0 [(V • E)E-Ex(VxE)] + — [(V • B)B -Bx(VxB)], 


Mo 


(1.51) 


(Gauss’ law and the modified Ampere’s law were used in Eq. 1.48, and Faraday’s law was used 
in Eq. 1.50, and V • B = 0 in Eq. 1.51). 
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Re-writing Eq. 1.47 with the terms involving time-derivatives on the left hand side 


d 

dt 


P part T 


4 


[ (eoE 

Jv 


x B) drr 


1 


e 0 [(V • E)E-Ex(VxE)] + — [(V • B)B-Bx(VxB)] d 3 r (1.52) 
v l Mo 


we can identify the left hand side as the rate of change of particle and field momentum in 
volume V, and so the electromagnetic field momentum density must be 


g = £ 0 E x B = /i 0 £oE x H = -Is. 


(1.53) 


The right had side of Eq. 1.52 must equal the rate of flow of momentum into volume V. We 
will see this more easily if we are able to re-write the right hand side as a flux integral with the 
help of Gauss’ theorem, but we must first write the integrand as the divergence of some field. 
But, since the rate of increase of momentum density is a vector field having 3 dimensions, the 
quantity we must take the divergence of will have 9 dimensions. This tensor field, T^, is called 
the Maxwell stress tensor , and momentum conservation in electrodynamics is expressed by 


d 

dt 


Ppart 4“ 



x B) d 3 r 



d 3 r, 

(1.54) 

dS, 

(1.55) 


where n is the outward normal unit vector at surface S which bounds volume V. In Appendix B 
I prove that 


Tij = £q 


EiEj + — -Sij(E • E + c 2 B • B) 


(1.56) 


Summary of important concepts and equations 


Electro-motive force £ (emf) 


— A non-conservative force per charge Eem which acts inside a source of emf to maintain a 
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voltage across its terminals 

— Examples include batteries and electrical generators 

— The total electric field includes any (conservative) electrostatic field Ees such that 

£= E E m -dr = (p E E m -dr = (p (E E m + E E s) ■ dr. (1-57) 

J inside source J via source J via source 


Faraday’s law 


— The emf around circuit T is minus the rate of change of magnetic flux through surface S 
bounded by the circuit, 

(f E • dr = — -f- f B • dS 

Jr dt J s 

— The negative sign is because the emf drives a current which tries to oppose change (Lenz’s 
law) 
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— In differential form V x E 


<9B _ <9A 

~dt’ ~ ~~dt 


— Including both conservative and non-conservative fields E 


-VI/ 


dA 

dt 


Inductance 


— Mutual inductance: magnetic flux in one circuit due to current in another M — &2/I1 

— Self inductance: magnetic flux in a circuit due to current in the same circuit L = <fq//i 

— Work done to produce magnetic field W = \LI 2 

— Energy density in magnetic field u = • B = \BB = ^fiH 2 = 

Maxwell’s equations 

— Maxwell modified Ampere’s Law to make V • (V x B) = 0, 

< 9 E 

V x B = /i 0 J + m o^ow^r 
Displacement current density: 3d = £0 

Gauss’ law, the no magnetic charge law V-B = 0, Ampere’s law (modified) and Faraday’s 
law are collectively referred to as Maxwell’s equations. 

— These four equations unify electricity and magnetism and together with the Lorentz force 
law encapsulate classical electrodynamics 

Energy and momentum conservation 



— Energy conservation is expressed by Poynting’s theorem 

— The flow of electromagnetic energy is described by the Poynting vector S = ExH whose 
direction gives the direction of energy flow and whose magnitude gives the energy crossing 
unit area per unit time 
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— Electromagnetic momentum density is g = soE x B = /io^oE x H = S/c 2 

— The flow of electromagnetic momentum into closed surface S is given by 



(1.58) 


where n is the outward normal unit vector and Tij is the Maxwell stress tensor. 


Exercises on Chapter 1 

1-1 A magnetic dipole of moment m = mz is located at the origin. A thin circular conducting 
ring of radius a vibrates such that the position of its centre is r = [zq + bcos(cjt)]z with 
6 < a C zo. The plane of the ring remains parallel to the x — y plane during the vibration. 
Find the emf around the ring in the </> direction. 

1-2 A thin disc of radius a and height h contains charge +q uniformly distributed throughout 
the disc. The disc is located with its centre at the origin, and rotates about the z-axis with 
angular velocity u — uoz. 

(a) Using cylindrical coordinates but with R being the cylindrical radius to avoid confusion 
with the charge density p( r), specify the current density J (i?, </>, z) as a function of position. 
In the limit h <C a find the magnetic dipole moment. 

(b) Consider a circular loop of radius Rq around the z-axis at height zq above the disc for 
the case Rq <C a <C zq. Find the magnetic flux through the loop, and hence find the vector 
potential at the loop. 

(c) If, due to friction in the axle, the disc’s angular velocity is decreasing exponentially 

with time t as uo{t) — where to is the decay time scale, find the electric field at 

the loop at time t — 0. 

1-3 A light rigid rectangular circuit with resistance R has mass m attached to the middle 
of the lower side (width 5), and the top side is suspended horizontally using frictionless 
bearings to form a simple pendulum of length h as shown in the diagram below. In the 
absence of a magnetic field the position of the pendulum mass would be described by 
r m (t) ~ h0Q cos(cjt) x where uo = \Jg/h. A uniform magnetic field B points in the verti¬ 
cally upward direction. 
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(a) Assuming the position of the pendulum mass is still described by r m (t) « h 0$ cos (cut) x, 
what is the magnetic flux <&#(£) through the circuit, and hence the emf as a function of 
time? Take the direction around the circuit indicated by the arrow to correspond to positive 
emfs and currents. 

(b) What is the force on the lower side of the circuit due to the magnetic field? What 
is the instantaneous work done by the pendulum against this force? Compare this with 
instantaneous power dissipated in the circuit? What are the consequences of the presence 
of the magnetic field for the motion of the pendulum? 

1-4 Consider the section of a two-wire transmission line shown below. Show that the self¬ 
inductance per unit length for the case where D a is given by 
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L = 


— In —. 
7 r a 


(1.59) 



1-5 Consider a coaxial cable as an infinite cylindrical inductor and find the inductance per unit 
length. 



1-6 The diagram shows a parallel plate capacitor. Find the current /(£), and the displacement 
current density J d between moving capacitor plates, and check that the total displacement 
current is loif) = /(£). Neglect fringing effects. 



1-7 Consider a straight piece of wire radius a and length Az, along which current / is flowing. 
The potential difference between the ends is AV. Find the Poynting vector at the surface of 
the wire and use it to determine the rate at which energy flows into the wire, and compare 
the result with Joules’s law. 

1-8 A long solenoid carrying a time-dependent current I(t) is wound on a hollow cylinder whose 
axis of symmetry is the z-axis. The solenoid’s radius is a, and it has n turns per metre. 
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(a) Write down the magnetic intensity H(r, t) and magnetic field B(r, t) everywhere. What 
is the energy density in the magnetic field inside the solenoid? 

(b) Find the electric field E(r ,£) everywhere using Faraday’s law in integral form. 

(c) Find the magnetic vector potential A(r, £) everywhere. 

(d) Find the Poynting vector S(r, £) inside the cylinder, and hence the energy flux into 
a section of the cylinder of length h and the rate of increase of energy density inside the 
cylinder. Compare this with the rate of increase of magnetic field energy inside length h 
of the cylinder. 

1-9 Using the Maxwell stress tensor find the pressure exerted on a perfectly absorbing screen 
by an electromagnetic plane wave at normal incidence. 
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2 Electromagnetic waves in empty space and linear dielectrics 

Learning objectives 

— To revise solutions of the ID wave equation, and to understand the meaning of the wave¬ 
length, frequency, phase, phase constant, angular frequency and phase velocity for it’s 
sinusoidal solutions, and how they are related. 

— To know that EM waves can propagate along transmission lines, to be able to obtain the 
wave equation for a lossless transmission line of given inductance per unit length L and 
capacitance per unit length C, and to find the characteristic impedance Z — yjL/C and 
phase velocity v p = 1 /\JLC — 1/ -y/JIe. 

Id 2 / 

— To know that the 3D wave equation V 2 / = has monochromatic plane wave so- 

v z ot z 

lutions /( r, t ) = /o exp[i(k • r — ujt)\ where k is the wave vector. 

d 2 E 

— To be able to derive the wave equation for EM waves in a linear medium V 2 E — /i£—-^-=0, 

ot z 

the phase velocity v p = 1 /y/jle, the wave impedance of the medium Z = y/s/ii, the 
solution E(r, t) = Eq exp[z(k • r — ut)\ and how it relates to B and k, and find the energy 
density, Poynting vector and intensity. 

— To understand the polarisation terminology needed for reflection and refraction, including 
the plane of incidence, and to be able to derive the laws of reflection and refraction, and 
to apply the amplitude reflection and transmission coefficients, calculating reflectance and 
transmittance, Brewster’s angle and the critical angle for total internal reflection. 


2.1 The wave equation and its monochromatic plane wave solutions 


I will start with a discussion of the one-dimensional (ID) wave equation and its general solu¬ 
tion before moving on to sinusoidal waves. Next, the amplitude, wavelength, frequency, wave 
number, phase velocity and angular frequency will be defined, and the relationships between 
them given. The concept of a complex wave function with its real part representing a physical 
quantity will then be introduced, the ID wave equation will be generalised to three-dimensions 
(3D), and the monochromatic plane wave solutions to the 3D wave equation will be given in 
complex form. 

In the simplest case (Fig. 2.1), a wave pulse travels through a non-dispersive medium at constant 
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speed v and with a fixed shape. The wave function f(x,t) represents a physical quantity, e.g. 
pressure, displacement, ^-component of magnetic field, etc. 




Figure 2.1: A wave pulse travelling (a) in the +x direction, and (b) in the — x direction, at 
speed v in a non-dispersive medium. 

In the above examples the wave functions are, 

(a) /OM) = f(x-vt,0) = ij)+{x-vt), (2.1) 

(b) f(x,t) — f(x + vt, 0) = + vt), (2.2) 


where and could be any differentiable functions. A linear superposition of these two 
functions provides the general solution of the one-dimensional wave equation 

&j_ = 1 ^/ 

dx 2 v 2 dt 2 

and represent waves travelling in the -\-x and — x directions, respectively. By partial differenti¬ 
ation with respect to x and t one can easily show that f(x,t) = — vt ) satisfies the wave 

equation: 

y = V ’+(x-vt)- = ^(-v) 2 *p'l(x - vt) = iPl(x-vt). (2.4) 

One can similarly prove that ^_(x + vt) is also a solution of the wave equation. 

The sinusoidal wave f(x,t) = focos[k(x — vt) + 6] shown in Fig. 2.2 has the form ^+(x — vt) 
and so it must satisfy the wave equation. Adjacent maxima are separated by one wavelength 
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A = 2n/k where k is the wave number. At any position x the wave function f(x,t) representing 
a physical quantity oscillates at frequency v — vj A. At time t the wave function maximises at 
positions where the phase , [k(x — vt) + 5], equals 2ni r where n is an integer and 6 is the phase 
constant. 



Figure 2.2: A sinusoidal wave with wavelength A travelling in the x direction at speed v in a 
non-dispersive medium. 

It is usual to write a sinusoidal wave in terms of its wave number k = 2 tt/X and angular frequency 
UJ — Vk — 27TZ/, 


f{x, t ) = /o cos {kx — (vt + 5). 


(2.5) 


In this case, if k is positive the wave travels in the direction, and if k is negative the wave 
travels in the — x direction. Usually engineers prefer to write /(#, t) = /o cos(cjt — kx + 5) for a 
wave travelling in the direction. 

In classical physics it is often convenient for the wave function to be complex, and to write it in 
terms of complex exponentials. Then its real part represents the physical quantity. (Caution: 
this is different to the wave function T(x,t) encountered in quantum mechanics whose modulus 
squared represents probability density.) Then the physical quantity is 

f(x,t) = |/o|Re{e^-^+ 5 )} = Re {/ 0 } = |/ 0 | cos (kx - ut + 5) (2.6) 

where i = \/~T an d /o = \fo\ elS is the complex amplitude. From now on, whenever we see a 
complex quantity the physical quantity is to be understood to be its real part. 
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The wave functions represented by Eqs. 2.5 and 2.6 are monochromatic (from the ancient 
Greek for single-colour - in general single-frequency) solutions of the wave equation, unlike the 
impulse of Fig. 2.1 which Fourier analysis would show to contain a broad range of frequencies. 
Note that engineers usually prefer to use j as the symbol representing \f—l, and would write 
f(x,t) = Re for a wave travelling in the +x direction. 


In reality waves travel in three-dimensions, and generalising the one-dimensional wave equation 
(Eq. 2.3) to three dimensions gives the three-dimensional wave equation 


&l_ = v 2 f = 

dx 2 dy 2 dz 2 v 2 dt 2 ’ v 2 dt 2 

It has monochromatic plane wave solutions 


(2.7) 


/(r, £) = /o exp[i(k • r — ut)] (2.8) 

where, as before, /o is in general complex and the (real) physical quantity is understood to be 
the real part of the complex solution. Notice that kx in the ID equation is replaced by k • r 
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where k is the wave vector. Its magnitude is just the wave number k = 2tt/X as before, and the 
wave propagates in the direction of k. The field /(r,£) is identical at all points at which the 
phase is 

(k • r — cat) = 2nn + const., (integer n), (2.9) 


i.e. in planes perpendicular to k separated by integer multiples of A. Lines drawn perpendicular 
to these planes of constant phase, i.e. parallel to k, are called rays (see Fig. 2.3a); a concept 
useful in ray optics. 



Figure 2.3: Wave crests and rays of (a) a plane wave, and (b) a spherical wave. 


The speed v actually gives the speed at which the phase propagates, the phase velocity , 


v p = 


uj 

k' 


( 2 . 10 ) 


If there is dispersion, i.e. different frequencies propagate at different velocities, then k = k(u) 
and in Chapter 3 we will need to distinguish the phase velocity from the velocity at which 
wave-packets propagate, the group velocity v g . For this reason, from now on I will use v p rather 
than v to represent the wave velocity in non-dispersive media. As we shall shortly see, v p — c 
for electromagnetic waves in a vacuum. 


Another important class of solutions of the 3D wave equation are the spherical waves (see 
Fig. 2.3b) 

/(r, t) — /o r~ l exp [i(kr — ut)\ (spherical wave travelling outwards from origin) (2.11) 
/(r, t) — /or -1 exp [i(—kr — cut)} (spherical wave travelling in towards origin). (2.12) 


37 


Download free eBooks at bookboon.com 























Essential Electrodynamics 


Electromagnetic waves in empty space and linear dielectrics 


The wave intensity is the energy transported per unit time through unit area perpendicular to 
the wave direction. This is proportional to square of wave function, and for spherical waves 
is proportional to r -2 . In the case of outgoing spherical waves this is the well known “inverse 
square law” which applies equally to the sound or light intensity from point sources. Examples 
of spherical waves travelling inwards include a convex lens or concave mirror focusing light 
waves to a point, and radio waves from a distant source being focused by a radio telescope dish 
onto a small antenna at its focus. 


2.2 Transmission lines 


Transmission lines are electrical cables which consist of two parallel conductors of uniform cross 
section, usually separated by a dielectric, and which are designed to transmit radio frequency 
(McHz to ^300 MHz) signals from one place to another. Examples include twin lead (two 
parallel wires), coaxial cable and the strip line. At lower frequencies (50 Hz or 60 Hz) overhead 
power lines are used to transmit electricity from a remote generating station to a city. At 
microwave frequencies (^ 300 MHz to ~ 300 GHz) waveguides are used, for infrared frequencies 
specialised waveguides are made of layers of different dielectrics, and at optical frequencies we 
see the application of optical fibres. In this section we shall only discuss transmission lines 
which comprise two parallel conductors. 

We will consider a section of a transmission line as if it were an electrical circuit, and so we need 
to briefly review the circuit rules due to German physicist Gustav Kirchhoff (1824-1887). His 
first rule, the junction rule, expresses conservation of charge - it requires that the rate at which 
charge flows into a junction must equal the rate at which charge flows out of that junction. 
This is conveniently expressed in terms of currents as: “at any point, the sum of all currents 
entering a junction must equal the sum of all currents leaving that junction”. For the junction 
shown in Fig. 2.4(a) this implies I 2 + I 3 + I 4 — h- 


(a) (b) (b) 



Figure 2.4: Circuits to illustrate (a) Kirchhoff’s junction rule, and (b) Kirchhoff’s loop rule. 
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Kirchhoff’s 2st rule, the loop rule, expresses energy conservation, and requires for DC circuits 
that the sum of the changes in potential around any closed path of a circuit must be zero. For 
the circuit shown in Fig. 2.4(b) which, from Ohm’s law, has current / = V/(R\ + R 2 ) flowing 
around it this implies 


V ba + V cb + V ac = 0. V- v Ri- v R 2 = 0 . (2.13) 

-Ttl + -Lb 2 -TL1 + -IL2 

In AC circuits which include inductors, the electric field is not conservative and so the potential 
is undefined. However, we can still use Kirchhoff’s 2nd rule if for the change across an inductor 
we use the emf. Hence, for the AC circuit in Fig. 2.4(c) we have 

£(t) - I(t)R- L^~ = 0. (2.14) 
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2.2.1 Lossless transmission line equations 


We will only discuss lossless transmission lines, i.e. where there is no resistance present between 
or along the conductors. Then the equivalent circuit of infinitesimal length dx is shown in 
Fig. 2.5. Note that here and in the following equations L is the inductance per unit length and 
C is the capacitance per unit length. 

I(x, t) l(x+dx, t) 

B . ^ 9 ^ D 


V(x+dx, t) 


E 

dx 

Figure 2.5: Equivalent circuit of infinitesimal length dx of a lossless transmission line. 


V(x, t) 



Applying KirchhofFs loop rule to loop ABDEA, 

V(x,t) — Ldx d — V{x + dx,t) — 0, (2.15) 

and KirchhofFs junction rule to junction C gives 

I(x,t) — C dx ^ — 7(x = 0- (2.16) 

From these we obtain the lossless transmission line equations also known as the telegraphers’ 
equations 


dV(x,t) _ dl(x,t) 

dx ~ dt ’ 

dl(x,t) _ dV(x,t) 

dx ~ dt 


(2.17) 

(2.18) 


Differentiating Eq. 2.17 with respect to x and Eq. 2.18 with respect to £, 

d 2 V(x 1 t) d 2 I(x,t) d 2 I(x,t) d 2 V(x,t ) 

dx 2 dxdt ’ dtdx dt 2 
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Then, eliminating d 2 I/dxdt , we get 

d 2 V(x,t ) d 2 V(x,t ) 

•' dx 2 dt 2 

which is the wave equation for V(x,t) on a lossless transmission line. 


( 2 . 20 ) 


If instead we differentiate Eq. 2.17 with respect to t and Eq. 2.18 with respect to x, and then 
eliminate d 2 V / dxdt, we get 


d 2 V (x, t ) d 2 /(x, t) d 2 /(x, t) 

dtdx dt 2 ’ dx 2 

d 2 I(x,t ) d 2 I(x,t) _ 

9x 2 dt 2 


dxdt 


which is the wave equation for I(x,t) on a lossless transmission line. 


( 2 . 21 ) 

( 2 . 22 ) 


Equations 2.20 and 2.22 are in the form of the ID wave equation (Eq. 2.3), and so we see that 
signals must travel along transmission lines as waves with phase velocity 


Vlc' 


(2.23) 


For harmonic time-dependence with angular frequency uo we can use complex exponentials 


V(x,t) = V f e i{ - kx ~^ + v b e ^- kx ~ ut) , J(x, t) = I f e^ kx ~ ut) + , 


(2.24) 


where k = cj/v p , and Vf and V& are the complex voltage amplitudes of waves propagating in 
the forward and backwards directions, i.e. in the positive and negative x-directions (we take 
the real parts for the actual voltage and current). Similarly, If and I 5 are the complex current 
amplitudes of waves propagating in the positive and negative x-directions. 


The characteristic impedance of the transmission line is 
7 _ V, _ Vi 

“ T, - V 


(2.25) 


We can obtain the characteristic impedance by substituting one of the harmonic solutions, say 
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for the forward wave in Eq. 2.24, into one of the transmission line equations, say Eq. 2.17. Then 


-ik v f e i{ - kx ~ ut) 


-iuLI f e i{ - kx ~ ut \ 


(2.26) 


Z = 


uL t L [L 

~k ~ Vp ~ 7W Vc' 


(2.27) 


As an example, we consider two types of transmission line, the two-wire transmission line 
(Fig. 2.6a) and the coaxial cable (Fig. 2.6b). 




Figure 2.6: Geometry of (a) two-wire transmission line, (b) coaxial cable. 
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For the two-wire transmission line 


L = — In 

TV 



c = 


TVS 




and so 


Vp — 


1 



For the coaxial cable line 


L = — In 
2 n 



C 


1 *(T 


and so 


v p = 


1 



(2.28) 


(2.29) 


(2.30) 


(2.31) 


Notice that in both cases v p is the same and equal to ( fie) ~ 1 ,/2 for the insulating material having 
permeability fi and permittivity £ — this is true for any lossless transmission line — but that 
Z depends on the geometry (type and dimensions) of the transmission line. 


2.2.2 Standing waves on a lossless transmission line 

Unless the line is terminated with the characteristic impedance there will also be reflected waves 
traveling in the —x direction, and the wave amplitude will be the sum of incident and reflected 
components. If the line is open-circuit, i.e. nothing connected to output terminals, the resistance 
is infinite. In this case the reflected wave will be inverted and the result will be a standing wave. 
To see this let’s apply the boundary condition on the current, namely at the open circuited end 
the current must be zero as the charge cannot move in or out of the open-circuit end of the 
transmission line. If we set the open-circuit end of the transmission line to be at x — 0, then 
we must have 

1(0, t) = I f e~ iu>t + I b e~ iojt = 0. (2.32) 
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This condition can only be met for all t if 1^ = —//, he. total reflection with the reflected wave 
inverted. For convenience, we can set A = /q/2z, and then we have 


Akx _ p —ikx 

I(x,t) = To-—-= To sin(kx) e~ lujt . (2.33) 

Taking the real part, we have a standing wave as in Fig. 2.7, 

/(x,t) = Iq sm(kx) cos(ujt). (2.34) 



Figure 2.7: The current as a function of position along an open circuited lossless transmission 
line is a standing wave: red solid curve — current at t = nT where n is an integer and T = 2i t/uj 
is the period; blue dashed curve — current at t = (n + \)T. Note that curent nodes (where 
1 = 0) occur at every half wavelength from the open-circuit end where x = 0; the wavelength 
is A = 2ir/k. The lower plot shows direction of current flow in each conductor at t = nT. 


2.3 EM waves in vacuum and linear media 


Starting with Maxwell’s equations, we shall shortly derive the wave equations for the electric and 
magnetic fields and show that they predict the existence of electromagnetic (EM) waves in free 
space travelling at speed c. The relationships between the fields, the wave vector and angular 
frequency for a valid monochromatic EM plane wave will be determined, and the energy density, 
Poynting vector and wave intensity will be derived. Since monochromatic EM plane waves are 
an idealisation, we discuss when this may be a good approximation for quasi-monochromatic 
EM waves through introducing the concepts of coherence. 

As I shall shortly show, Maxwell’s (1865) theory predicted the existence of electromagnetic 
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waves travelling at the speed of light. Experiments by Hertz in 1886, in which he produced and 
detected radio waves, provided proof of Maxwell’s theory. At the time it was thought that a 
medium was required for propagation of EM waves. This “luminiferous aether” was supposed 
to be at rest while the Earth moved through it - spinning around its axis once a day while 
orbiting the Sun once per year, as the Sun orbits around the centre of the Galaxy, etc. 

If the aether were fixed, while the Earth moved through it, then the light transit time over a 
fixed baseline on Earth would depend on the orientation of the baseline with respect to the 
“aether wind”. In 1887 Michelson and Morley, and later experimenters, found that there was 
actually no difference with orientation, proving that there is no aether and that EM waves can 
propagate without a medium. This is unique, as all other classical waves require a medium to 
propagate in (e.g. sound waves in air). 

To show that EM waves in a vacuum are indeed allowed, we need to derive their wave equation. 
We will actually derive the wave equation for a charge-free and current-free linear medium. The 
results we obtain for this case can easily be specialised to the vacuum case by setting /i /io 
and e £q. 
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First, take the curl of Faraday’s law, giving 


Vx(VxE) = Vx(-^) = ~(VxB) 


(2.35) 


Next, using the 2nd derivative rule V x (V x a) = V(V • a) — V 2 a we get 


V(V-E)-V 2 E = 




(2.36) 


In a charge-free region V • E = p/e o = 0, and we can use use Ampere’s law to replace V x B: 


0-V 2 E = - — 


d 

dt 


T , dE 


(2.37) 


Hence, for a charge-free (p = 0) and current-free region (J = 0) 


2 <9 2 E 

V E - ^ £ ~d/2 =Q 


(Wave Eqn. for E). 


(2.38) 


We can derive the wave equation for B in a similar way by taking the curl of Ampere’s law, 
using the same 2nd derivative rule, and finally using Faraday’s law to obtain 


v2b -m4?=° 


(Wave Eqn. for B). 


(2.39) 


Equations 2.38 and 2.39 are identical in form to three-dimensional wave equation (Eq. 2.7), and 
this shows that Maxwell’s equations predict that EM waves are able to propagate through a 
linear material with phase velocity 


1 




(2.40) 


We also see that EM waves can propagate in a vacuum (/i = /io, e = £o), and the phase velocity 
of electromagnetic waves in vacuum, i.e. the speed of light, is 


1 


c = 


v / M0^0 


= 299792458 m s 


-l 


(2.41) 
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The monochromatic plane wave solutions of Eqs. 2.38 and 2.39 are 

E(r, t) — Eg exp[z(k • r — cut)\ and B(r, t) — Bq exp[i(k • r — cut)\. (2.42) 

Although these are solutions of the wave equation, they are not valid EM waves unless they 
also satisfy Maxwell’s equations. What constitutes a valid EM wave is discussed next. 


2.3.1 Relationships between k, E, B and H 


The monochromatic plane wave solutions must have Eo, Bo, w and k such that E and B are 

consistent with Maxwell’s equations. To find the relationships between k, E and B we shall 

need the time derivative, divergence and curl of E = Eoe z ( k ' r-a; ^ and B = Boe^ k ' r_a; ^, as given 
below, which are derived in Appendix C, 

zcjE, V • E = ik • E, V x E = ik x E, (2.43) 

iuj B, V • B = ik • B, V x B = ik x B. (2.44) 

Putting E = Eoe z ( k ' r-wt ) and B = Boe z ( k ' r-a;t ) into Maxwell’s equations (with p — 0 and J = 0) 
we obtain 


— E - 
dt 


at 


B = 


V-E = 0 —> k-E = 0, (2.45) 

k B = 0, (2.46) 

k x E = cjB, (2.47) 

k x B = — /io^E. (2.48) 

We see immediately that k, E and B are orthogonal vectors ask-E = 0, k-B = 0 and 
E B = E (k x E)/w = 0. Hence, the EM field of a valid monochromatic EM plane wave is 
described by 


V-B = 0 


__ OB 

v * E = -ar 
VxB = IJ£ m 


E(r, t) — Eo exp[i(k • r — ut)\, (2.49) 

B(r,t) = c<; _1 k x E(r, £), (2.50) 
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subject to k • E 0 = 0. Finally, since v p — oo/k = {/is) 1 / 2 , 
B(r,t) = ~jJEE(r,t ) = —E(r,t), 

Vp 

H(r,t) = E(r,t) = ±E(r,t), 

fi v ^ 


(2.51) 

(2.52) 


where the constant Z is called the wave impedance of the medium, and for free space takes the 
value ~ 377fi. Fig. 2.8 is a cartoon showing the E, B, k and v p vectors for a 

monochromatic EM plane wave. 



Figure 2.8: The electric and magnetic fields of a monochromatic EM plane wave travelling 
in the +z direction with the electric field in the direction and magnetic field in the =by 
direction. 
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2.3.2 Energy and momentum in electromagnetic waves 

Since for an EM wave B — ^/Jls E, the electric and magnetic contributions to energy density 
are equal, 


eE 2 _ B 2 
~ 2 ~ ~ 2 /? 


(2.53) 


and so the total energy density in the EM wave is 


u — 


eE 2 B 2 \ 2 B 2 

~ + v) ~ eE ~ -j’ 


u(r,t) = eEq cos 2 (k • r — cut + 8 ). 


(2.54) 

(2.55) 


The energy flux is given by the Poynting vector S = E x H, and for EM waves in a linear 
non-dispersive medium this is 


S = Ex 


1 

—t 

Lm 


ui x k x E 




E 2 k = 


1 




eE 2 k = v p eE 2 \a. 


S(r ,t) = v p u(r,t)]t. 


(2.56) 

(2.57) 


The intensity (W m 2 ) of the wave is the time-average of the magnitude of the Poynting vector 


I = v p e |£’o| 2 (cos 2 (k ■ r — cut + 8)) = v p 


s\E 0 \ 2 


(2.58) 


Since the momentum density is g = S/c 2 (as discussed in Chapter 1), the momentum flux of an 
EM wave is v p g. Its magnitude is just the amount of momentum crossing unit area of a surface 
perpendicular to k per unit time, and we may use it to find the radiation pressure. If a parallel 
beam of radiation is perfectly absorbed then the radiation pressure is 

v v 2 

Prad,abs = ( u ) (2.59) 

whereas if the radiation is perfectly reflected then the radiation pressure is twice as large. 

To obtain the radiation pressure of isotropic radiation we need to integrate the momentum 


49 


Download free eBooks at bookboon.com 









Essential Electrodynamics 


Electromagnetic waves in empty space and linear dielectrics 


flux incident on area A of one side of a plane surface, taking only the component of momentum 
density perpendicular to the plane, i.e. geos#, and the projected area A cos# (Fig. 2.9), 


Prad,iso 


1 

A 



cos # A cos # d cos # 



(2.60) 


Note that the radiation pressure in empty space is p = ix/3, and this is the equation of state of 
a photon gas. 



Figure 2.9: Geometry used in calculating the radiation pressure of isotropic EM radiation. 


2.4 Coherence of EM waves 


The monochromatic plane wave solutions are an idealisation because they are of infinite extent 
in all directions and are pure sinusoidal waves, i.e. the spread in angular frequencies is Aoj = 0. 
The same is true of monochromatic spherical waves. Such waves are completely coherent Two 
identical copies of such waves, say each with amplitude Eq and peak intensity Iq = (£o/2)EqC, 
can be made to interfere in a predictable way and show interference fringes alternating between 
100% constructive interference with / max = (eq/2){Eq + Eq) 2 c = 4/q and 100% destructive 
interference with / m i n = (eo/2)(Eo — Eq) 2 c — 0, such that the visibility of the fringes is 

V = — ~ — = 1. (2.61) 

^max T -^min 


Real EM waves are only partially coherent, but they may display coherence at a particular 
location over a limited time called the coherence time £ c , due to a partially coherent wave train 
of length equal to the coherence length l c — v p t c passing through that location. The width 
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of the wave front across which coherence is maintained is called the coherence width d c . The 
criterion for determining these paramaters is the time, or distance, over which the phase differs 
significantly from that of a pure sinusoidal wave. The coherence time, length or width may be 
used to determine the conditions under which a monochromatic EM plane wave may be a good 
approximation to the EM field present. 

Spectral line emission can be considered quasi-monochromatic , but even then there will be a 
narrow range of frequencies about the nominal transition frequency. Thus any source of EM 
waves will in practice produce waves of finite extent and have a spread of frequencies. For 
example, for an atomic transition the natural line width is A uo ~ T where T is the spontaneous 
transition rate (Einstein A coefficient for spontaneous emission). 

Fourier theory shows that AxAk x ~ 1. Given that v p = cj/fc, and that for a wave travelling in 
the x direction k x = fc, it follows that Ak x = Auj/v p such that the wave trains associated with 
emitted photons will have a coherence length 

l c = Ax~^ ( 2 . 62 ) 
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as illustrated in Fig. 2.10(a). As an example, consider the Sodium D lines which have A ~ 
590 nm and T ^ 6.1 x 10 7 s _1 , for which we find that 4 ~ 5 m or about 8 million wavelengths. 


(a) Vp (b) 



Figure 2.10: (a) A finite EM wave train due to an atomic transition, (b) A phase shift resulting 
from a collision during the transition. The coherence length l c is the smallest coherent part of 
the wave train. 


In practice, the coherence time will be shorter and the observed width of a spectral line broader. 
If, for example, during the transition the atom collides with another atom then a phase shift 
is introduced, and the coherence length will be reduced accordingly (Fig. 2.10b). This effect 
becomes more important with increasing density and temperature and is referred to as pressure 
broadening. Also, because of thermal motion, the frequency of a photon emitted by an individual 
atom will be doppler shifted by an amount depending on the atom’s component of velocity 
toward the observer. This Doppler broadening produces a gaussian line shape with line-width 
proportional to y/T. 


To determine the coherence width consider a source of diameter represented as an aperture 
illuminated by a quasi-monochromatic plane wave in Fig. 2.11 and a screen a distance L away. 
Consider the EM field at point O on the screen. It has contributions from the entire source 
(all points between A and B). These contributions will have a definite phase at O provided the 
difference between paths AO and CO is much less than half a wavelength, and this is equivalent 
to the condition kD 2 <C L. 


quasi-monochromatic 
light 



Figure 2.11: Geometry used in estimating the coherence width. 


As we move across the screen away from point O to point P coherence will be lost as the 
difference between paths AP and CP becomes comparable to half a wavelength (phase difference 
Ac/) = 7r). The distance r c where this occurs is the radius of a patch of the wavefront called the 
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coherence area , A c = i rr^, over which the wavefront may be considered spatially coherent (the 
coherence width is d c = 2 r c ). The phase difference is Acf) = k(L+ — L_) ~ krD/L provided 
D < r <C L. Then, 


d c 


27 tL 
~kD 


XL 

~D' 


(2.63) 


There exist both artificial and natural sources of coherent EM waves. One example of a natural 
source of coherent radiation is an astrophysical maser where a molecular line in the microwave 
region of the spectrum is observed. Maser stands for microwave amplification by stimulated 
emission of radiation. Astrophysical masers occur in shocked regions in molecular clouds in the 
Galaxy. An upper long-lived level of a molecular species is “pumped” by collisional excitation 
or infrared photons giving rise to a population inversion (more molecules in the upper level 
than in the lower level they decay to). This population inversion causes stimulated emission to 
dominate over absorption, resulting in a negative absorption coefficient. 

A single photon from spontaneous decay can initiate maser action. As it propagates through 
the region of the molecular cloud having the population inversion the photon triggers repeated 
stimulated emission of essentially identical photons in phase with itself, as do the photons 
produced by stimulated emission. There is thus a cascade or avalanche in which the number of 
photons build up exponentially with distance, and because the path-lengths through the cloud 
are huge the resulting brightness (intensity per unit solid angle) of coherent microwave radiation 
is sufficiently high to be observable with radio telescopes. 

Artificial sources of coherent EM radiation include both masers and lasers. In the laboratory, 
artificial masers using hyperfine splitting of the 21cm line spin-flip transition of the ground state 
of neutral atomic hydrogen in a weak magnetic field (Zeeman effect discussed in Chapter 6 of 
“Essential Electromagnetism”) are used as a frequency standard. 

The laser (light amplification by stimulated emission of radiation) has many practical appli¬ 
cations and works on the same principle as the astrophysical maser just discussed. The main 
differences are that transitions typically in the optical or infrared are used, and various methods 
(depending on the medium) are used to pump the upper (metastable) level. Also, in the laser, 
the enormous path lengths needed are achieved by having a laser cavity containing the laser 
medium between two parallel mirrors in which the light is reflected back and forth, with some 
small fraction escaping from one of the mirrors (which is partially transmitting) to form a laser 
beam. 
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The subject of coherence (and the closely-related field of quantum optics) is important in 
physical optics and is a field in its own right of which we have just skimmed the surface, but 
suffice it to say here that the electromagnetic field in the centre of a laser beam is an excellent 
approximation to that of a monochromatic EM plane wave, and that many of the results I 
shall derive in the remainder of this chapter for monochromatic plane waves will apply also to 
partially coherent quasi-monochromatic beams. 


2.5 Polarisation of EM waves 


In this section linear, circular and elliptical polarisations will be described in terms of the electric 
field of a monochromatic plane wave propagating in the direction which can be written 

E(r, t) = (E 0 , x e iS *5t + E 0 , y e iS y y)e i{ ~ kg ~ ut \ (2.64) 

The actual field is 
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Re{E(r,£)} = xEo,^ cos(kz — uot + 6 X ) + y^o,y cos(kz — ut + 5 y ). (2.65) 

Depending on the ratio of Eo, y to Eo :X and the relative phase (S y — S x ) the monochromatic plane 
wave will be linearly, circularly or elliptically polarised. 


2.5.1 Linear polarisation 

In linear polarisation the two phase constants are equal, 5 y = 5 X = 5, and so 

Re{E(r, t)} — (xEq^ + y^o,y) cos (kz — ut + S). (2.66) 

The locus of the tip of the E vector performs simple harmonic motion along a line at angle 
6 = arctan^o^/Eo^) to the E^-axis in the E x -E y plane as shown in Fig. 2.12. This means 
that the EM wave is linearly polarised (plane polarised) such that the electric field direction 
is at angle 6 to the x-axis. Fig. 2.8 already showed E and B for a linearly polarised wave. 
Linearly polarised quasi-monochromatic plane waves may be produced by passing unpolarised 
light through a monochromatic filter and then through a polariser such as Polaroid film, or 
better still by passing a laser beam through a polariser. 



Figure 2.12: The locus of the tip of the E vector and the vector in the E x -E y plane for 
the case of linear polarisation. 
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2.5.2 Circular and elliptical polarisation 

In both circular and elliptical polarisation the electric and magnetic field vectors rotate with 
angular frequency u. In circular polarisation Eq^ = Eq^ x = Eq and (5 y — 6 X ) = ±7r/2, so 
that the locus of the tip of the E vector undergoes circular motion with “radius” Eq in the 
E x -E y plane and angular frequency c u, as shown in the Fig. 2.13(a). This means that the EM 
wave is circularly polarised. There are two possibilities depending on the sign of (5 y — 5 X ), 
usually referred to in optics textbooks as right or left polarisation, or to use particle physics 
nomenclature negative or positive helicity, respectively. Note that we are free to define S x = 0. 

When discussing Faraday rotation in Chapter 3 it will be useful to use new complex basis 
vectors for circularly polarised waves propagating in the +z direction, 

e# = (x-iy), e L = (x + zy). (2.67) 

Note that these basis vectors are not actually unit vectors as their magnitudes are both \f2 
rather than 1 (I have deliberately omitted the hat symbol), but we shall find them useful in 
their unnormalised form in Chapter 3. Then for right circular polarisation, 


E(r, t) = Eq en e i ( kz ~ ujt ) = e 0 (x + e“ i7r / 2 y) e i{ - kz ~^ . (2.68) 

Re{E(r,£)} = Eo[5icos(kz — ut + S x )+ysin(kz — ut + 5 X )], (2.69) 

and for left circular polarisation, 

E(r,t) = E 0 e L e^ kz -^ = E 0 (x + e +i7r / 2 y) . (2.70) 

Re{E(r,t)} = Eo[Slcos(kz — cut + 6 X )—ysin(kz — cut + 6 X )]. (2-71) 


In Fig. 2.13(a) we are located at some fixed z, and view the wave coming towards us. If you bend 
your fingers and point your thumbs in the direction of the source of the wave (i.e. thumbs point 
in the —z direction), and if the fingers of your right hand curl in the direction of rotation of the 
E vector then the polarisation is “right-handed”. Fig. 2.14 shows a wave which is right-circularly 
polarised. Note that this is the usual convention used in optics. 
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Figure 2.13: Loci of the tips of the E and VpB vectors in the E x -E y plane as observed at 
fixed z for (a) circular polarisation and (b) elliptical polarisation. (The wave is propagating in 
the +z direction.) 
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Figure 2.14: Electric field of an EM wave with right circular polarisation. Note that as viewed 
looking toward the source of the wave from a fixed point the E vector rotates clockwise. 

In elliptical polarisation the tip of the E vector rotates with angular velocity uo tracing an ellipse 
in the E x -E y plane. In the most general case the major axis of the ellipse is at some arbitrary 
angle with respect to the x or y direction. However without loss of generality we could carry 
out a rotation of coordinates such that the major axis is parallel to either the x or the y axis 
such that (S y — S x ) — ±7t/ 2 as in circular polarisation. Then 

Re{E(r, t)} = x.Eq :X cos (kz — ut + S x ) + y^o,y sin(A; 2 : — ut + 5 X ) (right), (2.72) 

Re{E(r, t)} — cos (kz — ut + 5 X ) — y^o,y sin (kz — ut + S x ) (left). (2.73) 

An example with Eq^ > Eq^ x is shown in Fig. 2.13(b). 


2.6 Reflection and transmission of EM waves at an interface between linear 
media 


In this section I will start by deriving the laws of reflection and refraction. The refractive index 
of a linear medium is defined by 


_ c ck I fie 

Vp U) V ^0^0 


where the final approximation /i ~ /iq applies to non-magnetic materials. 



Since k = uo/v p for monochromatic EM plane waves 


B(r ,t) 


| \ _ pi _ 

— k x E(r,f) — —kxE(r,t) — — kxE(r,t). 


(2.75) 
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2.6.1 Laws of reflection and refraction 

Consider a monochromatic plane EM wave incident on a plane boundary between two dielectric 
materials. Part of the wave is transmitted and part is reflected. From symmetry arguments, 
the reflected and transmitted rays must be in the same plane as the incident ray and the normal 
to the boundary. This plane is called the plane of incidence. Fig. 2.15 shows the ray geometry 
in the plane of incidence. 


hi *\ T 



Figure 2.15: Incident, reflected and transmitted rays shown in the plane of incidence. 

We use subscripts i, r and t to refer to quantities related to the incident, reflected and transmit¬ 
ted waves, respectively. We also use subscript i to indicate the material in which the incident 
wave (and reflected wave) propagates, and subscript t to indicate the material in which the 
transmitted wave propagates. 

The electric and magnetic fields of the incident, reflected and transmitted waves are 


Ej = E i} o exp{i(kj • r - urf)} 

(incident wave), 

(2.76) 

E r = E r fl exp{i(k r • r — u r t)} 

(reflected wave), 

(2.77) 

E* = E tfi exp{i(k 4 • r - u t t)} 

(transmitted wave), 

(2.78) 

n ^ 

B = — k x E (for each wave' 
c 

1. 

(2.79) 


The electromagnetic field must satisfy the boundary conditions for E,D,B and H at every 
point on the interface between the two materials at all times. For the boundary conditions to 
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hold at all times , uo must be the same for all three waves. 

For the boundary conditions to be satisfied at every point , e.g. simultaneously at points “(i)” 
and “(ii)” in Fig. 2.15, 

lq • Ar = k r • Ar = lq • Ar. (2.80) 

kiSmOi = k r sm9 r = kt sin 0 t . (2.81) 


Since ki — k r because the incident and reflected waves propagate in the same medium, we have 


0 i = Q r (law of reflection). (2.82) 

Similarly, since k = uj/v p = uon/c and co is the same for all waves, 

ni sin^ = rit sin 0t (law of refraction or “Snell’s Law”). (2.83) 
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We now know what the directions of the reflected and transmitted waves are, but how much 
of the incident intensity is reflected or transmitted? The reflectance , i?, also known as the 
reflectivity, is the fraction of the incident power that is reflected. The transmittance , T, also 
known as the transmissivity, is the fraction of the incident power that is transmitted. We 
shall see that the reflectance and transmittance are different if the orientation of the plane of 
polarisation with respect to the plane of incidence, is parallel or perpendicular to the plane of 
incidence. 

In general, the incident wave can be treated as a superposition of two linearly polarised waves 
whose planes of polarisation are orthogonal, with one wave polarised parallel to, and the other 
perpendicular to, the plane of incidence. This is equivalent in quantum mechanics to the wave 
function being a superposition of eigenstates of the system. The polarisation terminology to be 
used is discussed in the next section. 


2.6.2 Polarisation terminology for reflection and refraction 


Polarisation perpendicular to the plane of incidence is written as “cr”, “s” or “TE mode”, and 
in equations as “_L”. Polarisation parallel to the plane of incidence is written as u 7r” , “p” or 
“TM mode”, and in equations as “||” This confusing nomenclature comes from the German 
and Greek, s for senkrecht (perpendicular) and p for parallel (parallel), and their Greek letter 
equivalents. TE mode refers to transverse electric, i.e. in this case E is perpendicular to the 
plane of incidence , and TM mode (transverse magnetic mode) is when B is perpendicular to the 
plane of incidence. I shall not use this nomenclature here because the incident, reflected and 
transmitted waves of both polarisations are actually all TEM waves, i.e. transverse electromag¬ 
netic waves with both E and B perpendicular (transverse to) to k. I shall instead reserve this 
nomenclature for when I discuss waveguides where there are wave field components parallel to 
k, with TE waves having no component of E parallel to k, and TM waves having no component 
of B parallel to k. 


2.7 Fresnel equations 


In the following sections we shall obtain the Fresnel equations for the amplitude reflection 
coefficients, and the reflectance and transmittance vesus the angle of incidence 6 for the two 
polarisations. 
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2.7.1 Boundary conditions at interface 


As an electromagnetic wave intersects the interface between two dielectrics, it is the boundary 
conditions on the electromagnetic field that determine the properties of the reflected and trans¬ 
mitted waves. In the dielectric where the incident wave is present, the electric and magnetic 
fields are the vector sums of those of the incident and reflected waves. The boundary conditions 
that need satisfying can be derived from Maxwell’s equations in integral form and are the same 
as were derived in “Essential Electromagnetism”. The relevant geometry is shown in Fig. 2.16 
together with the Gaussian pillbox and Amperian loops to be used. 



Figure 2.16: Gaussian pillbox and Amperian loops straddling the interface to be used in 
determining the boundary conditions on the EM wave fields. 

The Gaussian pillbox is to be used with Gauss’ law § D • dS = q/ :e nc i. and the no-magnetic- 
charge law j> B • dS = 0. The pillbox is infinitesimally thin and so we are only concerned with 
free surface charge, but in our case we have erf = 0. Hence, 


£i(^i T B r ) ' ^1 — £t^t ' 
(B i + B r ) ei = B t • ei 


(B.C. onD), (2.84) 

(B.C. onB). (2.85) 


The Amperian loops are to be used with Ampere’s law (as modified) § H • dr — f (jy + ^) • dS 
and Faraday’s law § E • dr = — f ^ • dS. The Amperian loops are infinitesimally thin such that 
the contributions of ^ and ^ are negligible, and so we are only concerned with free surface 
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currents but in our case K/ = 0. Hence, 


(B^ + B r) / f^i 

• e 2 = B t /nt • e 2 

(B.C. on H), 

(2.86) 

(Bi + B r) / l^i 

• e 3 = B t /ij, t • e 3 

(B.C. on H), 

(2.87) 

(E i + E r ) 

• e 2 = E f • e 2 

(B.C. onE), 

(2.88) 

(Ei + E r ) 

• e 3 = E t • e 3 

(B.C. onE). 

(2.89) 


For the equations involving the magnetic field we can substitute B — ( n/c)E . Thus, for 
a monochromatic plane EM wave in a non-magnetic material, for which /i ^ fi o such that 
e ~ n 2 e o, we have a set of six simultaneous equations for E r and E t to be solved in terms of 
the (known) amplitude the electric field of the incident wave, E\, its angle of incidence 6i and 
the refractive indices of the two dielectrics. For the cases of polarisation perpendicular to, or 
parallel to, the plane of incidence, three of the six equations are identity equations leaving three 
equations (one of these three being redundant) which we will solve in the following sections. 
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2.7.2 Amplitude reflection and transmission coefficients for perpendicular polari¬ 
sation 



Figure 2.17: EM field geometry for perpendicular (cr, s) polarisation. The electric field is shown 
pointing out of the page. 


The EM field geometry for perpendicular (cr, s) polarisation is shown in Fig. 2.17. In this case, 
E • ei = 0, E • S 3 = 0 and B • e 2 = 0, and so the three simultaneous equations remaining are 


{El + Er)(ni/c) cos( 7 r /2 - 0j) = E^(n t /c) cos( 7 r /2 - 9 t ) 

(B.C. on B), 

(2.90) 

{-E^ + Er) [Tii /(c/i 0 )] COS 6i = -E^[n t / (c/i 0 )] cos Q t 

(B.C. on H), 

(2.91) 

{Et + E]r) = Et 

(B.C. on E). 

(2.92) 

Solving simultaneous Eqs. 2.91 and 2.92 we obtain 



_L _ cos 9i — (nt/rii) cos 9t ± 
r cos 9i + (nt/rii) cos9 t 1 


(2.93) 

E ± _ 2 c °sOi E± 

1 cos9i +(n t /rii) cosO t *’ 


(2.94) 
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Using Snell’s law nt sin Ot = n^sin^, from which 

cos 6t — \Jl — nf /n 2 sin 2 6i (2.95) 

we arrive at the amplitude reflection coefficient for perpendicular polarisation , 


r± = 


E^r cos Qi ~ \j - sin 2 

cos 9i + f frffnf 2 - sin 2 0* 




(2.96) 


and the amplitude transmission coefficient for perpendicular polarisation , 


E t 

t± = —r 

Efl 


2 cos 6i 


(2.97) 


cos Oi + a Unt/ni) - sin 2 
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2.7.3 Amplitude reflection and transmission coefficient for parallel polarisation 



Figure 2.18: EM field geometry for parallel ( 7 r, p) polarisation. The magnetic field is shown 
pointing in to the page. 


The EM field geometry for parallel ( 7 r, p) polarisation is shown in Fig. 2.18. In this case, 
B ei = 0, B • S 3 = 0 and E • e 2 = 0, and so the three simultaneous equations remaining are 


(eiE" + EiE" ) sin 0i = e t Ef sin 0 t 

(B.C. on D), 

(2.98) 

( E i + E r) n i = E t n t 

(B.C. on H), 

(2.99) 

{—E\ + Ef,) cos Oi = ~e\ cos0 t 

(B.C. on E). 

(2.100) 

Solving simultaneous Eqs. 2.99 and 2.100 we obtain 



_ {n t /Hi) cos 0i — cos 0 t || 
r (■ nt/rii ) cos 6i + cos 0 t 1 ’ 


(2.101) 

E \\ _ 2 cos 0i E \\ 

1 ( nt/rii) cos Oi + cos 0 t 1 


(2.102) 
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Using Snell’s law to replace 9 t we arrive at the amplitude reflection coefficient for parallel 
polarisation , 


r M = 


e} ( n t /«i) 2 cos Qi - \J (n t /nj ) 2 - sin 2 0* 

-^J (rat/ Hi) 2 cos Oi + (nj/n*) 2 - sin 2 #; 


and the amplitude transmission coefficient for parallel polarisation , 


t E[ 


2 (nt/nfl) cos 6t 


^1' (rit/ni) 2 cos 6t + y (nt/ni) 2 - sin 2 0* 


(2.103) 


(2.104) 


2.7.4 Amplitude coefficients and phase shift for external reflection 


By definition, for external reflection > U{. Hence, y () 2 — sin 2 6^ in Eqs. 2.96, 2.97, 
2.103 and 2.104 is real for all 0^, and so therefore are rj_, tj_, ry and £y plotted in Fig. 2.19(a) 
for light incident from air {uf ^ 1) on to flint glass (rq « 1.62). 


Notice that for the example in Fig. 2.19(a) as the angle of incidence varies from Of = 0° to 
90° that the amplitude reflection coefficient for parallel polarisation changes smoothly from 
ry = +0.24 to —1. The angle of incidence at which the amplitude reflection coefficient for 
parallel polarisation is zero is called Brewster’s angle 0#, he. ry(0#) = 0. At Brewster’s angle 
only perpendicular polarised waves are reflected, and this is one method of producing linearly 
polarised light. We may obtain Brewster’s angle by setting r y = 0 in Eq. 2.103 to obtain 


cos 9 b = 


V 1 + KM) 2 ’ 


, a n t 
tan 9b — —* 

Ui 


(2.105) 

(2.106) 


The amplitude reflection coefficients are in general complex, but in the case of external reflection, 
ry and r± are real for all 9f. If r > 0 the phase shift is zero, but if r < 0 there is a phase shift 
on reflection by n radians (180°) since e ±Z7r = —1. This is illustrated in Fig. 2.19(b). 
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Figure 2.19: Amplitude reflection and transmission coefficients (a) and phase shift (b) for 
external reflection off flint glass. 


2.7.5 Amplitude coefficients and phase shift for internal reflection 


The amplitude coefficients for the case of internal reflection r |_, £_l, ry and £y are plotted in 
Fig. 2.20(a) for flint glass to air. There is a critical angle of incidence 6 C such that if ^ 0 C 
Snell’s law cannot be satisfied as sin Qt would need to exceed 1. For Q{ > 0 C there will be no 
transmission, t — 0, but there will be total internal reflection , |r| = 1. From Snell’s law with 
Qt = 90° the critical angle is 

= sin_i (s) ■ (2 - io7) 


For internal reflection n t < ni , and so y ( nt/ui) 2 — sin 2 6^ in Eqs. 2.96-2.104 is either real if 
0i < 6 C or imaginary if 0i > 6 C . Hence, for 0i < 6 C the amplitude reflection coefficients are real, 
and so the phase shift is (j) — 0° if r > 0 and <$> = 180° if r < 0, as was the case for external 
reflection. However, for 6i > 6 C the amplitude reflection coefficients are complex, and what is 
plotted in Fig. 2.20(b) for > 6 C is actually \r±\ and |ry|. 

Adopting the convention used in optics texts for the sign of (j) (e.g. Pedrotti & Pedrotti, 2nd 
Ed., p. 414) the amplitude reflection coefficient is defined r = |r|e - ^. The phase shift (j) may 
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then be derived by noting that for 6i > 0 C Eqs. 2.96 and 2.103 may be written in the form 


r = 


— ib \Ja 2 + b 2 e ^ 


a + ib yja 2 + b 2 e+^ 


— e “2# — e ~i4> 


(2.108) 


and so the phase shift is 


(f) = 2 tan 


-l 


(2.109) 


where b = \J sin 2 6i — (n^/n^) 2 , and for perpendicular polarisation a — cos 6i and for parallel 
polarisation a — (nt/rii) 2 cos 6i. The resulting phase shifts are plotted in Fig. 2.20(b). 
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Figure 2.20: (a) Amplitude reflection and transmission coefficients, and (b) phase shift, for 
internal reflection within flint glass. (For 6^ > 0 C in part (a) \r±\ and |r||| is plotted.) 


2.7.6 Normal incidence 

For normal incidence the plane of incidence is undefined and the results we have obtained for the 
cases of parallel and perpendicular polarisations must therefore become identical as 6i —>> 0. For 
perpendicular polarisation the amplitude coefficients (Eqs. 2.96 and 2.97) for normal incidence 
are 


r±(0°) = 


1 - nt/rij 

1 + nt/rii ’ 


*l(0°) = 


1 + nt/rii' 


( 2 . 110 ) 


whereas for parallel polarisation Eqs. 2.103 and 2.104 become 


r,,(0 o ) = 


nt/rij - 1 

nt/rii + 1 ’ 


t||(0°) = 


nt/rii + 1 ’ 


( 2 . 111 ) 


From Eqs. 2.110 and 2.111 (and 6i = 0° in Figs. 2.19 and 2.20) we see that r||(0°) = —r±( 0°), 
but for normal incidence perpendicular and parallel polarisations should be identical! 


So what is going on? To understand, we must remember that these coefficients and phase shifts 
are defined relative to the electric field directions of the incident, reflected and transmitted 
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waves given in Figs. 2.17 and 2.18. When these amplitude reflection coefficients coefficients are 
applied to the electric fields of the reflected and transmitted waves in Figs. 2.17 and 2.18 for a 
small angle of incidence the fields are as in Fig. 2.21. 




Figure 2.21: Electric field directions for near-normal incidence for (a) external reflection off 
flint glass and (b) internal reflection within flint glass. 


For external reflection (Fig. 2.21a) at normal incidence the electric field direction of the reflected 
wave is opposite to that of the incident wave, whereas for internal reflection (Fig. 2.21b) the 
electric field direction is in the same direction as that of the incident wave. Hence, for reflection 
at normal incidence there is a phase shift of 180° on reflection for external reflection, but no 
phase shift for internal reflection irrespective of polarisation. 


2.7.7 Reflectance and transmittance 


The reflectance and transmittance are defined in terms of the rate of energy flow incident on 
unit area of the interface, • (—ei) = Sj cos Qj. Thus the reflectance and transmittance are 

( 2 - 

where 



R = 


Sf cos Of 
Si cos Qj ’ 


T — 


St cos Of 
Si cos Qj ’ 


S — eE 2 v p ps (n 2 £o)E 2 (c/n) 


(2.113) 
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is the magnitude of the Poynting vector, and the cos 0 factors take account of the projected 
area of the interface as seen by the incoming and outgoing waves. Hence 


R — \r 


rh cos 0 t ^2 

Tii COS Oi 


1 -R. 


(2.114) 


For the case of internal reflection (Fig. 2.20) we see that as Oi 0 C both amplitude reflection 
coefficients obey r 1, while the amplitude transmission coefficients remain finite. However, 
because cos 0t —> 0 as Oi 0 C while cos 6^ —>► cos 0 C which is finite, from Eq.2.114 we must 
have T —>► 0 as Oi 0 C . It is straightforward to show that R + T = 1, as required by energy 
conservation. The reflectance is shown in Fig. 2.22 for the cases of (a) external reflection and 
(b) internal reflection. 


Notice that for the case of external reflection shown in Fig. 2.22(a) and grazing incidence , i.e.. 
Oi 90°, the reflectance approaches 100% (R 1). This phenomenon is used in the design 
of imaging X-ray space telescopes, and a typical design based on XMM-Newton is sketched in 
Fig.2.23 
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Figure 2.22: Reflectance for (a) external reflection off flint glass and (b) internal reflection 
within flint glass. 


2.7.8 Evanescent waves and frustrated total internal reflection 

When Oi > 6 C we have total internal reflection, but there must still be an EM field present on 
the the other side of the interface to satisfy the boundary conditions. This transmitted EM 
wave will have wave vector 

k t = (sin 9 t es + cos^ei). (2.115) 

For 6i > 6 C , Snell’s law (derived from the boundary conditions) gives an imaginary value for 
cos 9 t 


cos Qt 


y/l — sin 2 0 t 




(2.116) 


Also from Snell’s law we have sin 6^ = ( rii/rit ) sin 6^, so 


kt 


umt 

c 


/ ni . ^ 

— sin Oi e 3 

\n t 


± 



(2.117) 
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X-rays from astrophysical source 



A 


0.7 m diameter 


V 


Figure 2.23: Cross-section through an imaging X-ray telescope (not to scale). Dimensions 
given are for the XMM-Newton X-ray telescope for which the 58 concentric mirrors are made 
of nickel coated with gold and are separated by 1-5 mm. 


Hence the transmitted EM wave is 


E* = S t exp{i(k t • r - cut)}, 


= S t exp < i 


un t I rii ^ 

- — sin 0i e 3 ± 

c \ n t 


= S t exp < i 


unt f nt . . _ 

- — sin Vi] e 3 • r 

c \n t 



x exp 



cut 


(2.118) 

(2.119) 


( 2 . 120 ) 


and we see that we have a wave propagating in the e 3 direction (i.e. parallel to the surface) and 
so is a surface wave. However, the wave is decaying exponentially in the ei direction, i.e. with 
distance away from the surface of the dielectric. Such a surface wave is called an evanescent 
wave. Examining Eq. 2.120 we see that the decay distance for the electric field is ^ c/cu, or of 
the order of magnitude of the wavelength, and that the wave cannot propagate into the region 
away from the dielectric as a monochromatic plane wave. 


Evanescent waves may cross a small gap (as in quantum tunnelling) if a slab of similar dielectric 


74 


Download free eBooks at bookboon.com 

























Essential Electrodynamics 


Electromagnetic waves in empty space and linear dielectrics 


is brought within a few wavelengths of the dielectric in which the wave is incident. In that case 
part of the incident wave tunnels through the gap and propagates in the slab as a monochromatic 
plane wave in the same direction as the original incident wave, but with reduced intensity 
depending on the width of the gap. This “frustrated total internal reflection” has practical 
applications in optics with beam-splitters and optical fibre junctions. 


Summary of important concepts and equations 


ID wave equation 


dx 2 


1 

v 2 dt 2 


has general solution f(x,t) = */j+(x — vt) + vp-(x + vt) 


— moves in -\-x direction 

— moves in —x direction 

— wave speed v 
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Sinusoidal wave travelling in +x direction 

— /(x, t) — fo cos (kx — cut + 5) 

— wavelength A = 2n/k 

— wave number k = 2ty/X 

— frequency v — v/ A 

— phase [k(x — vt) + 5} 

— phase constant S 

— angular frequency uo — vk = 2i tv 

— v above is actually the phase velocity v p — cu/k 

Complex form of sinusoidal wave 

— f(x,t) = Re{/oe^-^} = |/ 0 | Re { e <(**-^+«)} = |/ 0 | cos(kx - urt + 5) 

— fo = \fo\e tS is a complex amplitude. 


3D wave equation 


V 2 / 


1 

v 2 dt 2 


— monochromatic plane wave solutions /( r, t) = fo exp[z(k • r — cut)] 

— k = k k is the wave vector 


Lossless transmission lines 

— two parallel conductors of uniform cross section separated by an insulator (permittivity 
£, permeability /i) and designed to carry radio frequency signals 

— for inductance per unit length L and capacitance per unit length C, the wave equations 
are 


d 2 V{x,t ) d 2 V(x,t) d 2 I(x,t) d 2 I(x,t) 

dx 2 dt 2 _ ’ dx 2 dt 2 
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— characteristic impedance Z — y/L/C and phase velocity v p = 1 / y/LC = 1/ y/jls 


Wave equation for EM waves in linear medium 

— V’E-^O 

— phase velocity of EM waves v p = 1 /y/JIs] c = 1/y^o^o = 299792458 m/s 

— wave impedance of the medium Z — y/e/n, Zq = y/fio/so ~ 377 Q 

— EM fields: 


E(r, t) = Eo exp[z(k • r — cj£)], B(r, t) = u 1 k x E(r, £), 

B( r, i) = —#(r, i), R(r, t) = ^ E(r, i), 

v p Z 

subject to k • E 0 = 0 

— energy density i/(r, t) = eEq cos 2 (k • r — out + S) 

— Poynting vector S = v p eE 2 \l 

— intensity / = v p e\Eq\ 2 /2 


Polarisation 

— Linear polarisation 

E = (x£o,£ + y^o , y ) cos (kz — ut + 6) 

— Circular polarisation 

E = Eq [x cos (kz — ut + 5 X ) + y sin (kz — c ut + S x )] (right) 

E = Eq [x cos (kz — ut + S x ) — y sin(/cz — uot + 5 X )] (left) 


Reflection and refraction 

— plane of incidence contains incident (i) ray, normal to interface, reflected (r) and trans¬ 
mitted (t) rays 
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— law or reflection 6i = 6 r 

— refractive index of medium containing incident wave 

— refractive index of medium containing transmitted wave nt 

— Snell’s law sin 9i = nt sin 0 t 

— amplitude reflection and transmission coefficients depend on polarisation 

— parallel polarisation has E parallel to plane of incidence 

— perpendicular polarisation has E perpendicular to plane of incidence 

Fresnel equations 

— amplitude reflection and transmission coefficients 

E ~ L cos0; ~ V ( nt/rii ) 2 - sin 2 0; 

r± = ttl = - / — 

* COS0J + y ( nt/rii ) 2 - sin 2 0j 
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E t 

t± = -E = 

E/r 


2 cos Oi 


2 sin 2 Oi 


r ii = 


Er 


E'l 


E' 


cos 9i + \J {nt/ni)' 

\i) 2 cos Oi - (nt/ni)' 


(nt/ni 


2 — sin 2 Oi 


(n t /n. 


if cos Oi + yj ( nt/rii ) 2 - sin 2 Oi 


E \ 


2 (nt/rij) cos Oj 


( nt/rii ) 2 cos Oi + yj(n t /rii) 2 - sin 2 Oi 


— reflectance and transmittance R= |r| 2 and T = 1 — R 

— Brewster’s angle t\\(0b) — 0: only _L polarisation reflected, tan Ob — Tit/rii 

— total internal reflection for Oi > 0 C where sin 0 C = nt/rii 


Normal incidence 

— for external reflection the electric field direction of the reflected wave is opposite to the 
electric field direction of the incident wave (180° phase shift) 

— for internal reflection the direction of the electric field of the reflected wave is the same 
as that of the electric field of the incident wave (no phase shift) 


Exercises on Chapter 2 

2-1 Prove that the spherical waves given by 

/(r, t) = f 0 r -1 exp [i(±kr - ut)\ (2.121) 

are solutions of the 3D wave equation. 

2-2 A monochromatic plane wave E(r, t) = Eq exp[z(k • r — cot)] is travelling in the +z direction 
through a lossless linear medium with relative permittivity e r = 4 and relative permeability 
Hr ~ 1 and is polarised in the x direction. The frequency is is = 1 GHz and E has a 
maximum value of +10 -3 V/m at t — 5 ns and z = 1 m. 

(a) Find the angular frequency, phase velocity, wavenumber, wave vector, and wavelength. 
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(b) Obtain the instantaneous expression for E(r , £) valid for any position and time. 

(c) Obtain the instantaneous expression for H(r,£) valid for any position and time. 

(d) Find the Poynting vector and its time-averaged value. 

(e) Find the locations where E x is maximum when t — Os. 

2-3 Given the electric fields for the following polarisations, 

E(r, t) — (xEq,^ + yEo cos {kz — ujt + 5) (linear), (2.122) 

E(r,£) = Eo[5icos(kz — cut + S x ) — ysin (kz — uot + 5^)] (left circular), (2.123) 

E(r, t ) = cos (kz — cot + S x ) + y^o,y sin {kz — uot + S x ) (right elliptical) 

(2.124) 

with i?o,y > Eq^ x , find the instantaneous and time-averaged energy densities and Poynting 
vectors. In each case, assume the wave is propagating in a medium with permittivity e and 
permeability p. 

2-4 For a monochromatic EM plane wave incident on a plane interface between two dielectrics 
describe what is meant by perpendicular (a, s ) and parallel (7r, p ) polarisation states. 
Include a suitable diagram in your answer. 

2-5 Define what is meant by the amplitude reflection and transmission coefficients, and by the 
reflectance and transmittance. 

2-6 (a) Use the amplitude reflection coefficient for parallel polarisation 
El ( nt/rii ) 2 cos (9* - Jnt/m ) 2 - sin 2 

m(0O = nr = - / =• (2J25) 

E i (jit/nif cos 9i + y (n t / riif - sin 2 
to show that Brewster’s angle is given by both 

cos 6 b — , 1 - and tan^ —. (2.126) 

Vi + (n t /n ? ) 2 "i 

(b) Show that 

tan^ = — (2.127) 

Ui 

can also be derived from Snell’s law and by requiring the angle between the reflected and 
transmitted rays to be 7r/2. 
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2-7 Give a formula with explanations for the reflectance and transmittance in terms of the 
magnitude of the Poynting vectors of the various waves, and their angles with respect to 
the normal to the interface. How is the reflectance related to the amplitude reflection 
coefficient? How is the transmittance related to the reflectance? Why? 

2-8 Explain the physical meaning of the critical angle, and derive its formula. 
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3 Electromagnetic waves in dispersive media 

Learning objectives 

— To understand that if the wavenumber is complex the imaginary part leads to absorption, 


such that the intensity is I — I$e akr , with absorption coefficient a = 2 x Im{k(uS)}. 


— To learn that in the Lorentz oscillator model of a dielectric, for j oscillators with natural 
frequencies ujj and damping coefficients 7 j, and on average fj electrons with natural 
frequency uoj per atom, and N atoms per unit volume, the refractive index is 



— To know that the absorption coefficient near a resonance has a Lorentzian profile 



and that if (dRe{n}/dcu) < 0 (anomalous dispersion) the group velocity can exceed c, 
but that anomalous dispersion is accompanied by absorption such that causality is not 
violated. 


To be able to derive the wave equation of a good conductor, use it to obtain the dispersion 
relation, skin depth S(uj) and absorption coefficient a(uj) = 2/5(uj). 


To be able to derive the wave equation for a dilute plasma, the dispersion relation, plasma 


frequency, show that EM waves cannot propagate below the plasma frequency, and know 
that at frequencies above the highest resonant frequency a dielectric behaves like a dilute 
plasma. 

— to learn that in a magnetised plasma or dielectric, left and right circularly polarised waves 
travel at different phase velocities causing the plane of polarisation to rotate (Faraday 
rotation). 
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3.1 Dispersion and absorption 


In general, the refractive index n(uS) is complex, as is the wave number k(uS) — n{uo)/c. Hence, 
for a monochromatic plane wave in matter 


E(r, t) — Eq exp [i(k • r — ut)\ 


— Eq exp 


Re{/c}k • r — out^j exp Im{A:}k • r 


(3.1) 

(3.2) 


This means that the intensity /(r) = (S(r,t)) oc ([E(r,t)] 2 ) decays exponentially with distance, 
/(r) = Eq exp — a(o;)k • r , where a(uj) = 2 x lm{k(uj)} is the absorption coefficient. 


3.1.1 Lorentz oscillator model of dielectrics 


The molecular polarisability of atoms can be approximated by modelling the atom semi- 
classically as a nucleus of effective charge +e surrounded by a spherical electron “cloud” of 
radius ao ~ 10 -10 m corresponding to a single valence electron (charge -e). The justification 
being that the nucleus is screened by all but one of the electrons, and therefore has effective 
charge +e. 


When an electric field is applied, the electron cloud becomes displaced by r e with respect to the 
nucleus and experiences a restoring force — ftr e , where n ~ e 2 /47T£oaQ is the “spring constant” 
(see Fig. 3.1, and Chapter 4 of “Essential Electromagnetism”). The net force on the valence 
electron is 


m e 


d 2 r e 
dt 2 


—eE — 


(3.3) 


and in electrostatics, we set the net force to zero to find equilibrium displacement r e = — eE/tt, 
dipole moment p = e 2 E/^ and molecular polarisability <a po i = e 2 /n — 


In electrodynamics, E = E (t) oscillates, and the electrons oscillate in response. Because of 
various drag forces analogous to friction, we must add a “damping term” which is proportional 
to velocity, 


m e 


d 2 r e 
dt 2 


-eE (t) 


KT e - TUeJ- 


dXe 

dt 


(3.4) 
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No applied field Applied electric field 



Figure 3.1: Semi-classical model of an atom, without and with an applied electric field (dis¬ 
placement greatly exaggerated). 


The resonant frequency of simple harmonic oscillation for the valence electrons in our crude 
model of atoms is 


UJQ 



4:7T£oaQm e 


~ 2 x 10 16 rad s 1 


(3.5) 
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which corresponds to frequencies in the ultraviolet. If a monochromatic plane wave (photons) 
with frequency close to c^o were incident on the atom, the work done against the friction-type 
forces to allow the atom oscillate would cause strong absorption of photons. In reality an atom 
is not a classical system with only one resonance, but rather a quantum system with many 
resonant frequencies. Nevertheless we can use this semi-classical model to describe the essential 
dielectric properties of materials. 


The equation of motion of the electron is 


rn. 


d 2 o d 


d 2 v e o dr e 

= -eE (t) - m e uj 0 r e - m e 7 —, 


dt 2 


r e (t) = -eE (t). 


dt 


(3.6) 

(3.7) 


For the sinusoidal electric field of an EM wave E(f) oc exp (—iut), we would expect that 


r e(t) = r e;0 exp (-iut), 
(d/dt) r e = -iur e (t), 

(d 2 / dt 2 )r e = -uj 2 t e (t). 
m e [-LU 2 + Wq — i'ycu] r e (t) = —eE(t). 


Hence, the dipole moment is 


P(*) 


-e r e (t) = 


(e 2 /m e 


Wq — uj 2 — i'yu 


E(t). 


(3.8) 

(3.9) 

(3.10) 

(3.11) 


(3.12) 


For rij oscillators with natural frequencies loj and damping coefficients 7 j, and on average fj 
electrons with natural frequency ujj per atom, and N atoms per unit volume, the polarisation 
field is 


PM) 


Ne 2 

m e 



UJ- 


fj 

— uj 2 — i^jUJ 


EM). 


(3.13) 
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But, P = Xe^oE and e — (1 + Xe)^o 5 and the refractive index must be 

(3. 

We see that the refractive index is complex, and so attenuation of the wave will occur although, 
as we shall shortly show, this is only important close to resonant frequencies uoj. 


n(uj) 


eV' 2 

£o) 


= 1 + 


Ne 2 


s 0 m e ' w? 
3 =i 3 


E 


f, 


1/2 






For gases the 2nd term in brackets is small, and we can approximate the square root by the 1st 
two terms of the binomial expansion 


n(uo) 


i + ^E 


fj 


2som e ^ uj 2 — uj 2 — ’ 


Ne 2 ^ (uj ~ w 2 ) + i'YjU 
!+-- 

3 = 


2s 0 m e ^ J 3 (uj 2 - uj 2 ) 2 + 7 2 (jj 2 


nH - |l + ^rE/: 


(«? - ^ 2 ) 


2e 0 m e ^ (w 2 - uj2 ) 2 + 7 2 ^ 2 


+ 


(3.15) 


(3.16) 


(Er f li _ 

2 £o m e ^ / Vf - W 2 ) 2 + 7 7 2 w 2 ' 

J = 1 ^ J 

(3.17) 


Purely for illustration of the important features, I plot in Fig. 3.2 the refractive index for the 
rather unrealistic case where there are just two rather broad resonances. 


The wave number k — oon/c is complex, and the equation for the real part 

(3. 

is the the dispersion relation for the substance, and is used to find the phase velocity at a given 
frequency, v p (uS) = uj/k r (uS). 


k r (uo) 



Ne 2 




f: 


("1 - ^ 2 ) 


2sq m e ^ J 3 (uj 2 - uj 2 ) 2 + 7 2 lu 2 


The absorption coefficient a(u), which has units [m - 1 ] 


, is obtained from the imaginary part of 
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the wave number 


a(u) — 


2 ki(u) 

(3.19) 

2 n o 2 

N V f ■ 

£0^eC “ 3 (uJj ~ CO 2 ) 2 + 7?W 2 ’ 

(3.20) 


Provided the resonances are narrow such that (ujj + u) « 2cu, we may make the approximation 


2 'v 


a(w) * N i^rX IiL{ - u) 


3 = 1 


(3.21) 


where 

1 t7j/2) 

' ’ * ( Wj - w) 2 + (7,/2) 2 

is the Lorentzian function plotted in Fig. 3.3. 


(3.22) 


> Apply now 
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Figure 3.2: Refractive index and reciprocal group velocity in the Lorentz oscillator model for 
two resonances: ui = 3, (N e 2 / £om e ) f\ = 0 . 2 , 71 = 0 . 12 , and CU 2 = 5, (Ne 2 /som e )f2 = 
2 , 72 = 0.4. The yellow shaded regions show the frequency ranges where there is anomalous 
dispersion. The black curve shows c/v g where v g is the group velocity which I shall discuss later 
in this chapter. 


The Lorentzian function has full width at half maximum (FWHM) equal to 7 j and is normalised 
to unit area, and describes the shape of absorption and emission lines, and is found in many 
physical situations, both quantum and classical, involving driven resonant systems. Heavily 
damped oscillators have broad line-widths — respond to a wider range of driving frequencies 
around a resonant frequency 07 . The spectral line-width depends on the Q factor, which is a 
measure of the sharpness of the resonance and depends on 7 ^. 


If we wish to apply this model to real materials, we need to replace the classical oscillators with 
resonance frequencies Uj with quantum oscillators corresponding to quantum transitions from 
level i to level k with absorption of a photon with energy close to hojik = {E^ — E {), and with 
oscillator strengths f ^ and line-widths 7 ^ determined using quantum mechanics. 


88 


Download free eBooks at bookboon.com 












Essential Electrodynamics 


Electromagnetic waves in dispersive media 



Figure 3.3: The Lorentzian function. 

3.2 Dispersion 


For EM waves in a medium the phase velocity is 

/ N = ^ = c 

Vp U R e{k(uj)} R e{n(uj)} 

Since the refractive index depends on frequency, dispersion will occur. Dispersion causes dif¬ 
ferent colours present in white light to be separated when passed through a glass prism, and 
wave packets to travel at a velocity called the group velocity. 


(3.23) 


3.2.1 Wave packets and group velocity 

A wave packet is the superposition of a group of waves travelling together and localised in 
space at some initial time. For a wave packet travelling in the z direction with its electric field 
polarised in the y direction, say, and with E y (z, 0) being the ^/-component of electric field at 
t — 0, there would be a distribution of wavenumbers given by the Fourier transform 

E y (k) =[ Ey{z,0)e~ ikz dz. (3.24) 
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We can recover the original electric field by taking the inverse Fourier transform 

E y (z, 0) = lj E y (k)e ikz dk. (3.25) 

E y (z, 0) and E y {k ) constitute Fourier transform pairs, E y (z, 0) E y {k). The electric field of 
the wave packet at a later time would be 

E y (z,t) = T J E y {k)S kz -“^dk. (3.26) 


We shall now derive the velocity at which the wave packets travel. For a narrow distribution 
of wavenumbers centred on fco = k(uJo) we can perform a Taylor series expansion of u(k) about 
W 0 = w(fc 0 ), 


w(fc) = w 0 + (fc-fco)^ 

ah 


+ 


(fc — fco) 2 d 2 uj 




2! dk 2 


+ ..., 


CJO 


and substitute just the first two terms above into Eq. 3.26 to obtain 


(3.27) 
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1 f ~ / X 

( duo 

duo 

\1 

~ 2 ^ / Ey(k) exp 

i kz — uoot — k —- 
\ dk 

u t+h ik 

^0 / J 


dk , 


1 

. / _ duo 

\ 1 

« — exp 

2n 

'\ k °dk 

-w 0 M 
^0 / 


E y (k )exp 


_ ( duo 

\1 

ik z -— 

t 1 

dk 

/ _ 


(3.28) 

dk. (3.29) 


Taking the absolute values of both sides, 
\Ey(z,t)\ 

Hence, 


1 

r ~ 

/ duo 

y 


27r 

/ E y (k) exp 


q 

^0 / J 

dk 


(3.30) 


\Ey(Z,t)\ = 




£, 0 


CJO 


(3.31) 


which is in the form of a wave travelling with speed 


duo 


^ = dk 


(3.32) 


CJO 


which is the group velocity at angular frequency uoq. 


The classic example is the gaussian wave packet which is obtained by multiplying a monochro¬ 
matic plane wave by a gaussian envelope function 


E y (z, 0) = e 


— p-Z 2 K 2 (j2 z) e ~ ik 0 Z 


(3.33) 


as in Fig. 3.4(a). Generally, for a wave packet spread out by Az in z, the spread in wavenumbers 
is A k > 1/Az. The Fourier transform of a gaussian is itself a gaussian 

e- 2 /^ 2 ) <=► V^a z e-^l 2 (3.34) 

and for a monochromatic plane wave modulated by a gaussian envelope 

e -z 2 /(2a 2 z ) e -iko z ^ (T z e - a ^ k ~ k °) 2 / 2 . ( 3 . 35 ) 
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Hence, the electric field at an arbitrary time t > 0 is given by 

E y (z,t) = j ^ e -^(fc-fco) 2 /2 e iN-u;(fc)i] dk. (3.36) 

Provided Az 1/ko the wave packet propagates at the group velocity with an envelope of 
constant shape as illustrated by Figs. 3.4(b) and (c) for the example of a gaussian wave packet. 



Figure 3.4: (a) Wave packet constructed by multiplying a monochromatic plane wave at t = 0 
(black cosine curve) by a gaussian envelope function (dotted red curve), (b) The wave packet 
at t = 0. (c) The wave packet at t > 0. The black line joining black dots connects positions of 
constant phase at the two times, while the red line joining hollow red dots connects positions 
of maximum amplitude where the envelope function peaks; in this example v v > v g . 


3.2.2 Normal and anomalous dispersion 

From looking at Figs. 3.4(b) and (c) it is clear that energy is transported with the wave packet 
at the group velocity, and provided v g < c this satisfies special relativity and is termed normal 
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dispersion. In terms of the refractive index, 


Vg(Uj) = 



-l 


d fuo~Re{n(uo)} 
duo V c 



v g (uo) — c ^Re{n(<x;)} + uo 


d~Re{n} 

duo 


-l 


(3.37) 

(3.38) 

(3.39) 


The speed of light divided by the group velocity has been added to Fig. 3.2 and we see that 
provided (dRe{n}/duo) > 0 there is normal dispersion, but for (dKe{n}/duo) < 0 the group 
velocity can exceed c and this is termed anomalous dispersion. Anomalous dispersion is always 
accompanied by absorption, and so the group velocity is then no longer a meaningful concept 
- in this case it is not the speed of transportation of energy or information - and so special 
relativity is not violated. The regions of anomalous dispersion are shown by the yellow shaded 
regions in Fig. 3.2. 
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3.3 Refractive index of a conductor 


For EM wave propagation in a conductor we must include the current density J = a E, where 
a is the conductivity, when deriving the wave equation from Maxwell’s equations in matter. 
Taking the curl of Faraday’s law we derive the wave equation for the electric field 


V x (V x E) = V x (-dB/dt) 


V(V • E) - V 2 E = -|-(V x B 


9 9 

0 - V 2 E = - — 


dt 

d 

dt 


fdEX 


„ 2t , <9E <9 2 E . 

VE ' fUT w ~ IJ£ w =0 - 


(3.40) 

(3-41) 

(3.42) 

(3.43) 


Similarly, from Ampere’s law we can derive the wave equation for the magnetic field 


„ 2t , <9B <9 2 B „ 

vb _^__, i£ _ = 0 , 


(3.44) 


which we see has identical form. Substituting the monochromatic plane wave solution E = 
E 0 exp{z(k • r — ut)} into the wave equation we obtain 


[— k 2 + tiaiuj + fieuj 2 ]Ej — 0, 


(3.45) 


and hence the dispersion relation for EM waves in a conductor 
k(u) = (jiaiuj + jisuo 2 ) 1 / 2 . 


(3.46) 


For a good conductor a euj, and so we have 


k(cu) ~ i 1 / 2 (/icra;) 1//2 (3.47) 

k(u) « (1 + i)2~ 1 ^ 2 (/j J a(jj) 1 ^ 2 (3.48) 


where we have used z 1 / 2 = (e Z 7 r / 2 )V 2 = e Z7r / 4 = (cos | + i sin |) = (1 + z)2 1 / 2 . 


94 


Download free eBooks at bookboon.com 






Essential Electrodynamics 


Electromagnetic waves in dispersive media 


The wave number is complex and so the amplitude will decrease exponentially with distance 


(E(r,t)) oc exp Im{/c(c<;)}k • oc exp k • r/5(cc;)^ 


(3.49) 


where 


5(oj) = 



(3.50) 


is called the skin depth and is plotted against frequency v — oj/(2n) in Fig. 3.5 for soft iron, 
silver, lead and sea water. As an EM wave propagates through a conductor, intensity decreases 
as 


/ = To exp 2k • r/5(cc;)^ 
so that the absorption coefficient a(uo) = 2/5 (cj). 


(3.51) 



Figure 3.5: Skin depth of soft iron (cr = 11.2 x 10 6 m _1 , fi r = 50,000), silver (cr = 

62.1 x 10 6 m -1 ), lead (cr = 4.7 x 10 6 m -1 ), sea water at 20°C with 3% salinity 

(cr = 4.1 m -1 ). 

As we have seen, at high frequencies EM waves cannot penetrate far inside a conductor. One 
example is sea water which has a skin depth of about 30 cm at 1 MHz and 10 m at 1 kHz, 
depending on temperature and salinity (see upper curve in Fig. 3.5). This causes radio com- 
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munication problems with submarines, and so VLF (3-30 kHz) or even ELF (300 Hz-3 kHz) 
frequencies are used instead. 

A related phenomenon concerns the transmission of high frequency oscillating currents in con¬ 
ductors. These currents generate oscillating magnetic fields which induce electric fields, which 
in turn produce eddy currents and associated resistive power loss. These induced electric fields, 
which oppose the change in the magnetic field, are strongest in the centre of the conductor and 
effectively increase the resistivity there, with the upshot that the conduction current becomes 
confined to a thin outer layer of thickness comparable to the skin depth. Since the inside of 
a conducting wire is not used, a hollow wire works just as well and is lighter. Alternatively, 
plating a copper wire with a thin layer of silver (a better conductor) reduces power loss. Finally, 
EM waves (and associated energy) can be transmitted inside waveguides which, at microwave 
frequencies, are suitably designed hollow metal conducting pipes. 

The reflectance of good conductors is high, and so a thin coating ^0.1 /im of silver or aluminium 
is deposited on the front surface of a concave parabolic or spherical glass surface in a reflecting 
telescope. The refractive index of a conductor is complex. However, we can still use the 
amplitude reflection coefficient for external reflection (air/metal) at normal incidence (for either 
polarisation), e.g. from Eq. 2.110, to obtain the reflectance (Exercise 3-1) 
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( D) = (1 - Re{n}) 2 + (Im{n}) 2 
{ ’ (l + Re{n}) 2 + (Im{n}) 2 ‘ 


(3.52) 


Data on R e{n}(uo) and I m{n\(uj) are available and have been used together with Eq. 3.52 to 
calculate the reflectance as a function of wavelength for aluminium, silver, gold and copper 
shown in Fig. 3.6. 


We are also able to estimate the refractive index of a good conductor from its conductivity 
using the formula we derived for the wave number (Eq. 3.48), and then use Eq. 3.52 to obtain 
an estimate of the reflectance, which has been added to Fig. 3.6 (dashed curves). As can be 
seen, clearly reality is much more complicated than our simple classical theory, and requires 
condensed matter theory for a proper understanding. Nevertheless, there is good agreement 
above 2/im (2000 nm) in copper, and the agreement is acceptable above 10/im (10,000 nm) in 
aluminium, gold and silver. 


3.4 Wave propagation in dilute plasmas 


In this section we consider wave propagation in an ionised plasma comprising n e electrons and 
rii ions per unit volume. In a dilute plasma we can ignore the (typically distant) ions as they 
provide little damping to the motion of the electrons, and because of their much higher mass 
they contribute little to the current density. 


We start by deriving the electrical conductivity from the motion of individual electrons which 
undergo acceleration in response to the electric field of the wave 

m e ^ = —eE(t). (3.53) 

The current density arises from the bulk flow of electrons in the plasma 


J (t) = n e (-e)v(t), 


and so 


dJ_ 

dt 


n e e 2 


m e 


E (t). 


(3.54) 


(3.55) 
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log(/l///m) log(4//*m) 




log(A///m) 


log(4//*m) 


Figure 3.6: Reflectances of aluminium, silver, gold and copper obtained using data 
on the complex refractive index (Palik 1985, Rakic 1995, Refractive Index Database 
http://refractiveindex.info) and Eq. 3.52 are shown by the solid curves. The dashed 
curves show the result for a good conductor using Eq. 3.48 and n = ck/uj in Eq. 3.52 
and resistivities at 0°C of aluminium, silver, gold and copper p = 1/a = 2.42 x 
10 -8 , 1.47 x 10 —8 , 2.05 x 10 -8 and 1.54 x 10 -8 ft m, respectively (Kaye & Laby 
http://www.kayelaby.npl.co.uk/general_physics/2_6/2_6_l.html). 
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From this expression we can work out the conductivity for 

an oscillating electric field Eoe lujt 

J = neg2 E 0 [ e~ iujt dt , 

(3.56) 

m e J 

n e e 2 1 iojt 

= hiQ . e , 

(3.57) 

m e —lu 

, ,n e e 2 „ 

J = i E. 

m e uj 

(3.58) 


We see that the conductivity is purely imaginary, 
n e e 2 

<7 = z- 5 —, 

m e (jj 

so that the current density and electric field are 90° out of phase. Consequently, the time- 
averaged power will be zero. This is analogous to the current and voltage in an LC circuit 
where energy is continually being exchanged between electric field energy of the capacitor and 
magnetic field energy of the inductor. 


(3.59) 
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3.4.1 EM waves in dilute plasmas 


In a dilute (low density) neutral (p = 0) plasma we may assume p — po and s = £o, and obtain 
the wave equation for E in the usual way from the curl of Faraday’s law, and substituting for 
dJ/dt using Eq. 3.55, 


V(V • E) - V 2 E = - — (V x B), 


dt 

0 -v 2 e = -| 

dt 


T / 3E 

MoJ + Mo^o ( 


2tti _ ,, / n e e- \ , „ / d 2 E\ 


V 2 E — po 


n P e 


m P 


E 


V E = po 
1 - 0 


m P 


E + Mo^o 


dt 2 J 


c 2 dt 2 


(3.60) 

(3.61) 

(3.62) 

(3.63) 


Substituting the electric field of a monochromatic EM wave E = Eq exp{i(k • r — cot)} into the 
wave equation we obtain 


-k 2 - 


n e e 2 (j 2 


+ 


c 2 £om e c 2 


E = 0, 


(3.64) 


and hence the dispersion relation 


? 2 2 2 2 
fee = UJ — UJ* 


(3.65) 


where 


(jJnr) - 


n e e " 
£: 0 m e 


2 \ 1/2 


(3.66) 


is the plasma frequency. EM waves cannot propagate if uo < uo p as k would be purely imaginary. 


100 


Download free eBooks at bookboon.com 












Essential Electrodynamics 


Electromagnetic waves in dispersive media 


3.4.2 Phase and group velocity in dilute plasmas 

From the dispersion relation (Eq. 3.65) we can derive the refractive index, phase and group 
velocities 


n(cj ) 

kc 

= — 

UJ 

v p (u) 

uj 

k 

Vg(Uj) 

~dk 

du) 






i -i 


1/2 


- 1/2 


or 


1/2 


- 


(3.67) 

(3.68) 

(3.69) 


Note that for the case of a dilute plasma v p v g — c 2 , and so 


c v Q 


n{(jj) = — —. 


Vr> C 


(3.70) 


The phase and group velocities, and refractive index are plotted in Fig. 3.7. 



Figure 3.7: Phase velocity, group velocity and refractive index of a dilute plasma. 
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No EM wave propagation can occur below the plasma frequency which is the cut-off frequency 
of a dilute plasma. Above the cut-off frequency the phase velocity exceeds c but the group 
velocity is less than c. 


3.4.3 Comparison of plasmas and dielectrics at high frequency 


For frequencies well above the plasma frequency uo p the refractive index of a dilute plasma 
(Eq. 3.67) becomes 


n(uj) 



(3.71) 


We shall show that the high frequency behaviour of the refractive index of a dielectric has the 
same form, and so may be approximated by the refractive index of a dilute plasma. 
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For frequencies well above the highest resonant frequency, i.e. uo ujj for all j, we can neglect 
absorption (set 7 j = 0). Then the refractive index of a dielectric in the Lorentz oscillator model 
(Eq. 3.17) becomes 


n{uj) ~ 



Ne 2 ^ (^-^ 2 ) ) 

2 som e (w 2 - w 2 ) 2 + 7 2 w 2 y 


2som e ^ J uj 2 

j=i 


n(u;) « 1 


2 w 2 


where 


(N'e 2 \ l/2 
Up ~ \£om e ) 


(3.72) 


(3.73) 


(3.74) 


(3.75) 


and N' = N Y^jLi fj the total number density of electron oscillators. This is identical to the 

refractive index of a plasma at high frequency (Eq. 3.71). 


3.4.4 Dispersion of pulses from a pulsar 

Pulsars (magnetised neutron stars) emit one or more narrow radio pulses every rotation period, 
rather like a lighthouse. The interstellar medium is partly ionised and the electron number 
density varies throughout the Galaxy - the average value in the Galactic plane is n e ~ 3 x 10 -7 
m -3 giving u p = uj p /2tt ~ 50 kHz. Pulses travel at v g (v), so high frequencies arrive earlier than 
low frequencies, with pulses at extremely high frequencies travelling at c. 

Pulsars are often observed by radio telescopes at two or more frequencies, typically at 650 MHz 
to 700 MHz frequencies, well above the plasma frequency. The time delay of a pulse at angular 
frequency uo from a pulsar at distance L with respect to the same pulse observed at angular 
frequency uo oo is 


m = f 1 

Jo 


dr 


v g (r,uj) 


L 

c 


(3.76) 
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where, for uo cj p as is the case for radio observations, we obtain from Eq. 3.69 and Eq. 3.66 


v g (r,uj) 


1 L _ KK)] : 


2 \ - 1/2 




1/1 [w p (n e )] 2 
2 w 2 


Hence, 


1 L n e ( r)e 2 1 
c l 2som e uj 2 


(3.77) 


t(w) = 


2£ora e c 


DMw“ 2 = 


8TT 2 £QfTl e C 


DM v 


-2 


(3.78) 


where the dispersion measure is 


DM = / n e (r)dr. 
Jo 


(3.79) 


The dispersion measure can be found from the pulse delays at two or more frequencies, giving 
the pulsar’s distance L if n e ( r) along the line of sight to the pulsar is known. If L is known, 
e.g. by parallax or other methods, the observed dispersion measures of pulsars throughout the 
Galaxy can be used to refine models of n e (r). 


Identical physics applies to radio pulses propagating through the ionosphere between, e.g. Global 
Positioning System (GPS) satellites, and the Earth. In this case the relevant quantities are the 
vertical total electron content 

POO 

VTEC = / n e (h)dh , (3.80) 

J o 

where h is the altitude, and the slant total electron content STEC=VTECx sec 6 where 9 is 
the zenith angle of the satellite as seen by the tracking station. VTEC and STEC are usually 
measured in units of 1 TECU= 10 16 electrons m -2 . 


3.4.5 Faraday rotation and applications to radio astronomy 

When a linearly polarised wave propagates through a magnetised plasma (or dielectric), its 
plane of polarisation rotates. This effect was observed by Michael Faraday for polarised light 
passing through a dielectric in the presence of a magnetic field. It is caused by the phase velocity 
being different for right and left circularly polarised waves. We can most easily show this for 
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the case of a dilute plasma in which there is a component of magnetic field present along the 
direction of wave propagation. 

For convenience let’s assume that B = B yz and that the wave propagates in the +z direction. 
Then using the (unnormalised) basis vectors for circular polarisation (Eq. 2.67), at z = 0 the 
electric field of a right-circularly polarised wave is 

E(x, y, 0, t) = E r e R e~ lwt = Er (x — iy)e~ lut . (3.81) 

Furthermore, let’s assume that electrons located at z = 0 will undergo circular motion in the 
constant magnetic field in response to the electric field of the wave. For convenience consider 
an electron orbiting around the origin. Its position will be 

r (t) = r R e R e~ lujt = r R (x - iy)e~ lut , (3.82) 

Re{r(£)} = tr cos (o;t)x — tr sin (u;t)y (3.83) 

which clearly describes circular motion with radius tr. 
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Applying the Lorentz force law m e r = — e (E + r x B) to an electron orbiting in the xy plane 
around the origin in the presence of magnetic field B = B nz and a right-circularly polarised 


(3.84) 

(3.85) 

(3.86) 

(3.87) 

(3.88) 


wave propagating in 

the 5 

5 direction (Eq. 

3.81). 

> 


—m e 

Lo 2 r R e R e 

-iujt 

= —e (Er 

e~ iojt 

— iouvReRe lut x Byz) 

m e 

uj 2 r R (x - 

*y) 

= e [Er(x- 

iy) ~ 

!wr fl (x 

- iy) x B|,z] , 

m e 

uj 2 r R (x - 

*y) 

1 

<X, 

II 

iy) - 

iur R (- 

y-ix)P||] , 

m e 

uj 2 r R (x - 

*y) 

1 

<X, 

II 

iy) - 

ur R (St- 

- *y)E||] • 




eE R 






r R 

m e (jj 2 + ecu 

iBu' 




Using this result we can immediately write down the dipole moment pr — —erR, polarisation 
field Pr — XR^oEr , electric susceptibility xr > permittivity sr = (1 + Xr) 8 o> refractive index 
tir ~ y] sr/ <Sq and wave number Pr = lj tir/c for number density n e of electrons, 


PR = -e 


cEr 


m e (jj 2 + ecj Bn ’ 


_ n e e 9 
Pr = ——W 2 


m e 

,2 


( 2 + v) 


Xr 


£r 


«R = 


V?) 

- (i + —) 
uP V w / 


cu 2 


£0> 
1 1/2 


7 W 

kr = — 
c 


l ~ ~2 ( 1 + —) 1 
a ;" 1 V cj / 


1 1/2 


(3.89) 

(3.90) 

(3.91) 

(3.92) 

(3.93) 

(3.94) 


where uo c = eB\\/m e is the electron cyclotron frequency and ui p is the plasma frequency (Eq. 3.66). 


Provided w and ui ui c 


k R » — 

c 


i- l 4( 1 + - 

2 uj z \ co j 



CJ 


r ^ _ 

/ 

c 


1 

2 uT 2 


1 - - fl - — 

UJ 


(ko + Ak), (3.95) 
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uo 1 oof, \ 

where ko « — 1 —- 77 is the high-frequency approximation for the wave number in an 

c \ 2 u z J 

unmagnetised plasma and 


A k = 

2 uj z c 


(3.96) 


Similarly, for left-circular polarisation 


k L 


uo 

c 


l- 1 ^ (i+ — 


2 w 2 


(^) 


{ko- Ak). 


(3.97) 


Now consider a plane-polarised wave propagating in the z direction and being initially polarised 
in the x direction, i.e. E(x,y, 0,t) = Eo9.e~ lujt — (£'o/2)(ei / + After propagating 

distance 2 : the electric field is 


E(x, y, z,t) = f ( e L e lkhZ + e R e 


ik R z\ e ~uot 


Ec\ r 


= ^T (X + iy)e~ iAkz + (Z-iy)e 


AkAk z 


2 L 


O i(koz—ujt ) 


^iAk z _ ^—iAkz 


= T [ S ( eiAkZ + e~ iAkz ) - iy ( 

= E 0 [x cos(A k z)+ y sin(A k z)\ e i{ - koZ ~^. 


)] 


D i(koz—ujt) 


(3.98) 

(3.99) 

(3.100) 

(3.101) 


Hence the plane of polarisation at frequency v has rotated with respect to that at v —> oo by 
angle 


9 = A k z 


1 , ,2 

1 (dp L>J C 

2 c j 2 c 


1 n e e 3 1 

2 SqUI^C (jO 2 


2_ A 2 R„ a 

87T 2 SqUI^C 3 


(3.102) 


Since the high-frequency behaviour of a dielectric is essentially identical to that of a plasma, 
the formula above applies equally to dilute plasmas or dielectrics at high frequency. 


Faraday rotation is a valuable tool which radio astronomers use for measuring the magnetic 
field in the interstellar medium or in clusters of galaxies. By observing the polarisation angle 
of radio emission from a polarised background source such as pulsar or quasar at two or more 
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frequencies, one can solve for the the rotation measure 
0 1 77 p^ C 

RM = — = / n e (r)B(r)-dr (3.103) 

8n z £om z c 6 J 

= 2.63 x 10 -13 J n e (r)B(r) • dr (rad m -2 ), (3.104) 

along the line of sight. In the case of pulsars, dividing J n e {r) B(r) • dr obtained from the RM 
by the dispersion measure DM = J n e (r ) dr (Eq. 3.79) gives the weighted average value of the 
component of magnetic field along the line of sight (weighted by the free electron density). 


Summary of important concepts and equations 


Complex wavenumber and absorption 

— Electric field, E(r,t) = E 0 e i ( Re{fc}K ' r “ wt )e- Im ^> £ r 
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— Intensity, / = I^e °^' r 

— Absorption coefficient, a = 2 x I m{k(uo)} 


Lorentz Oscillator Model of Dielectric 


— Nucleus surrounded by electron “cloud” acting as simple harmonic oscillator: resonant 
frequency cjo; damping constant 7 

~ d 2 


— Equation of motion of the electron m e 


W + + 7 ^ 


r e (t) = —eE(t) 


— For fj electrons with resonant frequency coj and damping coefficient 7 j , and N molecules 
per unit volume, 


PM) 

n(uS) 


Ne 2 


m- 


E 


1 + 


, J'=1 3 

Ne 

J 


00 1 - — UJ 2 — i^jUJ 


E(r,i). 


E 


fj 


1/2 


£o m e ^ co 2 - - co 2 - i'yjbJ 


— Absorption coefficient near a resonance has a Lorentzian profile 


M) 


1 (7j/ 2 ) 

t. ( Uj - co) 2 + (7i/2) 2 


Dispersion 


— Phase propagates at the phase velocity v p (uj) 


uj 


k 


UJ 


— Wave packets travel at the group velocity v g (u) 

— Information and energy travel at v g 


du 

dk 


UJ 


— If (dKe{n}/dcj) < 0 (anomalous dispersion) the group velocity can exceed c, but anoma¬ 
lous dispersion is accompanied by absorption, so causality is not violated 


Refractive index of a conductor 
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<9E 


d 2 E 


— Wave equation: Y E — lmj— -= 0 

at ot z 

— Dispersion relation for good conductor: k{uS) (l + i)2- 1 /2( /xaw )l/2 

2 \ 1/2 


Skin depth S(cu) = l - ) 

— Absorption coefficient a(cu) = 2 /S(cj) 


Wave propagation in dilute plasmas 

— n e electrons and rii ions per unit volume - ignore ions in dilute plasma 

o 2 \ i q2 E 


Wave equation V 2 E — - I E 

m P 


— Dispersion relation k 2 c 2 — uo 2 — uo 2 


z 2 dt 2 


= 0 


Plasma frequency uo p = 


n e e^ 


£o m e 


2 \ 1/2 


EM waves cannot propagate if uo < uo p (cut-off frequency) 

Above u p : v p > c and v g < c 

At frequencies above the highest resonance frequency a dielectric behaves like a dilute 
plasma 


Faraday rotation 


— In a magnetised plasma or dielectric, left and right circularly polarised waves travel at 
different phase velocities 

— For plane-polarised waves this causes the plane of polarisation to rotate by angle 

2 m~c uj z 11 


where Bn is the component of the magnetic field parallel to the wave direction 
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Exercises on Chapter 3 

3-1 The resistivity of silver is p = 1.6 x 10 s Q m, and its permeability is //, = 0.9998/to, find 
the reflectance for light of wavelength 500 nm at normal incidence from air. 

3-2 Show that the time-averaged power density (E • J) in a dilute plasma is zero. 

3-3 At midnight and noon, at a certain location and date, the electron number density in the 
ionosphere was as given in the plot below. 

(a) Label the x-axis at the top of the plot in terms of the plasma frequency v v (MHz). 

(b) Find the height of the reflecting layer at 1 MHz and 3 MHz at midnight and noon. 

(c) Find the minimum frequencies for communication with orbiting satellites at midnight 
and noon. 

(d) Estimate the total electron content in (electrons m -2 ) at midnight and noon. 

(e) An instantaneous broad-band pulse of radio-frequency interference (RFI) is emitted 
overhead by an orbiting satellite and observed by an terrestrial detector with a bandwidth 
from 1.2—1.8 GHz. What is the duration of the observed pulse? 
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log[A7 e /m' 3 ] 

3-4 For pulsar PSR J1817-3837 the observed position angle (PA) of linear polarisation with 
respect to the North Celestial Pole is plotted below vs. frequency. The dispersion measure 
is 102.85 pc cm -3 , where 1 pc = 3.09 x 10 16 m (1 parsec) is the distance unit used by 
astronomers. Find the average value of the parallel component of the magnetic field along 
the line of sight to this pulsar. 


PSR J1817-3837 



3-5 We are receiving radio signals from an interplanetary space probe which is far from the 
Sun. It is currently on a trajectory such that it has been eclipsed by the Sun, but it is now 
emerging from behind the Sun. The space probe regularly broadcasts instantaneous broad¬ 
band pulses of radio-frequency emission every second, and it also broadcasts spacecraft 
experiment instrument data over a narrow frequency band of width 100 kHz extending from 
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10,000.0 MHz to 10,000.1 MHz. It transmits using a dipole antenna aligned perpendicular 
to the plane of the ecliptic (the plane containing the Earth’s orbit). We are just renewing 
radio contact as the space probe is emerging from behind the solar limb. 

The Sun’s corona is an extremely hot (> 10 6 K) plasma which is highly variable, and 
has dynamic coronal loops of magnetised plasma. The corona is located above the Sun’s 
photosphere. At the time we are observing, the base of the corona in front of the space probe 
as it emerges from behind the Sun has an electron number density of n e (i?©) ^ 3x 10 15 m -3 . 
The corona’s a scale height is about H ~ 10 8 m, such that the electron number density 
decreases with height as n e {r ) = n e (i?©) exp[—(r — Rq)/H]. Imagine that, at the time we 
are observing, the Sun’s magnetic dipole moment is m — 2 x 10 29 A m 2 and it happens to 
be pointing directly towards Earth. 

(a) What is the minimum radio frequency that we are able to receive from the space probe? 

(b) What is the additional delay of the pulses due to propagation through the corona? 

(c) By how much has the plane of polarisation of the EM wave carrying the data signal 
rotated due to propagation through the corona? 

[Make what you think are reasonable approximations - do not attempt a rigorous calcula¬ 
tion. The radius of the solar photosphere is i? 0 = 6.955 x 10 8 m.] 
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4 Waveguides 






Figure 4.1: An assortment of waveguide components. (Image reproduced with permission, 
courtesy of Altair Technologies, Inc., http://www.altairusa.com/index.php). 


Learning objectives 

— To be able to show that the boundary conditions at the surface of a perfect conductor 
forbid transverse electromagnetic (TEM) waves propagating inside a waveguide, but allow 
transverse electric (TE) waves and transverse magnetic (TM) waves. 

— To understand that the field of a monochromatic EM wave propagating inside a rectan¬ 
gular waveguide with k = kz can be written 

E (x,y,z,t) = E 0 (x,y)e l< ' kz ~ ut \ B (x,y,z,t) = B 0 (x,y)e^ kz ~ ujt \ 

where the field amplitude functions Eo (x,y) and Bo (x,y) can be found by substitution 
into Maxwell’s equations in differential form and applying the boundary conditions. 

— To know that by using Faraday’s and Ampere’s laws we get equations for the x and y 
components of Eq (x,y) and Bq (x,y) in terms of the z components, and that by subse¬ 
quently using V • E = 0 and V • B = 0 we obtain partial differential equations for the 
z components of Eq (x,y) and Bq (x,y) which can be solved by separation of variables 
leading to various discrete TE and TM modes. 
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— To be able to apply the boundary conditions to obtain the dispersion relation, phase and 
group velocities for each of the TM mn and TE mri modes, and show that each mode has a 
cut-off frequency the lowest being for the TEio mode. 


Waveguides are devices used to guide electromagnetic waves. At optical frequencies these would 
be optical fibres. At microwave frequencies these are empty hollow metal pipes and it is these 
we shall consider in this section. Examples of waveguide components are shown in Fig. 4.1. 


We shall approximate the metal by a perfect conductor, and inside a perfect conductor E = 0 
and B = 0, so that the boundary conditions on the electromagnetic fields at the metal surface 
just outside the conductor are simplified to 

E 11 =0, B- 1 = 0, (4.1) 


i.e. the electric field at the surface is normal to the surface, and the magnetic field at the surface 
is parallel to the surface. 


We shall show that the boundary conditions do not allow transverse electromagnetic ( TEM ) 
waves such as those present in empty space to propagate inside a long (length L A) narrow 
(width W < A) hollow metal pipe. Either E or B must have a longitudinal component for an 
EM wave to propagate. We can easily prove this by contradiction, i.e. by trying to find a TEM 
solution for, say, the case of a hollow metal pipe of circular cross section. Remembering that 
magnetic field lines have no beginning or end, and that the boundary conditions require that 
the magnetic field at the surface of a (perfect) conductor to be parallel to the surface, I sketch 
the simplest magnetic field line allowed both by the boundary conditions and by the cylindrical 
symmetry, and having B transverse to the axis of the pipe (Fig. 4.2a). 



Figure 4.2: A waveguide of circular cross section showing attempts to construct the fields of 
a TEM wave: (a) a magnetic field allowed by the boundary condition and symmetry, (b) an 
electric field allowed by the boundary condition and symmetry but ruled out unless there is a 
central condutor as in part (c). 
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Remembering that electric field lines start on positive charge, the pipe is empty and the only 
possible location of charge is on the surface of the conductor, and that the boundary conditions 
requires the electric field at the surface of the (perfect) conductor to be normal to the surface, I 
sketch the simplest electric field lines allowed by the boundary conditions and by the cylindrical 
symmetry while having E transverse to both B and the axis of the pipe (Fig. 4.2b). There is 
clearly a problem because electric field lines must end on negative charge but the pipe contains 
no charge other than on the surface. 

The only way to solve the problem while maintaining both E and B transverse to the pipe’s 
axis is to introduce a central conducting wire (Fig. 4.2c), but then we no longer have an empty 
pipe - instead we have a coaxial cable. A coaxial cable does act to guide TEM waves, and is 
therefore strictly speaking a waveguide, but is not what we usually mean by a “waveguide”, 
i.e. a hollow conducting pipe. Of course one can shine a narrow laser beam along the axis of a 
metal pipe of large diameter, but in this case the pipe is not actually guiding the laser’s TEM 
wave. What we are really talking about is the impossibility of a TEM wave with coherence 
width comparable to or larger than the width of the pipe, so that the pipe could be considered 
as genuinely guiding a monochromatic TEM plane wave along the pipe. 

Although our proof was for a pipe of circular cross section, we shall shortly prove that TEM 
waves can’t propagate down a rectangular cross-section pipe, and it can also be proved for long 
hollow metal pipes of any shape. There must be a component of either the magnetic field or 
of the electric field along the pipe. If B has a component along the waveguide then only E is 
transverse, i.e. we have a transverse electric (TE) wave. Similarly, if E has a component along 
the waveguide then only B is transverse, i.e. a transverse magnetic (TM) wave. 


4.1 Waves in rectangular waveguides 


Rectangular cross-section metal pipes are used to guide EM waves at microwave frequencies, 
say to or from a transmitting or receiving antenna. We shall consider monochromatic plane 
waves propagating down a such rectangular cross-section waveguides of infinite length in the 
^-direction (Fig. 4.3). The electromagnetic field of a monochromatic EM wave inside the waveg¬ 
uide is of the form 

E (x,y,z,t) = E 0 {x,y)e i{ ~ kz ~^\ B (x,y,z,t) = B 0 (x, y)e i{ - kz ~^\ (4.2) 
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z x 

Figure 4.3: Waveguide geometry. 


where 


E 0 (x,y) = E°(x,y)5i + E°(x,y)y + E Q z {x,y)z, 
B o(x,y) = B°(x,y)x + B°(x,y)y + B° z {x,y)z. 


(4.3) 

(4.4) 
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We need to find what we might call the “field amplitude functions”, E®(x, y), Ey(x, y), E®(x, y), 
B®(x,y), By(x,y) and S^(x,y), that specify the EM field at (x,y,z = 0) at time t — 0 that 
satisfy Maxwell’s equations inside the waveguide as well as the boundary conditions on the EM 
field at the conducting surface. 


Putting Eqs. 4.2 into Faraday’s law, 


^ _ dB 

V X E= ~~9t’ 


(4.5) 


for e lujt time-dependence and e lkz ^-dependence Faraday’s law gives 


a ikz 


dE z 

dy 

ii 

1 

dB x 

dt 


o 

1 

ikEy 

= iuB° x , 

(4.6) 

dE x 

dz 

dE z _ 
dx 

dBy 

dt 

=> 

ikE° x - 

9E° Z 

dx 

= iuB°, 

(4.7) 

dEy 

dx 

dE x _ 
dy 

dB z 

dt 


9E° y 

dx 

9E° X 

dy 

= wB° z . 

(4.8) 


Putting Eqs. 4.2 into Ampere’s law for empty space, 


„ „ i dE 

VxB = ?aP 


(4.9) 


for e lujt time-dependence and e lkz z-dependence Ampere’s law gives 


a ikz 


dB z 

dBy 

1 dE x 


9B° Z 

— ikBy = 

~ iuj E o 

(4.10) 

dy 

dz 

c 2 dt 

dy 

c 2 x: 

dB x 

dz 

dB z 

dx 

1 dE y 
c 2 dt 


ikB° x 

II 

~ iUJ E° 

c 2 ^V’ 

(4.11) 

dBy 

dB x 

1 dE z 

=> 


9B° X 

~ i0J E° 

(4.12) 

dx 

dy 

c 2 dt 

dx 

dy 

c 2 hz ' 
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Equations 4.6 and 4.11 give 


BS. = I ^ ) (ifc)- 1 , 


f! = I ikB! - f “ 


dx 


-1 


from which we can eliminate E® to obtain B 


0 = » /.ag _ “dEl\ 

x ( uj/c ) 2 — k 2 \ cte c 2 9?/ / 


Similarly, one can show that 


„o = i " Ml 

27 ( u/c ) 2 — k 2 \ c 2 dx ) ' 


(u;/c) 2 — fc 2 y dx ^ dy ) 


, <9M> 


£% = 


{oj/c ) 2 — k 2 \ dy 


, <9M> dB° z 
k—-?- -oj 2 


cte 


Thus, if we know the longitudinal fields we can obtain the transverse fields. 


(4.13) 

(4.14) 


(4.15) 


(4.16) 

(4.17) 

(4.18) 


We have already shown that TEM waves cannot occur in a hollow circular cross-section waveg¬ 
uide. We are now in a position to prove that TEM waves are also impossible in a rectan¬ 
gular cross-section waveguide. Putting both E® — 0 and B® = 0 into Eqs. 4.15-4.18 gives 
B° x = B° y = 0 and E° x = E° y = 0, proving that unless there is a longitudinal component of E or 
B that there will be no EM wave inside the waveguide. 


4.2 Waveguide modes 


Putting Eqs. 4.2 into Gauss’s law in empty space, V • E = 0, we obtain 


SEl 9E$ 

dx dy 


+ ikE° z = 0 


(4.19) 
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Next, substituting Eqs. 4.17 and 4.18 we obtain 

d 2 E2 d 2 B ?' 


(w/c) 2 — k 2 dx 2 dxdy 

i / d 2 £° _ q2b z 
( uj/c) 2 — k 2 \ dy 2 U dydx 

d 2 d 2 . . , 2 

+ v + ( “ /c) 


-\- ikE g — 0 , 

(4.20) 

tq 

iN! O 

II 

o 

(4.21) 


Similarly, putting Eqs. 4.2 into V • B = 0 and substituting Eqs. 4.15 and 4.16 yields 


^2 + ^2 + W C ) 2 - ^ ) B * ~ °- 


(4.22) 


Only two types of EM wave are possible in a waveguide: TM (“transverse magnetic”) with 
E z = 0, and TE (“transverse electric”) with B z = 0 — to find their allowed EM fields we will 
need to solve Eqs. 4.21 and 4.22, respectively. 
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4.2.1 TM modes 

TM waves have B z — 0 and E z ^ 0. We can solve Eq. 4.21 for E®{x,y) by the method of 
separation of variables. We seek a solution of the form 


E° z (x,y) = X(x)Y(y) 

(4.23) 

and substitute it into Eq. 4.22 to get 


r d J +x Z + ^~ k2 ' x r = o, 

(4.24) 


(4.25) 

Since the 1st term depends only on x, and the 2nd term depends only on y, 

constant such that 

each term must be 

1 d 2 X 2 

X dx 2 ~ Kx ’ 

(4.26) 

1 d 2 Y 2 

Y dy 2 ~ Ky ' 

(4.27) 


and 


~ K l + [( w / c ) 2 - k 2 } = 0. (4.28) 

Equations 4.26 and 4.27 have sinusoidal solutions for X(x) and Y(y) of the form 

X(x) — F sin(K x x) + Gcos(k x x), (4.29) 

Y(y) = E sin {n y y) + Q cos (K y y) . (4.30) 


We now need to find allowed values of F, G, J 7 , Q , k x and K y from the boundary conditions on 
the electric field at the surface of the conductor Ell = 0 (Eqs. 4.1), i.e. at x = 0, x — a, y = 0 
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a 


Figure 4.4: Cross-section of rectangular waveguide with width a and height b. 

and y — b in Fig. 4.4. This condition is satisfied by G = Q — 0 and 

E®(x,y) E mn sin ^ m7VX ^j s [ n for m = 1, 2, 3,... and n = 1,2, 3,_ (4.31) 

This solution gives the amplitude function of E z for the TM mn mode; the 1st index m by 
convention is associated with the larger dimension and the 2nd n with the smaller dimension. 
We have now obtained the solution of Eq. 4.21 for TM waves. Note that the lowest possible 
TM mode is TMn. 


4.2.2 TE modes 


TE waves have E z — 0 and B z / 0. We can solve Eq. 4.22 for B®(x,y) by the method of 
separation of variables 

B° z (x,y)=X(x)Y(y) (4.32) 

as for the case of the TM modes, and with the same form for the solutions X(x ) and Y(y) (as 
given by Eqs. 4.29 and 4.30). We now need to find allowed values of F, G, F, <5, k x and K y for 
TE modes from the boundary conditions on the magnetic field at the surface of the conductor 
= 0 (Eqs. 4.1), i.e. at x — 0, x = a, y = 0 and y — b in Fig. 4.4. 


Since the waveguide is aligned along the z direction the ^-component of B is always parallel to 
the conducting surface and so B^ = 0 is always satisfied by B z . To get a useful constraint on 
X(x ) we must instead use B x with F>^(0,?/) = B®(a,y) = 0 as required by B^ = 0. Eq. 4.15 
gives B® in terms of B® and but since for TE waves E z — 0 we obtain 


B x( x ,y) 


f ik \ dB^jx^y) 
\(cc;/c) 2 — k 2 ) dx 


(4.33) 
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Hence, 


dB z (x, y ) _ q (^at x = 0 and x = a), (4.34) 

ox 

and differentiating Eq. 4.29 with respect to x we obtain 

k x F cos(k x x) — K x Gsin.(K x x) = 0 (at x = 0 and x — a) (4.35) 

which requires F = 0 and k x — run/a for m = 0,1, 2 ,.... 

Similarly, To get a useful constraint on Y(y) we must instead use B y with 5^(0, y) — S^(6, y) = 
0, and Eq. 4.16 which gives By in terms of B z and E z (0 for TE waves). The result is that 
T — 0 and k v — nn/b for n = 0,1, 2,.... We have now obtained the solution of Eq. 4.22 for 
TE waves 

B z (x,y) = B mn cos ^ m7rx 'j cos ^ for m = 0,1, 2,... and n — 0,1,2,_ (4.36) 
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Note that either m > 0 or n > 0 otherwise B®(x,y) is constant and Eqs. 4.15-4.18 yield 
E x = Ey = B x = By = 0. Eq. 4.36 gives the amplitude function of B z for the TE mri mode. We 
have now obtained the solution of Eq. 4.22 for TE waves. 


4.3 Dispersion relation, phase and group velocities 


We can find the dispersion relation from Eq. 4.28 

k 2 c 2 = u 2 — k 2 x c 2 — k 2 c 2 (4.37) 

kc = (4.38) 

where 

^ --V , 0 2+ © 2 ' <4 - 39> 

If uo < (jJmn the wave number is imaginary, and the TE mn and TM mn waves do not propagate 
(they die off exponentially), so cc mri is the cut-off frequency for the TE mn and TM mn modes. 
Since there is no TMio, the lowest cut-off frequency occurs for the TEio mode and is cjio = cn/a. 
The complete EM fields for the TE and TM modes are given in Appendix 4.36. 


The phase and group velocities for TE and TM modes are 


uj 


k i/l - (w mn /w ) 2 


> c, 


b- = c\f\ — (w mn /w ) 2 < c, 


dk/dui 


(4.40) 

(4.41) 


and they are plotted in Fig. 4.5. Note that v g /c = c/v p (similar to a dilute plasma). 


When using a waveguide to transmit a signal at a given frequency, it is important that only one 
waveguide mode is excited otherwise more than one copy of the signal will propagate along the 
waveguide, each with a different velocity. From Fig. 4.5 we see that this can be accomplished if 
the frequency is between zqo and z/20- Waveguides with different widths are therefore designed 
for use at different frequencies, e.g. a standard waveguide for the “X-band” (8.2-12.4 GHz) has 
internal cross section 2.286 cm by 1.143 cm. 
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Figure 4.5: The first few modes of a rectangular waveguide with a = 2b. 

4.4 The EM field and power transmitted for TEi 0 mode 


For all the TE modes the amplitude function of the longitudinal component of the magnetic 
field is given by Eq. 4.36. We can get all the other field amplitude functions from B®(x,y) 
and E®(x, y) = 0 using Eqs. 4.15—4.18, and then multiplying by exp [i (kz(cu) — cut)} to get the 
(complex) EM field, and finally taking the real part. In this way one can find the non-zero EM 
field components of a TEio wave travelling in the ^-direction: 


En, - — 


( CL\ / 7VX \ / 7TX \ 

—J sin y—J sin (kz — cut) = E max sin y—J sin (kz — cut), 


B x 

B 7 


£>io £ (—1 sin f) sin (kz — cut) — — E max — sin f) sin (kz — cut), 
\7T/ V a J cu V a J 


( 7VX \ 71/ 7VX \ 

— ) cos (kz — cut) = — E max — cos ( — ) cos (kz — cut), 
a J acu \ a J 


where 


k = yw 2 /c 2 - (7T/a) 2 , E ma , x = -B w —. 

7r 


(4.42) 

(4.43) 

(4.44) 


(4.45) 
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The electric and magnetic fields for the other TE modes, and for the TM modes, can be obtained 
in a similar way from E®(x,y) and B® — 0. The electric and magnetic fields for all TE and 
TM modes are given in Appendix D, and for example, Fig. 4.6 shows the field configuration for 
TEio and TE 20 waves travelling in the z direction. It is straightforward to show (Exercise 4-2) 
that the energy density per unit length of waveguide for the TE 10 mode is 


U 

L 


Sq ab 

~T~ 


Ei 


(4.46) 



Figure 4.6: The electric field and magnetic field configurations of TE10 and TE 2 o waves 
travelling in the z direction. The electric field is in the =by direction and the colour image shows 
it’s amplitude (black corresponds to E = —E max y and white to E = +E max y). The magnetic 
field direction and strength is shown by the vectors. 
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The Poynting vector is 


S(r ,t) = E(r ,t) x H(r ,t), 

= [E y (r,t)y] x — [B x (r,t)Sc + B z (r,t)z ], 

Mo 

= — [~E y (r,t)B x (r,t)z + E y (r, t)B z (r, f)x], 

Mo 


1 

Mo 


v, 


f-^maxsin 2 ^ sin 2 (fc2! — cut)' 1 . 


7 r 

2acc; 


sin 


27TX 


sin[ 2 (kz — ut)\ x 


S(r, t) = ^£o^max s i n2 ) sin 2 (/c2: — cot)\ 


1 7F n o . (21TX\ 

»2^ £ ”“ S,11 (~l S,ll[2(fa ' a,f)lx - 


(4.47) 

(4.48) 

(4.49) 


(4.50) 


(4.51) 
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z 


(4.52) 


with v g = c 2 k/u: and k = y/cc 2 /c 2 — ( n/a ) 2 . Note also that we have used sin 9 cos 6 — \ sin(20) 
above. 


The power transmitted is the integral of the flux of (S) over the cross-sectional area of the 
waveguide 



•6 


—) dydx , 


P = 


a / 


(4.53) 



dx. 


(4.54) 


P = abv g ^El 


max* 


(4.55) 


Notice that this is just the energy per unit length (Eq. 4.46) multiplied by the group velocity, 
and we see that the energy propagates along waveguide at v g , as expected. 

Summary of important concepts and equations 

Boundary conditions at the surface of a perfect conductor (Ey = 0, = 0) 

— Forbid transverse electromagnetic (TEM) waves which have B • k = 0 and E • k = 0 

— Allow transverse electric (TE) waves which have E • k = 0 and B • k ^ 0 

— Allow transverse magnetic (TM) waves which have B • k = 0 and E • k ^ 0 

Field of a monochromatic EM wave inside rectangular waveguide (k = kz ) 



— Field amplitude functions Eq(x, y ) and Bq(x, y ) can be found by substitution into Maxwell’s 
equations in differential form and applying the boundary conditions. 
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— Using Faraday’s and Ampere’s laws we get equations for the x and y components of 
Eq (x,y) and Bq (x,y) in terms of the 2 : components. 

— Using V-E = 0 and V B = 0 we obtain partial differential equations for the z components 
of Eo (x,y) and Bo (x,y) which are solved by separation of variables. 

— Applying boundary conditions give the solutions for TM mn and TE mn modes 

TM: E° z (x,y) = E mn sm ) sin , (m = 1,2,3,...; n = 1,2,3,...), 

TE: B° z (x,y) = B mn cos ) cos , (m = 0,1,2,...; n = 0,1,2,...). 

Dispersion relation, phase and group velocities 



— uj mn is the cut-off frequency for the TE mn and TM mn modes. 
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— The lowest cut-off frequency occurs for the TEio mode and is cj 10 = cn/a. 

— The phase and group velocities for TE and TM modes are 

uo c 

P k \/l - kn/w) 2 

^ = dCfcLo = c Vi - (“W ") 2 < c. 

Exercises on Chapter 4 

4-1 A standard waveguide for the “X-band” (8.2-12.4 GHz) has internal cross section 2.286 cm 
by 1.143 cm. Find the the first two cut-off frequencies, as these will give the frequency 
range for which only the TEio mode will propagate. 

4-2 Find the electromagnetic field energy per unit length of waveguide for TEio mode. Express 
your result in terms of afeo^max- Start with the fields for the TEio mode 

( 7TX \ 

—J sin {kz - ut), 

/c / 7 tx \ 

B x = — E/ max — sin I— )sm{kz-ut), 
uj V a / 

7T / 7 YX \ 

B z —E max —cos ( — ) cost kz — cut) 

au V a J 

where k — \Zcd 2 /c 2 — ( 7 r/a) 2 . 

4-3 Consider a hollow standard WR-159 F-band waveguide which has an internal cross section 
of 40.386 mm x 20.193 mm. 

(a) A 2 GHz signal is fed into the waveguide. Will this cause a wave to propagate along 
the waveguide? 

(b) Assuming the wave would be in the TEio mode, what is the wavenumber? Discuss the 
implications of your answer. 

4-4 Consider the case of launching an EM wave with frequency uo < cjio into a rectangular 
waveguide. Assuming the wave will attempt to propagate as a TEio wave travelling in 
along the waveguide (in the £ direction), 

(a) Derive all the components of the electromagnetic field for the case where k is purely 
imaginary. 

(b) Obtain the Poynting vector and discuss the energy flow. 


(4.57) 

(4.58) 


(4.56) 


130 


Download free eBooks at bookboon.com 







Essential Electrodynamics 


Radiation and scattering 


5 Radiation and scattering 

Learning objectives 

— To understand that an accelerated charge radiates, and to be able to obtain Larmor’s 

a 2 q 2 

formula which gives instantaneous power P =--. 

6n£o c 6 

— To be able to derive the time-averaged power radiated by an oscillating dipole p (t) = 

dP 
dQ 

sin 2 6 dependence. 

— To be able to show that a short centre-fed dipole antenna (total length d<C A) behaves 
as an oscillating electric dipole, but know that half-wave antennas with d — A/2 are more 
efficient radiators and have radiation patterns which resemble the standard dipole pattern. 

— To understand that EM waves incident on a particle (or small object) will produce os¬ 
cillating charges or currents causing some of the incident radiation to be scattered (i.e. 
absorbed and re-radiated). 

— To learn that the total cross section for scattering of EM waves is the time-averaged 
power re-radiated divided by the time-averaged incident energy flux, and that for free 
electrons acceleration and re-radiation according to Larmor’s formula leads to Thomson 
cross section for scattering of EM radiation or photons with frequencies v <C m e c 2 /h by 
free electrons 

— To know that for a molecule, an EM wave induces an oscillating molecular dipole <a m E(t) 
and a cross section which is proportional to A -4 thereby explaining why the sky is plue 
and sunsets are red. 

— To be able to derive the inhomogeneous wave equations for the potentials, their retarded 
potential solutions and to derive the vector potential for sinusoidal sources. 

— To be able to apply the retarded potential solutions in the far zone to obtain the vector 
potential and EM field for a long antenna. 


po cos(c<;£)z, and its radiation pattern ( 


Pq lj 4 sin 2 6 
327r 2 £o c 3 


showing the characteristic 


Whenever you accelerate a charge it radiates EM radiation. Electrons oscillating up and down 
an antenna are being accelerated and radiate radio waves. In an atom in transition from one 
quantum state to another there is a superposition of the two states with the electron’s probability 
density oscillating and changing shape at a frequency equal to the difference in energy levels 
divided by Planck’s constant, the result being the emission or absorption of a photon. 
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5.1 Larmor's formula 


Consider a stationary isolated point charge +q at the origin. Its electric field lines extend 
radially outward isotropically to infinity as shown in Fig. 5.1(a). Suppose at time t — 0 we 
accelerated that charge with acceleration a. — az for an infinitesimal time interval At such that 
it acquires velocity v = vz much less than the speed of light w«c. 




Figure 5.1: (a) Electric field lines of stationary point charge +q located at the origin, (b) 
Electric field lines of the charge at time t after it has been accelerated during an infinitesimal 
time interval from t = 0 to t = At, after which it reached velocity v = vz where v«c. 


Because information cannot be transmitted faster than light, if we look at the charge at time 
£, when it has moved to position r = vtz (Fig. 5.1b), we would see that at a distances r > ct 
from the origin the electric field lines would be identical to those of the charge before it had 
been accelerated. However, because v <C c, within distance r < c{t — At) from the origin, the 
electric field lines at time t would closely resemble those of a stationary charge +q located at 
r = vtz. The field lines inside radius c(t — At) must join up with those outside radius ct. This 
takes place in the spherical shell between c(t — At) < r < ct where there will be a kink in the 
field lines which propagates outwards from the origin as a spherical wave travelling at speed c. 


By taking a close-up look at one of the field lines in the vicinity of the kink (Fig. 5.2), i.e. in 
the spherical shell between c(t — At) < r < ct, we find that 


Eq vt smO (aAt)(r/c)sm9 r 

— = --— = --- = a^r smO. 

E r cAt cAt c 2 


(5.1) 


132 


Download free eBooks at bookboon.com 

















Essential Electrodynamics 


Radiation and scattering 



Figure 5.2: Enlargement of one of the electric field lines of the charge at time t after it has 
been accelerated for a short time interval from t — 0 to (t-\- At) after which it reached velocity 
v = vz where ^<c. 


where a = v/At is the acceleration. The radial component of E is given by Coulomb’s law, 
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E r = q/Aire^r 2 . Hence 

(5. 

Notice that Eq oc 1/r dominates over E r oc 1/r 2 at large r. For cAt <C vt, E « EqQ where 
the field is kinked, and so the electric field is perpendicular to the direction of propagation of 
the kink which is in the r direction. This is characteristic of electromagnetic waves which have 

k • E = 0. 

There will also be an azimuthal magnetic field in the spherical shell corresponding to c(t — At) < 
r < ct that drops off like oc 1/r and is due to the current density associated with the motion 
of the charge which before t — 0 had v = 0 and after t — At had velocity v = aAtz. In fact, 
within the kink region, 

E = B — —<j> 

4'K£§rC L c 

Since E = E6 the wave is linearly polarised in the 6 direction, i.e. the plane containing the 
acceleration direction (z) and the wave vector’s direction (?). 


(5.3) 


Eq = a^sinO 


Atysqt 2 


aq sin 9 
4:7T£o rc 2 


5.1.1 Total power 


The energy flux is given by the Poynting vector which for EM waves in empty space is given by 


S 


£<) Eq r 


£0 


aq sin 9 
4nsorc 2 


2 

r 


a 2 q 2 sin 2 9 ^ 
(4n) 2 £o r 2 c 3 ’ 


(5.4) 


Note the dependence on a 2 , g 2 , sin 2 9 and the 1/r 2 inverse square law dependence. The solid 
angle subtended at the origin by d(area) perpendicular to r is dQ = d(area)/r 2 . Hence, the 
power emitted during the acceleration per unit solid angle is 


dP a 2 q 2 sin 2 9 
dQ (4n) 2 £o c 3 
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The total instantaneous radiation power is 


P = 


f dP [ 1 a 2 q 2 sin 2 6 

= (b-^dn = / ,, , 2 o 




0 2 „2 


(5.6) 


which is Larmor’s radiation formula, named after British physicist Sir Joseph Larmor (1857-1942) 
who derived it in 1897, and is valid for v <C c. The restriction on v is not usually a problem as 
one can often Lorentz transform to a frame in which v = 0. 


Interestingly, this formula was derived using essentially this method by British physicist Sir 
Joseph J. Thomson (1856-1940) in (Thomson 1904, page 58) at a time when it was still thought 
EM waves propagated through a luminiferous aether! Thomson was actually explaining the pro¬ 
duction of “Rontgen rays” (what we now call X-rays) when “cathode rays” (energetic electrons) 
impinge on a metal target by a process we now describe as “bremsstrahlung” (meaning braking 
radiation). In that case the charges (electrons) initially had high velocities and were abruptly 
decelerated producing the EM radiation (X-rays). [Thomson, J J, 1904 Electricity and Mat¬ 
ter , page 58, Archibald Constable & Co., Westminster. Available at California Digital Library 
http: / / archive.org/details / electricityandmaOOthomiala]. 
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5.2 Electric dipole radiation 


Consider a charge q undergoing simple harmonic motion in the ^-direction between z — —d and 
z = +d. Its electric dipole moment is 

p (t) = qdcos(cot)z = po cos(cot)z. (5.7) 


The charge q experiences acceleration 

(5.8) 

(5.9) 


a (t) = —duo 2 cos(c ot) z, 

q 2 [a(t )} 2 = (qd) 2 uo 4 cos 2 (uot) = Pq uo 4 cos 2 (cot ). 


Using Larmor’s formula (Eq. 5.6) we obtain the total instantaneous and time-averaged electric 
dipole radiation power 


m = 


6^0 c 3 


cos 2 (c<;t), ( P) — 


Po u 4 
12tT£()C 3 


(5.10) 


5.2.1 Dipole radiation pattern 


Using Eqs. 5.5 and 5.9, the time-averaged power emitted per unit solid angle is 

(5. 

To illustrate the shape of the dipole radiation pattern, sin 2 6 is plotted against 6 in Fig. 5.3(a). 
To better visualise the radiation pattern the same curve is plotted with polar coordinates in 
Fig. 5.3(b) where the vertical axis is the z-axis and corresponds to the direction of the oscillating 
dipole, i.e. charge is being accelerated in the ±z direction. 

Note that there is no radiation in the direction in which charge is accelerated, and the radi¬ 
ation power per unit solid angle is maximum in all directions perpendicular to the direction 
of acceleration. As there is no dependence on (j) there is azimuthal symmetry, and the dipole 
radiation pattern could be better visualised as a 3D surface traced out by rotating the curve in 
Fig. 5.3(b) around the z-axis. 


(—) 

\d(l/ 


Pq cj 4 sim 6 

327r 2 £o c 3 


= ~(P)s in 2 e. 
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0 30 60 90 120 150 180 

6 (degrees) 


Figure 5.3: The power emitted per unit solid angle is proportional to sin 2 0 which is plotted 
against 6 in (a) with standard cartesian axes, and (b) polar axes. 

5.2.2 Short centre-fed antennas 

If an AC source I(t) = Iq cos uot is connected to a short (i.e. total length d « A) centre-fed 
dipole antenna (Fig. 5.4) it becomes an oscillating physical electric dipole, and so will radiate 
electromagnetic waves at frequency v — The current will be maximum at the centre 

(z = 0) of the antenna, and zero at its ends (z = ±d/ 2) because it has nowhere to go. If the 
antenna is short then the current will decrease linearly between the centre and the two ends, 
and the charge in each arm will be uniformly distributed along its length. 
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Z A 



d/2- 


0 - 


-d! 2- 


Figure 5.4: Short centre-fed dipole antenna. 


The charges in the upper and lower arms are given by 


Q 


/j 

I 0 cos (ut)dt = — sin (ut); Qi ow (t) = -Q up (t), 

id 


(5.12) 


2 J 0 2 1 0 

so that the line charge density on the upper and lower arms is H-- sin (out) and-- sin(cdt), 

ddd cod 

respectively, such that the dipole moment is 
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p (t) = posm(wt)z, po 


dlo 
2 w' 


(5.13) 


The total time-averaged radiation power is 

<p> - S- </2> ^ <5 ' i4> 

where (I 2 ) = /q/2 and, by analogy with resistance R where I 2 R gives the power dissipated as 
heat, 

Riad = 24 ^? ~ 197 (a) (f2) (5-15) 

is the “radiation resistance” of the antenna (A is wavelength) such that (I 2 )R Ta d gives the power 
radiated in EM waves. 


5.2.3 Half-wave antennas 


The half-wave antenna, i.e. full length d — A/2, has the advantage of having nearly zero re¬ 
actance, and so is an efficient radiator, as well as having a radiation pattern which closely 
resembles the sin 2 6 dipole radiation pattern. It also has a radiation resistance of R ra d = 7312 
which matches well the 7512 impedance of coaxial cable commonly used in radio-frequency (RF) 
applications. We shall discuss long antennas in more detail at the end of this chapter. 

A related type of antenna, which is effectively a dipole antenna, is the ground-plane vertical 
monopole antenna (Fig. 5.5a). It comprises a single vertical conductor insulated from the ground 
plane. The signal is connected between the base of the vertical and the ground plane near the 
base. The electric field of a charge placed above a grounded plane is identical (above the plane) 
to that of the charge plus a negative image charge at the mirror-image point below the plane. In 
a similar way, the ground-plane vertical monopole antenna acts like a dipole antenna. Usually 
a A/4 ground-plane vertical antenna is used as it is equivalent to a A/2 dipole antenna, at half 
the cost! The ground-plane need not necessarily be literally the ground. In fact the earth is 
far from being a perfect conductor and is often supplemented by laying a conducting surface 
such as wire netting on the ground. Car antennas are often A/4 ground-plane vertical monopole 
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antennas, with the “ground plane” beeng the car roof. 

To transmit or receive at microwave (GHz) frequencies, a A/4 ground-plane vertical monopole 
antenna is placed inside a microwave feed-horn located at the focus of a parabolic “dish”, and 
pointing to “illuminate” the dish (Fig. 5.5b). Such an arrangement with an appropriate “back¬ 
end” (receiver/transmitter and amplifiers) can be used to transmit and/or receive microwave 
signals. It is also the basic design of a radio telescope. 


(a) (b) 



Figure 5.5: (a) Vertical ground-plane monopole antenna, (b) Microwave feed horn (at focus of 
parabolic reflector) with a A/4 ground-plane vertical monopole antenna fed by a coaxial cable. 
(Not to scale.) 


5.3 Scattering of EM radiation 


We can understand scattering of electromagnetic radiation in terms of acceleration of charges by 
the electric field of an EM wave, or by the electric field of an EM wave inducing an oscillating 
electric dipole in a particle, conductor or dielectric, in either case resulting in scattered (re¬ 
radiated) EM radiation. The strength of the scattering is determined by the total scattering 
cross section a (m 2 ). The angular distribution of scattered photons is determined by the 
differential cross section da/d£L (m 2 sr -1 ). 
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5.3.1 Scattering of EM radiation by free electrons 


In the presence of a monochromatic plane wave E = Eq cos(k • r — cot), an electron located at 
r = 0 experiences electric field E = Eq cos(cut) and acceleration 


, N —eE (t) e ^ , % . 9 . e 2 E% , „ . 

a (t) - =-Eocos(cut), (a 2 ) = - (5.16) 

m e m e 2 

As the electron oscillates at the frequency of the incident radiation, the scattered radiation will 
have the same frequency. Using Larmor’s formula we see that the electron emits radiation with 
time-averaged power 


(p) 


127rrag£oc 3 


(W). 


(5.17) 


The power radiated by the electron is taken from the energy flux of the incident monochromatic 
plane wave E = Eq cos(k • r — cot), which at r = 0 is 
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S = 6 qEqC cos 2 ( c <;£), (S) = c (W m 2 ). (5.18) 

It is almost as if the electron absorbed all the radiation hitting area equal to (. P ) / (S) and re¬ 
radiated it with power (. P ). This area is called the Thomson cross section after J.J. Thomson 


Cj 1 — 


(_p± 

(S) 




= 6.65 x 10 


-29 


(m 2 


(5.19) 


and is the cross section for scattering of electromagnetic radiation by free electrons, and applies 
provided hu <C m e c 2 . 


5.3.2 Scattering of EM radiation by molecules 


A molecule acquires an electric dipole moment in the presence of an electric field p = <a m E. 
For the case of a dilute gas the Clausius-Mossotti equation yields 


a™ = — (£r- 1) = —{n 2 — 1) (5.20) 

n v n v 

where n v is the number density of molecules, e r the relative permittivity and n is the refractive 
index. In the presence of an oscillating electric field due to an incident monochromatic plane 
wave E 0 cos(k • r — cut) the molecule becomes an oscillating dipole which radiates with time- 
averaged power 


(p) 


12neo c 3 


£o \ 2 (n 2 - 1 ) 2 E ^ 4 
n v ) 127t^oc 3 


(5.21) 


Dividing ( P) by the time-averaged energy flux of the incident wave ( S) = we obtain 

the cross section 

_ (n 2 - l ) 2 4 _ (n 2 - l ) 2 /27tz/\ 4 _ 8n 3 (n 2 - l ) 2 4 

mo1 67 m^c 4 67 m^ \ c ) 3 n% 

Note that short wavelength light (blue) is scattered much more strongly than longer wavelengths 
(red), and this explains why the sky is blue and sunrises and sunsets are red. 
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5.4 Formal treatment of radiation by charges and currents 


So far I have discussed in a reasonable amount of detail the propagation of electromagnetic 
waves, but their generation has only been dealt with in an approximate fashion when we derived 
Larmor’s formula in a heuristic manner. A more rigorous approach involves deriving and solving 
the wave equations for the potentials V(r , £) and A(r,£) which include the charges p(r,£) and 
currents J(r, t) which are their sources, i.e. the inhomogeneous wave equations for the potentials. 


5.4.1 Inhomogeneous wave equations for the potentials 


These equations are obtained from Maxwell’s equations by substituting E 
into Gauss’ law 


V- 


-W 


dA 

dt 


P_ 

5 

£0 


vV + 


<9V • A 

dt 


P_ 

5 

£0 


- W - dA/dt 


(5.23) 


and into the modified Ampere’s law for B = V x A 


V x (V x A) 


Mo J + M o^o 


d_ 

dt 



V(V • A) — V 2 A 


ho J - V 



<9 2 A 


v 2 a 


d 2 A 


—/UqJ — V 


(v • A + /i q£q 



(5.24) 

(5.25) 

(5.26) 


Equations 5.23 and 5.26 can be written in a convenient form if we use the Lorentz gauge , i.e. if 
we set 


VT A 9V 

V A = —g,Q£Q — . 

dt 


Then, we get the inhomogeneous wave equations 


V 2 E 

v 2 a 


d 2 V 

- two -pr 

d 2 A 


p_ 

5 

£0 


-MoJ- 


(5.27) 


(5.28) 

(5.29) 
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5.4.2 Retarded potentials 


Imagine there is a variable point charge q(t) at the origin. Strictly, this is not allowed by 
conservation of charge, but if it were we would expect the potential outside the source to be a 
spherical wave solution of the homogeneous wave equation for V. Since q is at the origin, the 
solutions would be of the form V(r, t) = g(t + r/c)/r (incoming wave) and F(r, £) = g(t — r/c)/r 
(outgoing wave), where g(x) is any differentiable function, as may be verified by substitution. 
For unbounded space only the outgoing wave is the physical solution. Comparing this solution 
with the potential of a point charge q at the origin, the solution of the inhomogeneous wave 
equation for V must be 


V(r,t) 


g(* ~ r/c) 
4:7T£Qr 


(5.30) 


This makes sense as it takes time r/c for information about changes in the value of q at the 
origin to reach radius r. We can immediately extend these ideas to obtain the solution to the 
inhomogeneous wave equations (Eqs. 5.29). 
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For time dependent sources p(r,t) and J(r, £) we need to take account of the time taken for 
information to propagate from the source points r' to the field point r, which in vacuum takes 
a time |r — r'|/c. This means that the fields at the field point at time t depend on the charge 
density and the current density at r' at the retarded time t Te t = t— \r — r'\/c. The solutions of 
the inhomogeneous wave equations, the so-called retarded potentials , are then 


V{T,t) 


1 f p(/,t-\r-r'\/c) 3 , 

47T£ 0 J |r — r'| 


_ Mo f J(r ; ,t - |r - r'|/c) 
47 t J |r — r'| 


(5.31) 


5.4.3 Sinusoidal time dependence 

For sinusoidal time dependence of the sources 

p(r',f ret ) = p(r') e ~ iuitTet = p(r') e -MHr-r'|/c) = p ( r ') e ifc|r-r'| e -iurt ( 5 . 32 ) 

J(r',t ret ) = J(r / )e ifc l r - r, l e -^ t . (5.33) 

If we are only interested in the radiation fields we do not need V(r ,£) as we can obtain E(r,t) 
from B(r,t), and B(r,t) from A(r,t). 

The vector potential for sinusoidal sources is 


A < r ’ f > = £ 


«> r 


r — r 


(5.34) 


Approximations may be made at both small and large different distances from the source, but 
at intermediate distances no approximations are possible. 


In the near zone where k\r — r'\ <C 1 we have the quasi-static limit where 




•Hr') 


r — r 


d 3 r'e~ iujt . 


(5.35) 


In the far zone or radiation zone , where r' <C r, referring to Fig. 5.6 it is clear that the error is 
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negligible if we replace |r — r'| 


D ikr 




by (r — r • r'). Then 

e ~ikr r' ^ g -iart 


(5.36) 


0 


Figure 5.6: Justification of the approximation |r — r'| ~ (r — r • r') when r r'. 



5.5 Antenna theory and multipole expansion in the far zone 


We can make a multipole expansion for e lkr ' r ' using a Taylor series expansion 


^-ikr-r' 


E 

m =0 


(-ikr -r') m 
ml 


1 — ikr r' — ... 


(5.37) 


The 1st term in the series corresponds to electric dipole radiation and the 2nd term corresponds 
to magnetic dipole radiation. Physical radiating systems for these two cases are the short centre- 
fed antenna and magnetic dipole antenna shown in Fig. 5.7. 



Figure 5.7: (a) A short centre-fed antenna 


has d A. (b) A magnetic dipole antenna. 
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5.5.1 Electric dipole radiation term 


The vector potential for this term is 


A(r ,i) 



(5.38) 


and if we knew the velocities of all the N individual charges making up the current we could 
replace the integral with a sum 



N 

2=1 


N 


IN 


drj 

dt 


d 

dt 


N 

2=1 


dp 

dt 


(5.39) 


where p is the electric dipole moment. Now, we can write for the time dependence of the dipole 
moment p (t) = p^e~ lujt z such that dp/dt = —i(jopoe- l(Jjt z, giving the vector potential for an 
oscillating electric dipole 
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A(r, t) 


,/J-ou 
1 4t r 


Po 


gi(kr-ujt) 


r 


z 


1 4tt 


Po 


gi(kr-ujt) 

-(cos0r 

r 


sin 66). 


(5.40) 

(5.41) 


The magnetic field is V x A and for periodic radiation fields Vx —> zkx where, for the 
spherical wave, k = kr. Then 


B(r, i) = 
B(r,i) = 
H(r ,*) = 


fiQ uj e 


i(kr—ujt ) 


ikr x (cos 9 r — sin 9 6 ), 


flokuj e^kr-ut) 

■—;—Po-sin 6/0. 

47T r 

kuj e ^r-ujt) 


--J-P0 - 

47t r 


sin# 0. 


(5.42) 

(5.43) 

(5.44) 


We can obtain the electric field from Ampere’s law for a current free region, V x B = 
po£o(dE/dt), which is equivalent to zk x B = —icjpo^oE for periodic radiation fields. Hence, 


E(r,t) 


( — ) 

\UJ/J, 0 s 0 j 


( (M)koj \ 

V 47r ) 


g i{kr—uot ) ^ 

Po-sin#r x 0, 



gi i(kr—ujt ) ^ 

po-sin 9 6 . 

r 


(5.45) 

(5.46) 


Thus, for an electric dipole antenna as shown in Fig. 5.7(a) with the antenna oriented in the 
z (vertical) direction the radiation will be polarised with E in the 0 direction, i.e. vertical if 
viewed from a horizontal plane containing the antenna. 


The Poynting vector is S = Re{E} x Re{H} 


i(M) =(££) p§ 


cos 2 (kr — ujt) o ^ 

-\—o-- sin 9 r, 

4 nr z 


,4^2 


U cos 2 (kr — ujt) sin 2 9 r. 
Ioitsoc 6 ^ 


(5.47) 

(5.48) 
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The time-averaged power radiated per unit solid angle is r 2 (*S), so that 


dP u A pl . 2 

W = sw sm 9 


in agreement with the result in Eq. 5.11 obtained from a heuristic treatment. 


(5.49) 


5.5.2 Magnetic dipole radiation term 


A real example of an oscillating magnetic dipole is the magnetic dipole antenna which is sketched 
in Fig. 5.7(b). The magnetic dipole term in the expansion for the vector potential is 


A(r, £) 


—ik 


V oe 


ikr 


47 rr 


J (Vr')J(r') 


—iujt t3 / 


dV, 


(5.50) 


and in Section 5.6 of “Essential Electromagnetism” I showed that 


J J(r') r • r'd 3 r = m x r 


where 



x J(r') d 3 r' 


(5.51) 


(5.52) 


is the magnetic dipole moment. Hence, setting m = m^z = mo(cos0r — sin 6 6 ) the vector 
potential of an oscillating magnetic dipole moment is 


A(r,t) ik—r x (m^cosOr — rriQsmO 0) e %ujt 
47 rr 

_ _ Mo^o sin 9 


(5.53) 

(5.54) 


The magnetic field is V x A = zk x A for periodic radiation fields, where k = k r, i.e. 

B(r, t) = i(kr) x A(r, t) = -k 2 ^°T amd (5.55) 

47rr 

As we did for electric dipole radiation, we can obtain the electric field from Ampere’s law which, 
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for periodic radiation fields, is equivalent to E = —k x B/cj/io^o- Hence, the electric field is 
k 

E(r,i) = - rxB(r,i), 

Cdfl()£0 

k 3 nom 0 smd ikr _ iut ~ 

A ^ t 5 

(jj/ioeo 47 rr 

fc 2 m 0 sin jkr-iojt ~2 
ceo 47rr ^ 

Thus, for a magnetic dipole antenna as shown in Fig. 5.7(b) with the loop antenna located in 
the horizontal plane the radiation will be polarised with E horizontal. 

The Poynting vector is S = Re{E} x Re{H} = Re{E} x Re{B}/^o 

S(r, t) = — cos 2 (kr — ut) sin 2 dr. (5.59) 

16n z c 6 r z 

The time-averaged power radiated per unit solid angle is r 2 (*S), so that 


(5.56) 

(5.57) 

(5.58) 
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( 



4 2 

LJ IIqUIq 

327r 2 c 3 


sin 2 6 


(5.60) 


which has the same form as for electric dipole radiation with the substitution of /io^o for Po/^o, 
and has the same sin 2 6 dipole radiation pattern. The total time-averaged power is 


(P) = 


(j^/iorriQ sin 2 a 
327r 2 c 3 


sin 2 9dQ, 


. sin 2 a 2 

— 2?r * 11 


£ 


327r 2 c 3 


(1 — cos 2 6)d cos 6, 


4 2 

/ipiirmg 

127TC 3 


(5.61) 

(5.62) 

(5.63) 


5.6 Radiation from quadrupole antennas 


Magnetic quadrupole radiation may be produced by a pair of oppositely directed magnetic 
dipole antennas, as sketched in Fig. 5.8. We may add the vector potentials of the two magnetic 
dipoles taking account of the angle-dependent phase shift for each dipole. Looking at the paths 
from the centres of the two AC current loops which comprise the radiating system, the wave 
generated by the upper dipole travels a distance shorter by 5/2 cos 0, and the wave generated by 
the lower dipole travels a distance 5/2 cos 6 longer, compared to (fictitious) waves from midway 
between the centres of the current loops. 



Figure 5.8: Far zone approximation applied to a physical magnetic quadrupole. 

We can thus adapt the magnetic field of the magnetic dipole antenna (Eq. 5.55) to obtain the 
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magnetic field of a magnetic quadrupole 


g/p _ _ j„2 sift 0 ikr—iojt / — ikb cos 6/2 _ ^ikb cos 0/2\ 'q 

' ’ ~ 47rr V ) ’ 

2 /igmo sin 0 %kr—iojt t 


= -k 


4i rr 

zfc 2 /iomo e ikr ~ iujt 
2nr 


2isin (kb cos 6/2) 6 , 
sin0 sin (kb cos 6/2) 6. 


The electric field, Poynting vector and radiation pattern are then 


EM) = 

SM) = 

/—\ 
XdSi/ 


ifc 2 mo e' 


ikr—iut 


sin 6 sin (kb cos 6/2) (/>, 


So c27rr 

fc 4 mg[-sin 2 (fcr-w0] . 2 « . 2 


sq c47r 2 r 2 


sim 0 sim (kb cos 0/2) r, 


fc 4 mg . 2 


So c8n 2 


sin 2 6 sin 2 (fc6 cos 0/2) . 


(5.64) 

(5.65) 

(5.66) 


(5.67) 

(5.68) 

(5.69) 


The radiation pattern for such an oscillating magnetic quadrupole is shown in Fig. 5.9 for three 
values of b. For b — A/2 we have the equivalent of a half-wave antenna and we see a fairly 
standard quadrupole radiation pattern having two lobes at 50° and 130°. For larger b values 
the radiation pattern becomes more complicated with the addition of “side-lobes”. Electric 
quadrupole radiation can be dealt with by similar methods, i.e. the superposition of a pair of 
oppositely pointing oscillating electric dipoles. 


5.7 Radiation from long antennas 


In Section 2.2 we discussed standing waves on open-circuited transmission lines (Fig. 2.7). 
Imagine pulling apart the final quarter-wavelength of the line and stretching the ends apart 
to become an antenna of length equal to half a wavelength, i.e. a “half-wave antenna” as in 
Fig. 5.10. The current distribution on each line will be essentially unchanged, and so there will 
be a standing sinusoidal wave on the antenna as shown. The only justification of this procedure 
is that for these current distributions on the antenna one can predict the radiation pattern 
which is in agreement with that observed. 
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Figure 5.9: The radiation pattern of an oscillating magnetic quadrupole is shown for (a) 
b = A/2, (b) b = 3A/2, (c) b = 3A. The dotted curve in part (a) gives a dipole radiation 
pattern. 


For an antenna of arbitrary full-length d the current along the antenna is 




I 0 sin [k {d/2 + z)] e~ iut {-d/2 < z < 0) 


I 0 sin [k {d/2 - z)] e~ iut {0 < z < d/2) 


(5.70) 


where 2 is the distance along the antenna from its centre. We can use the far zone approximation 
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(Eq. 5.36) to obtain the the vector potential, 

(5.71) 

(5.72) 

where we have used r • r' = z'cos6 (see Fig. 5.11). 


A(r ,t) = 


Moe 


ikr rd/2 


Anr 


rd/2 

J-d/2 


I(z\ t)e 


—ikz 'cos 0 


dz , 



Figure 5.11: A centre-fed antenna. 
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Then, 


J—d/2 


4-7T r 

rd/2 


e -ikz’cose dz f + 


ra/Z 

/ sin [k (d/2 — 2 /)] 

Jo 


e -ikz'cos 0 dz ' 


(5.73) 


= —z 


- 


i(kr—ujt) 


47rr 


-ikz' cose ( cos i k ( d / 2 + z ')] + icos(9sin[/c (d/2 + z')] 


k sin 2 0 


0 

J —i/2 


+ 


= — (r cos 6 — 6 sin 6) 


-ikz'cose ( ~ cos[fc (d/2 - z')\ + i cos 6 sin[fc (d/2 - z:') } 

k sin 2 6 

cos(kd/2 ) — cos(cos0 kd/2) 


f\] \ 1 V 2 


J 0 


27rr 


k sin 2 0 


, (5.74) 


(5.75) 


For radiation fields the magnetic field is ik x A where, for the spherical wave, k = k r. Then 


B(r,£) = sin 0 0 
H(r,£) = sin 0 0 


lioIoe l ( kr cos(fcd/2) — cos(cos0 kd/2) 
2nr sin 2 0 

I oe i(k r ~ujt) 2 ) — cos(cos0 kd/2) 


2i rr sin 2 0 

and the electric field is E = —k x B/cj/io£o> giving 


E(r, t) — k sin0 6 


Ioe l ( kr cos(kd/2) — cos(cos0 kd/2) 


27T£Q(jJr 


sin 2 0 


(5.76) 

(5.77) 


(5.78) 


Finally the Poynting vector S = ExH and radiation pattern (dP/dQ) = r 2 (|S|) are 


S(r,i) 


I$k cos 2 (kr — cot) [cos(fcd/2) — cos(cos0 kd/2)] 2 
47r 2 £o (jor 2 r 2 sin 2 0 


(5.79) 


dP\ (I 2 ) [cos(fcd/2) — cos(cos0 kd/2)} 2 
d£l / 4n 2 £oc sin 2 0 


(5.80) 
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Current distributions and radiation patterns are given for standing-wave antennas for various 
lengths in Fig. 5.12. Of particular interest is the half-wave antenna shown in Fig. 5.12(a) whose 
radiation pattern can be approximated by that of an oscillating point electric dipole (dotted 
curve). As we move to longer dipoles (or shorter wavelengths) the antenna beam pattern 
becomes more complicated and no longer resembles a dipole pattern, having only a weak lobe 
at 6 = 0 but having strong side-lobes. 


(a) current 





Figure 5.12: Standing wave current distributions (left) and antenna radiation patterns (right) 
are shown for antennas of various lengths. The current distributions are shown at t = 0 (red) 
and at half a period later t = T /2 (blue dashed). The radiation patterns are normalised to 
the maximum values, (a) d = A/2 is the so-called half-wave antenna whose radiation pattern 
resembles the sin 2 6 dipole radiation pattern (dotted curve); (b) d = 3A/2. 


Integrating Eq. 5.80 over solid angle gives the time-averaged power 


(p) = 


(I 2 ) f 1 [cos(kd/2) — cos(cos6 kd/2)Y 


47 T 2 £qC 


2n 


i: 


1 — cos 2 6 


d( cos#), 


(5.81) 


where we have used (cos 2 (fcr — uot)) = 1/2 and ([ I(t )] 2 ) = ig/2. This is in the form of Ohm’s law 
which, for an AC current, is ( P) = ( I 2 )R , and so the radiation resistance of a long centre-fed 
antenna is 


A!rad — 


1 


^ [cos(fcd/2) — cos(cos# kd/2)Y 


2n£oc J_i 


1 — cos 2 6 


d(cos6). 


(5.82) 


This is plotted in Fig. 5.13(b)-(c) as the solid red curves. 
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(a) (b) 




Figure 5.13: (a) Radiation resistance versus antenna length calculated using the formula for a 
long antenna (solid red curve, Eq. 5.82) and for the approximate formula for the short centre- 
fed antenna (dotted blue curve, Eq. 5.15); (b) Current distribution on a centre-fed antenna of 
length d = 0.05A (blue), d = 0.5A (red) and d = 0.85A (green). 


The radiation resistance of a short centre-fed dipole antenna (Eq. 5.15) has been added to 
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Fig. 5.13a) as the dotted blue curves, and is very different to that calculated using Eq. 5.82 


(red curves). The reason for this is in the use of a different definition of the current appearing 


in Ohm’s law. For the case of a short center-fed antenna the maximum current / max is much 
less than Iq which is the peak of the assumed sinusoidal standing wave — compare the peak 
current / max <C Iq for d = 0.05A with the peak currents / max = Iq for d — 0.5A and d = 0.85A 
shown in Fig. 5.13(b). 

The antenna design will depend on the application, but to maximise the power radiated its 
impedance Z = R + iX should ideally be purely resistive, i.e. the reactance X should be 
negligible. The reactance X = ujL — 1/ujC of an antenna depends on its dimensions and the 
frequency, and is generally capacitive for short antennas and inductive for long antennas. In the 
case where the reactance is non-zero the antenna can be made to radiate more inefficiently by 
using impedance-matching techniques. Antenna design is an important topic in its own right 
and contains many subtleties which we cannot go into in this book — here I have only skimmed 
the surface. 

Summary of important concepts and equations 

An accelerated charge radiates: Larmor’s formula 

— Instantaneous power 


dP a 2 q 2 sin 2 6 
(47t) 2 £oC 3 ’ 



Time-averaged electric dipole radiation by oscillating dipole p (t) = po cos(u;t)z 



— Dipole radiation pattern has (dP/dQ) oc sin 2 6. 

— A short centre-fed dipole antenna (total length d « A) fed with current I = Iq cos (ujt) 
has dipole moment 



(5.83) 
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— For a short centre-fed dipole antenna resistance is 


(p) = 


d 2 (I 2 )(jj 2 

24?reoc 3 


— (/ )-Rrad! -Rrad — 


d 2 U 2 
24 tT£qC 3 


197 'a 


(Si). 


Half-wave antennas with d — A/2 are more efficient radiators (i? ra d = 73.13 fl) and have 
radiation patterns almost identical to the standard dipole radiation pattern. 


Scattering 


— EM waves incident on a particle (or small object) will produce oscillating charges or 
currents causing some of the incident radiation to be scattered (i.e. absorbed and re¬ 
radiated) . 

— The total cross section for scattering of EM waves is the time-averaged power re-radiated 
divided by the time-averaged incident energy flux. 
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— For free electrons an EM wave causes acceleration —eE (t)/m e and re-radiation according 
to Larmor’s formula. 

— The classical formula for scattering by free electrons is the Thomson cross section 


(JT — 


(5) 


67rmg£QC 3 


6.65 x 10- 29 (m 2 ). 


— For molecules, the EM wave induces an oscillating molecular dipole a m E(£) so that the 
scattering cross section is 


^mol 


(n 2 — l) 2 
67rn^c 4 


U) 


(n 2 — l) 2 / 2 - kv \ 4 

6nn% \ c J 


87r 3 (n 2 — l) 2 
3 n 2 v 


A” 4 . 


Inhomogeneous wave equations for the potentials 


V 2 V 


d 2 V 


£_ 

5 

£0 


V 2 A 


d 2 A 


~Mo J 


(5.84) 


Retarded potential solutions 


1 [ I r_ r'|/c) 3 , 


^(r,t) = t~~ / d 6 r', 

4tt£ 0 J |r — r'| 

Mo f J( r 7 ,t- |r — r'l/c) 




r — r' 


The vector potential for sinusoidal sources is 


^ = rj 


Mo /• J(rO e^lr-r'l 


d?r'e~ iut . 


r — r' 


(5.85) 

(5.86) 


(5.87) 


— In the far zone or radiation zone where / < r we can replace |r — r'| by r — r ■ r' 


A(r, t) = 


Moe 


,ikr 


47rr 


/ J (r') 


e -ikr-r’ d 3 r ' e -iut' 


(5.88) 
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— For a long antenna along the z axis 


A(r ,t) 



(5.89) 


Exercises on Chapter 5 


5-1 A walkie-talkie has a short centre-fed dipole antenna 12 cm long and radiates 3 W at 
450 MHz. Assume the antenna has been balanced by a suitable inductor such that its 
impedance is purely resistive, and that the transmission line from the transceiver is matched 
to the impedance of the antenna. 

(a) What is the radiation resistance of the antenna? 

(b) What is the value of the peak current in the antenna? 

(c) What is the peak electric field at the location of a receiving antenna at distance r? 

(d) With the transmitting walkie-talkie held vertically, an identical receiving walkie-talkie 
also held vertically at a distance of 2 km away can barely receive the transmission, i.e. 
the range is 2 km. What would be the range if: (i) the transmitting antenna were tilted 
at 60° towards the receiving antenna, (ii) the transmitting antenna was vertical but the 
receiving antenna was tilted at 60° towards the transmitting antenna, (iii) the transmitting 
antenna was vertical but the receiving antenna was tilted at 60° to the vertical in a plane 
perpendicular to the line of sight to the transmitting walkie-talkie? 

5-2 Starting with the Clausius-Mossotti formula and the time-averaged dipole radiation power 


= -{e r - 1), (P) 

n v 


127T£oc 3 ’ 


(5.90) 


find the cross section for scattering of light at A 500 x 500 nm wavelength by air molecules, 
and its associated mean free path l/(n v a). Dry air at STP has p = 1.30 kg m -3 , and 
n = 1.00029, and its molecular weight is 29. 

5-3 In the presence of an applied uniform electric field E = Eqz the radial component of the 
electric field outside a perfectly-conducting sphere of radius a centred at the origin becomes 


E r (r,d,<t>) = E 0 



COS#. 


(5.91) 


Find the cross section for scattering of monochromatic EM waves by a perfectly-conducting 
sphere of radius: (a) a « A, and (b) a^> A. Assume electric dipole radiation is responsible. 
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5-4 In the presence of an applied uniform magnetic field B = Bqz a perfectly-conducting sphere 
(radius a) acquires a magnetic dipole moment mo. 

(a) What are the boundary conditions at the surface of the sphere? Find the magnetic 
dipole moment that, together with the applied magnetic field, satisfies the boundary con¬ 
ditions. 

(b) Find the contribution of magnetic dipole radiation to the cross section for scattering 
of monochromatic EM waves by a perfectly-conducting sphere of radius A. 


5-5 In bremsstrahlung an energetic electron is accelerated in the electric field of an atomic 
nucleus and emits electromagnetic radiation. Consider the case of an electron of kinetic 
energy 10 keV travelling on a trajectory parallel to the z-axis, and that there is a lead 
(Pb) nucleus located at impact parameter b = 10 -10 m (distance of closest approach to 
the initial trajectory). Assume that for this impact parameter the electron’s deflection is 
small, and that a calculation of the energy lost to electromagnetic radiation can be done 
classically. 

(a) By comparing b with the de Broglie wavelength h/mv , reassure yourself that a classical 
calculation is valid for this impact parameter and velocity. 
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(b) By checking whether or not the deflection angle is small, reassure yourself that approx¬ 
imating the trajectory by a straight-line is valid. 

(c) Estimate the total energy radiated by the electron. [Neglect atomic electrons.] 

5-6 A pulsar is a rapidly-spinning magnetised neutron star observed at radio and sometimes 
higher (optical, X-ray and gamma-ray) frequencies. Neutron stars have masses of typically 
M ^1.4 Mq (solar mass 1 M Q w 2 x 10 30 kg), radius R ~ 10 km, moment of inertia 
ImoI — \MR 2 ~ 10 38 kg m 2 and a wide range of surface magnetic fields. Outside the 
neutron star the magnetic field may be approximated by that of a magnetic dipole at its 
centre making angle a to the spin axis. Pulsars typically emit one or more narrow pulses 
at the same time into each spin period P, and measurements at different epochs generally 
show that the spin is slowing down as characterised by the time-derivative of the period P. 

(a) Derive formulae in terms of P and P for the rotational kinetic energy P ro t, the rate 
of loss of rotational kinetic energy E rot , the pulsar’s characteristic age r = E rot /E rot and, 
assuming the slow-down is due to conversion of rotational kinetic energy to magnetic dipole 
radiation, the pulsar’s minimum equatorial surface magnetic field P m in- 

(b) The radio pulsar PSR J0157+6212 has P = 2.355 s and P = 1.89 x 10 -13 (“The 
Australia Telescope National Facility Pulsar Catalogue”, Manchester, R.N., et ah, 2005, 
Astron. J., 129, 1993, and http://www.atnf.csiro.au/research/pulsar/psrcat/). Find E r G t, 
Prot i T — Prot/Prot ^md P m i n - 
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A SI units and dimensions 


The SI base units are: the metre (m) which is the unit of length, the kilogram (kg) which is the 
unit of mass, the amp (A) which is the unit of electrical current, the second (s) which is the 
unit of time, the kelvin (K) which is the unit of thermodynamic temperature, the mole (mol) 
which is the unit of the amount of substance, the candela (cd) which is the unit of luminous 
intensity. 

In any derivation it is important to check that the dimensions are the same on both sides of an 
equation. For this we will need the dimensions of electromagnetic quantities in terms of mass, 
distance, time and amps, or equivalently the units in terms of the SI base units. 

The units and dimensions of the major electromagnetic quantities are tabulated below. [See 
e.g. NIST Special Publication 330, http://www.nist.gov/pml/pubs/index.cfm.] 


Quantity 

Symbol 

Unit 

Dimensions 

work 

w 

joule (J) 

m 2 kg s -2 

power 

p 

watt (W) 

m 2 kg s -3 

electric charge 

Q 

coulomb (C) 

s A 

electrostatic potential 

V 

volt (V) 

m 2 kg s -3 A -1 

electric field 

E 

volt per metre (V m -1 ) 

m kg s -3 A -1 

electric displacement 

D 

coulomb per square metre (C m -2 ) 

m 2 s A 

electrical resistance 

R 

ohm (Q) 

m 2 kg s -3 A -2 

capacitance 

C 

farad (F) 

m -2 kg -1 s 4 A 2 

permittivity 

£ 

farad per metre (F m -1 ) 

m -3 kg -1 s 4 A 2 

magnetic field 

B 

tesla (T) 

kg s -2 A -1 

magnetic intensity 

H 

amp-metre (A m) 

m A 

magnetic flux 


weber (Wb) 

m 2 kg s -2 A -1 

inductance 

L 

henry (H) 

m 2 kg s -2 A -2 

permeability 


henrys per metre (H m -1 ) 

m kg s -2 A -2 


166 


Download free eBooks at bookboon.com 






Essential Electrodynamics 


Appendices 


B Derivation of the Maxwell stress tensor 


Momentum conservation in electrodynamics is discussed in Section 1.4. Here we derive the 
formula for the Maxwell stress tensor, i.e. T{j in Eq. 1.54, which must therefore satisfy 


N >( fe 

* \ 3 


ax , T ” 1 A 


1 


e 0 [(V • E)E-Ex(VxE)] + — [(V • B)B-Bx(VxB)] d A r. (B.l) 
l/l Mo 


Using index notation with the Einstein summation convention, 


[(V • E)E - E x (V x E)] • = EiVjEj - £ij k Ej£ k i m \7iE m (B.2) 

= EjS jEj EfcijSkimEj^iEm (B.3) 

— EiS7 jEj (fiilfijm ^im^jl)Ej^iE m (B.4) 

= EiS7 jEj — (E m V iE m — E{S7 \Ei) (B.5) 

— EiV jEj + EjX7 jEi — —V iEmEm (B.6) 

[(V • E)E - E x (V x E)] i = VjEiEj - 1 Vj-% E • E. (B.7) 

Similarly, [(V • B)B — B x (V x B)]^ = VjBiBj — iVjV/B • B, (B.8) 

and so 


s 0 [(V • E)E -Ex(VxE)] + — [(V • B)B - Bx(VxB)] 


Mo 


E giV 


i V j Sq 


EiEj + — —5ij(E • E + (?B • B) 


d 




EiEj (?BiBj — — Sij(E • E c 2 B • B) 
2 


(B.9) 

(B.10) 


Hence, 


T{j — sq 


EiEj -|- c?B % B 3 — — Sij (E • E c 2 B • B) 


(B.ll) 
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C Time-derivative, divergence and curl of the fields of a monochro¬ 
matic EM plane wave 


Derived here for Eoe^ k r ojt ' 1 using index notation; results for B(je !<k r ojt ' 1 have identical form. 


V • E = V* 


E 0 e HK ' 


,z(k-r —ujt) 


= \7l E 0:l e^ 


,z(k-r— ut) 


= E 0 j e i( k -r-^) Vj [ i(fcmrm _ wt )] 
= ik m E 0 )je i( k - r -^)( Vj r m ) 

= ik m E^e^ r -^8 lrn 

= ik m E 0 , m e i( - k - r -^ 

= ik ■ Eoe i( - k ' r ~^ 


V • E = tk • E 


(C.l) 
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[V X E] h = ShlrNl 


Eoe i( k -r-^) 

J n 

= e Wn V,^,, n e <(k - r - wt) 

= £wn-So,ne* (k ' r_a;t) Vi[i(k • r - ut)] 
= e w „^ ) ne i(k ' r - wt) V { [i(Avnr m - cut)] 
= ifc TO e Wn ^o,ne i(k ' r “ a;t) V i r m 

= ifc m£Wn E 0 ,ne i(k ' r “ a;t) ^ m 

= ie Wn ^£;o, n e i ( k - r -^ ) 


-]fe =»{l 


fV x Ek = i<kx 


E 0 e : 


i(k-r— ujt) 


]}, 


V x E = ik x E 


(C.2) 


^ T? ^ TT ^(k-r—< ujt) 

m E m Eoe 

= -iuEoe^-^ 


—E = —mo E 

dt 


(C.3) 


D EM field of a rectangular cross-section waveguide 


Eqs. 4.15-4.18 are used to find E(r, t) and B(r, t) from B®(x,y) (Eq. 4.36) for the TE modes, 
and from E®(x,y) (Eq. 4.31) for the TM modes. The following formulae based on Eq. 4.38 will 
be useful 


^mn 





(D.l) 
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For the TE modes m or n or both must be non-zero: 


E x = 


En. = 


E, = 


B x = 


B q , — 


l o c 




nn / mnx\ . fnny\ 

-y B mn cos — — J sin \^-y- J sin {k mn z - ut) , 


u c mn 


B mn sin 


o -‘-'mn 

^mn ® 


0, 

kmn(u)c 2 mu 


fmnx\ fnny\ 

m ^—J cos y — — J sin (k mn z - ut) , 


/ mnx\ (nny\ 

Bmn Sin (-) cos — J sin - c ut ), 


V a 


k m n(u)c 2 nn 

—— ^ 


^mn b 


(mnx\ (nny\ 

cos ^-J sin y—yj sm \k mn z — cut) , 


B z = 


fmnx\ fnny\ 

cos ^- J cos y~y~J cos \kmnZ — cut) . 

For the TM modes both m and n must be non-zero: 
k mn (u)c 2 mn 


E x = - 


' En 


Enj - 


E z = 


k mn (uj)c 2 nn 


^mn ^ 


— Em,n sin 


_ fmnx\ fnny\ 

Emn sm ^—J sm — j cos ( k mn z - ut ), 


B x = — 


u c nn _ 

9 T~ *^mn 
0 

,2 


Bn, — 


mn 

(jjc 2 mn 


Emm, sin 


9 J - v mn 

^mn ® 


(mnx \ (nny \ N 

m ^—J cos — J sm ( k mn z ~ ut ), 


(D.2) 

(D.3) 

(D.4) 

(D.5) 

(D.6) 

(D.7) 


B z = 0. 


sin ( k mn z - tot ), 

(D.8) 

1 (kmnZ lot) , 

(D.9) 


(D.10) 

nn z - ut), 

(D.ll) 

lZ - lot ), 

(D.12) 


(D.13) 
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E Summary of vector calculus identities 


E.l Integral theorems 


Gauss’ theorem: 


[ V x F • dS = l F • dr (L bounds S ) 
JS JL 


Corollary to Gauss’ Theorem: 


[ V x F d 3 r = - <f F x dS (S bounds V) 

Jv Js 


Stokes’ theorem: 


[ V • F d 3 r = <f F • dS (S bounds V) 

Jv Js 


(E.l) 


(E.2) 


(E.3) 
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E.2 Product rules 

V(itv) = u(Vv) + v(Vu) (E.4) 

V(A • B) = A x (V x B) + B x (V x A) + (A • V)B + (B • V)A (E.5) 

V • (A x B) = (V x A) • B - (V x B) • A (E.6) 

V • (uA) = (Vw) • A + u(V • A) (E.7) 

V x («A) = (V«) x A + «(V x A) (E.8) 

V x (A x B) = (V • B)A - (V • A)B + (B • V)A - (A • V)B (E.9) 

E.3 Second derivatives 

v • (VVO = v 2 ^ (E. 10 ) 

V x (VV>) = 0 (E.ll) 

V(V • A) = V x (V x A) + V 2 A (E.12) 

V • (V x A) = 0 (E.13) 

V x (V x A) = V(V • A) - V 2 A (E.14) 
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E.4 Vector operations in cartesian coordinates 


dr = dx x + dy y + dz z 
d 3 r = dx d?/ dz 

v/ = ai x + S y+ & z 

V • A = 


5a: 


5y 


V x A = 


5A Z 5 A 


dy 


dz 


dz 


x + 


5 A,- 

5z 


V 2 .av av aV 

5a; 2 dy 2 dz 2 


dA z \ ^ t fdAy 0A X \ ^ 

^J y+ Ur-^J Z 


(E.15) 

(E.16) 

(E.17) 

(E.18) 

(E.19) 

(E.20) 


E.5 Vector operations in spherical polar coordinates (r,0,<p) 


dr = dr r + r d6 6 + r sin 0 d<p<p 
d 3 r = r 2 dr sin 6 dO dxp 


„, 5/ ^ 15/- 1 5/ - 

^/ — "^77 # H •—77 aT 0 

5r r 5tf r sintf 5(p 


_ . 1 5 . 2 . . 1 5 . . . . 1 5A> 

^ ' A — ~2 ^-( r A r ) H 7 - 777, (sni 6Ag) + 7 - 

dr r sm# dO r siny dq> 


V x A = 


1 


r sind 


d 8 Aq 

M {sm<> A*) - 


r + 


1 


1 9A r _ d_ ' 
sin 6 d<p dr ^ 


e 


i 

+ — 
r 


5 

dr {TAl) 


dA r 

~W 


<t> 


V 2 /= — 0 ^ 1 + 


r<2 dr 


dr 


1 


d 


r 2 sin 6 d6 


— sind —- + 


df\ 1 d 2 f 
09 ) + r 2 sin 2 6 dPp 2 


(E.21) 

(E.22) 

(E.23) 

(E.24) 


(E.25) 

(E.26) 
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E.6 Vector operations in cylindrical coordinates (p,<j),z) 


dr = dp'p + pd(f> <j> + dzz 
d 3 r = pdpdcfi dz 

df ^ 1 df-2 df ^ 

J dp p d<j) V dz 


__ . Id, . 

^ • A — - ^“(P-^p) “I oT 
pop p dcp 


+ 


__ . .'IdA* 

V X ~~ ( p d<f> 


0A { , 


dz 


P + 


dAz 

dz 

dA c 


1 d 


v v- p.; 

p dp V <9p 


df 


dz 

j ^dV dV 
p 2 dcj ) 2 dz 2 


dA, 

dp 


-7 1 / d . dA p 

<# ’ + p\Bp ip At) ~^ 


(E.27) 

(E.28) 

(E.29) 

(E.30) 

(E.31) 

(E.32) 
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F Maxwell's Equations 

Differential form: 


Vacuum: 


Matter: 



V • E 

P_ 

£o 

V D 

= Pf 

(Gauss’ law) 

V B 

= 0 

V B 

= 0 

(no magnetic charge) 

V x E 

dB 

- ~dt 

V x E 

5B 

~ ~~dt 

(Faraday’s law) 

V x B 

, dE 

= no J + ho^o 

V x H 

, m 

- J/+ aF 

(modifiedAmpere’s law) 


(F.l) 


Integral form: 



Note: in the 1st row S bounds V; in the 2nd row S is any closed surface; in the 3rd and 4th 
rows L bounds S. 
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Index 

absorption, 83 
absorption coefficient, 86 
aether, luminiferous, 45 
Ampere’s law, modified, 22 
Ampere’s law, unmodified, 20 
amplitude reflection coefficient for parallel po¬ 
larisation, 67 

amplitude reflection coefficient for perpendic¬ 
ular polarisation, 65 

amplitude transmission coefficient for parallel 
polarisation, 67 

amplitude transmission coefficient for perpen¬ 
dicular polarisation, 65 
angular frequency, 35 
anomalous dispersion, 92, 93 
antenna theory, 146 
antenna, ground-plane, 139 
antenna, half-wave, 139 
antenna, short centre-fed, 137 
antennas, long, 152 

boundary conditions at interface, 62 
Brewster’s angle, 67 

cartesian coordinates, 173 
coaxial cable, 42, 116 
coherence, 50 
coherence area, 53 
coherence length, 50 
coherence time, 50 
coherence width, 51 
conductivity, dilute plasma, 99 
conductor, dispersion relation, 94 
conductor, refractive index, 94 
conductor, wave equation, 94 
continuity equation, 20 
cross section, Thomson, 142 


cut-off frequency, 102 
cylindrical coordinates, 174 

dilute plasmas, 97 

dilute plasmas, group velocity, 101 

dilute plasmas, phase velocity, 101 

dimensions, SI units, 166 

dipole radiation pattern, 136 

dispersion, 83, 89 

dispersion relation, 86 

dispersion relation, conductor, 94 

dispersion relation, dilute plasmas, 100 

dispersion relation, waveguide, 124 

dispersion, anomalous, 92 

dispersion, normal, 92 

dispersive media, 82 

displacement current, 21 

Doppler broadening, 52 

electric dipole radiation, 147 
electric field, induced, 14 
electro-motive force, 11 
electromagnetic waves, 44, 82 
electromagnetic waves, energy, 49 
electromagnetic waves, momentum, 49 
emf, motional, 12 
energy, 18 

energy in electromagnetic waves, 49 
energy, conservation, 23 
energy, magnetic field, 18 
evanescent wave, 73, 74 

far zone, 145 
Faraday rotation, 104 
Faraday’s law, 13 

Faraday’s law, vector potential form, 14 
Fourier transform, 89 
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Fourier transform of gaussian, 91 
frequency, 35 
frequency, angular, 35 
frequency, cut-off, 102 
frequency, plasma, 100 
Fresnel equations, 61 

gaussian wave packet, 91 

grazing incidence, 72 

group velocity, 89, 91 

group velocity in dilute plasmas, 101 

group velocity, waveguide, 124 

induced electric field, 14 
inductance, 17 
inductance, mutual, 15 
inductor, 18 

inhomogeneous wave equations for the poten¬ 
tials, 143 

integral theorems, 171 
intensity, wave, 49 

Kirchhoff’s circuit rules, 38 

Larmor’s radiation formula, 132, 135 
laser, 53 

law of reflection, 59 
Lenz’ law, 13 
line width, natural, 51 
Lorentz gauge, 143 

Lorentz oscillator model of dielectrics, 83 
Lorentzian function, 87 
lossless transmission line equations, 40 
lossless transmission line, standing waves, 43 
luminiferous aether, 45 

magnetic dipole radiation, 149 
magnetic quadrupole radiation, 151 
maser, 53 

Maxwell stress tensor, 26 


Maxwell stress tensor, derivation, 167 

Maxwell’s equations, 20, 23, 175 

momentum density, 26 

momentum in electromagnetic waves, 49 

monochromatic plane wave solutions, 47 

monochromatic wave, 36 

motional emf, 12 

multipole expansion, vector potential, 146 
mutual inductance, 15 

natural line width, 51 
near zone, 145 
normal dispersion, 92 

one-dimensional wave equation, 34 
oscillator strengths, 88 

phase, 35 

phase constant, 35 
phase velocity, 37 

phase velocity in dilute plasmas, 101 

phase velocity, linear medium, 46 

phase velocity, lossless transmission line, 41 

phase velocity, waveguide, 124 

plane of incidence, 59 

plane wave, 36 

plasma frequency, 100 

plasmas, dilute , 97 

polarisation, 54 

polarisation, circular, 56 

polarisation, elliptical, 56 

polarisation, linear, 55 

Poynting vector, 25, 49 

Poynting vector, waveguide, 127 

Poynting’s theorem, 24 

pressure broadening, 52 

pressure, radiation, 49 

product rules, 172 

pulsar, dispersion measure, 104 
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pulsar, dispersion of pulses from, 103 
pulsar, rotation measure, 108 

quasi-monochromatic, 51 
quasi-static limit, 145 

radiation pressure, isotropic, 49 
radiation resistance, 139, 157 
radiation zone, 145 
radiation, electric dipole, 147 
radiation, magnetic dipole, 149 
radiation, magnetic quadrupole, 151 
rays, 37 
reflectance, 61 

reflectance of good conductors, 96 
refractive index, 58 
refractive index of conductor, 94 
retarded potentials, 145 
retarded time, 145 
rotation measure, 108 

scattering, 140 

scattering, by free electrons, 141 

scattering, by molecules, 142 

second derivatives, 172 

self-inductance, 17 

SI units and dimensions, 166 

skin depth, 95 

Snell’s Law, 60 

speed of light, 46 

spherical coordinates, 173 

spherical waves, 37 

standing waves, lossless transmission line, 43 
TECU, 104 

telegraphers’ equations, 40 
TEM wave, 115 
Thomson cross section, 142 
total electron content, vertical, 104 
total internal reflection, 68, 73 


transformer, 16 
transmission line, 38 

transmission line, characteristic impedance, 41 
transmission line, two-wire, 42 
transmittance, 61 

transverse electromagnetic wave, 115 
two-wire transmission line, 42 

units, SI, and dimensions, 166 
universal flux law, 13 

velocity, group, 89 
velocity, phase, 37 
visibility, 50 
VTEC, 104 

wave equation of conductor, 94 
wave equation, linear medium, 46 
wave equation, dilute plasma, 100 
wave equation, lossless transmission line, 41 
wave equation, one-dimensional, 34 
wave equation, vacuum, 46 
wave equations for the potentials, inhomoge¬ 
neous, 143 
wave intensity, 49 
wave number, 35 
wave packet, 89 
wave packet, gaussian, 91 
wave propagation in dilute plasmas, 97 
wave vector, 37 
wave, complex form, 36 
wave, monochromatic, 36 
wave, monochromatic plane wave, 47 
wave, plane, 36 
wave, TEM, 115 

wave, transverse electromagnetic, 115 
waveguide modes, 119 
waveguide TE modes, 122 
waveguide TM modes, 121 
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waveguide, dispersion relation, 124 

waveguide, group velocity, 124 

waveguide, phase velocity, 124 

waveguide, Poynting vector, 127 

waveguide, TE wave, 116 

waveguide, TM wave, 116 

waveguide, transverse electric wave (TE), 116 

waveguide, transverse magnetic wave (TM), 116 

waveguides, 114 

wavelength, 34 

waves, electromagnetic, 44 

waves, spherical, 37 
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